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CHAPTER  1 


Summary 


A  series  of  diagnostic  experiments  have  been  compiled  to  investigate  rf  and  microwave 
discharges  in  N2,  such  as  those  used  for  plasma  processing  of  electronic  materials.  To 
develop  a  picture,  as  complete  as  possible,  of  the  discharge  chemistry,  several  reactive 
intermediates,  e.g.  N,  NaCA),  N2(B)  and  N2(C),  were  monitored.  An  alternative  two- 
photon  laser  induced  fluorescence  (TALIF)  scheme  for  detection  of  atomic  nitrogen 
has  been  characterized  and  used  to  monitor  and  evaluate  various  atomic  nitrogen 
sources.  Two-photon  excitation  at  207  nm  to  the  N(3p)  ^83/3  state  was  demonstrated 
to  be  superior  to  the  traditional  211  nm  excitation  as  a  TALIF  scheme  to  monitor 
the  ground  N  atom  state.  Most  striking  was  the  low  quenching  rate  of  the  upper  (3p) 
^S|/2  state  by  N2  compared  to  the  traditional  technique.  The  two-photon  excitation 
rate  at  207  nm  was  also  measured  to  be  a  factor  of  3.5  greater  than  the  traditional 
scheme.  The  results  suggest  that  the  207  nm  excitation  scheme  should  be  employed 
in  future  N  atom  diagnostic  work. 

Photodissociation  of  N2O  and  subsequent  detection  of  N  atom  photofragments  by 
TALIF  was  also  performed  using  a  single  207  nm  laser  pulse.  Measurements  of  the 
TALIF  upper  state  quenching  by  N2O  indicated  that  the  “hot”  N  atom  photofrag¬ 
ments  had  an  average  kinetic  energy  much  greater  than  that  of  N  atoms  produced  in 
a  microwave  discharge. 

Prom  measurements  within  the  pulsed  rf  discharge,  a  kinetic  model  was  developed 
which  predicted  absolute  concentrations  of  various  discharge  species.  The  behavior 
of  several  discharge  species  were  also  monitored  in  the  N2  microwave  discharge  and 
compared  to  the  model  for  the  addition  of  small  amounts  of  H2  and  NH3. 
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Atomic  nitrogen  destruction  by  surface  and  volume  reaction  was  investigated  by 
temporal  N  atom  decay  within  a  pulsed  rf  discharge.  Among  the  surface  materials 
studied,  stainless  steel  displayed  the  highest  N  atom  surface  loss  rate  while  boron 
nitride  proved  to  be  the  most  resistant  to  surface  recombination.  The  volume  loss 
rate  of  atomic  nitrogen  was  also  quantified  to  be  substantially  increased  with  the 
addition  of  as  little  as  0.1%  O2  in  an  N2  rf  discharge. 
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CHAPTER  2 


Introduction 


2.1  Atomic  nitrogen  in  a  glow  discharge  and  its  application 

The  study  of  reactive  species  in  glow  discharges  can  produce  physical  and  chemi¬ 
cal  information  that  is  used  to  optimize  plasma,  processing  techniques  for  electronic 
materials.  An  absolute  understanding  of  the  interactions  within  a  particular  plasma 
reactor  would  be  ideal,  but  is  most  often  elusive  due  to  the  non-equilibrium  nature 
of  the  low  pressure  glow  discharge.  The  discharge  gives  rise  to  the  existence  of  highly 
energetic  electrons  that,  through  coUisional  processes,  can  potentially  generate  any 
species  that  may  arise  from  dissociation  and  recombination  of  the  feed  gas  mixture. 
Although  purely  empirical  optimization  methods  for  plasma  deposition/etching  pro¬ 
cesses  have  been  used  in  the  past,  the  knowledge  of  gas  phase  chemistry  as  a  function 
of  process  parameters  would  be  invaluable  to  the  researcher  charged  with  developing 
an  optimized  recipe. 

The  need  for  atomic  nitrogen  plasma  sources  for  advanced  electronic  material 
processing  has  rekindled  an  interest  in  the  nitrogen  glow  discharge  in  recent  years. 
In  most  glow  discharges  through  Nj,  the  reactive  species  with  the  greatest  num¬ 
ber  density  is  atomic  nitrogen  in  the  ground  (2p®)  ^83/3  state.  Furthermore,  the 
atomic  concentration  decays  slowly  if  a  low  surface  recombination  rate  is  maintained. 
As  atomic  nitrogen  is  an  essential  ingredient  in  plasma  processing  recipes,  such  as 
deposition  of  Si3N4  dielectric  material[l],  epitaxial  growth  of  semiconducting  III-N 
compounds[2,  3],  doping  of  ZnSe[4],  and  nitride  hardening  of  steel  surfaces,[5]  an  un¬ 
derstanding  of  the  production  and  decay  of  N  atoms  is  important  to  the  advancement 
of  plasma  processes. 
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2.2  Detection  of  atomic  nitrogen  produced  in  a  discharge 

A  first  step  in  developing  a  comprehensive  model  of  atomic  nitrogen  production  and 
subsequent  interactions  is  to  obtain  a  mapping  of  the  N  atom  concentration  produced 
within  the  discharge  under  various  conditions.  In  the  1980s,  laser  based  techniques 
were  developed  to  measure  concentrations  of  several  light  atomic  species,  which  com¬ 
bined  laser  Raman-shifting  with  two-photon  absorption  laser  induced  fluorescence 
(TALIF)[6].  More  recently,  advanced  non-linear  optical  crystals,  such  as  the  BBO, 
have  replaced  the  less  efficient  long  path  Raman  cells  to  generate  the  required  ultra¬ 
violet  laser  light,  down  to  nearly  200  nm,  for  atomic  TALIF  methods[7].  The  TALIF 
technique  has  been  demonstrated  for  several  small  atoms[6,  8,  9],  but  can  present  a 
challenge  in  achieving  good  detection  sensitivity  since,  in  general,  a  two-photon  exci¬ 
tation  cross  section  will  be  vastly  smaller  than  that  of  a  typical  one  photon  process. 

A  TALIF  scheme  for  atomic  nitrogen,  first  demonstrated  by  Bischel  et  ol.[6],  re¬ 
lying  on  absorption  from  the  ground  state  to  the  (3p)  ^D°  state  is  well  known.  The 
(3p)  ^D°  state  is  the  lowest  energy  state  in  which  two-photon  laser  absorption  is 
allowed  from  the  ground  state.  The  strongest  two-photon  transition  from  the  groimd 
(2p®)  ^83^2  state  to  the  (3p)  multiplet  has  been  found  to  involve  the  J=7/2 
level[10]  at  94883.1  cm“^  above  the  ground  state.  This  corresponds  to  the  absorption 
of  two  laser  photons  with  wavelengths  of  211  nm.  Once  the  (3p)  state  has  been 
populated  by  two-photon  absorption,  the  observed  fluorescence  line  results  from  the 
radiative  transition  (3p)  (3s)  ^P5/2  with  a  wavelength  of  868  nm.  This  is  the 

TALIF  sdieme  that  is  prevalent  in  the  literature  and  has  been  used  in  the  analysis 
of  flame  combustion [7],  low-pressure  post  discharge  systems[5],  and  in  a  free-burning 
arc  discharge[ll]. 

Notwithstanding  the  popularity  of  this  scheme,  it  has  several  limitations.  The 
fluorescence  emission  wavelength  of  868  nm  is  so  red,  that  many  common  photomul¬ 
tiplier  detectors  have  poor  sensitivity  here.  Additionally  the  emission  coincides  with 
strong  molecular  N2  emission  that  is  often  present  in  systems  for  which  one  would 
like  to  use  the  TALIF  approach  to  measure  N  atom  concentration.  This  coincidence 
lowers  significantly  the  sensitivity  of  the  TALIF  diagnostic.  Finally  the  ^D°  state  is 
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known  to  have  a  fairly  large  quenching  rate.  This  both  limits  the  sensitivity  of  the  di¬ 
agnostic  and  ultimately  the  accuracy  of  N  atom  concentration  measurements  because 
of  uncertainties  in  the  relatively  large  corrections  that  must  be  made  to  account  for 
quendiing. 

In  Chapter  4,  the  performance  of  an  alternative  atomic  nitrogen  TALIF  scheme 
with  laser  excitation  at  207  nm  is  compared  to  the  well  known  TALIF  scheme  with 
211  nm  excitation  described  above.  Comparisons  involving  several  parameters  were 
done  to  determine  the  merit  of  each  TALIF  sdieme  under  a  variety  of  discharge 
conditions.  The  alternative  TALIF  technique  demonstrated  marked  improvements  in 
the  form  of  a  higher  two-photon  absorption  rate,  a  lower  quenching  rate  of  the  TALIF 
upper  state,  and  a  more  practical  range  of  fluorescence  detection  wavelength  compared 
to  the  traditional  TALIF  scheme.  The  alternative  TALIF  scheme  also  demonstrated 
at  least  equivalent  merit  in  comparisons  of  the  generation  of  laser  power  and  the 
background  emission  intensity  in  the  spectral  region  of  detection.  An  excellent  case  is 
made  for  the  replacement  of  the  traditional  atomic  nitrogen  TALIF  sdieme  of  choice 
with  this  alternative  scheme. 

2.3  Photodissociation  and  TALIF  detection  of  N  atoms  from 
N2O. 

In  the  study  of  atomic  nitrogen,  there  may  be  situations  where  it  is  useful  to  conduct 
TALIF  measurements  on  N  atoms  produced  by  a  source  other  than  a  discharge.  A 
discharge  source  can  be  a  very  eflicient  method  of  N  atom  production,  but  the  various 
radical  species  accompanying  the  discharge  source  as  well  as  the  excessive  background 
radiation  can  be  prohibitive  for  certain  types  of  measurements.  Production  of  atomic 
nitrogen  by  laser  photodissociation  of  a  nitrogen  containing  molecule  could  become 
a  powerful  experimental  tool  by  eliminating  these  radical  species  and  background 
radiation  concerns.  The  photolysis  technique  would  be  especially  useful  if  a  molecule 
was  used  with  a  dissociation  energy  that  matches  the  N  atom  TALIF  laser  photon 
energy.  The  molecule  could  thereby  be  photodissociated  and  then  detected  by  TALIF 
within  the  same  laser  pulse. 
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Both  H  atoms  and  O  atoms  have  previously  been  produced  and  probed  by  TALIF 
in  this  manner.  The  H  and  0  atoms  were  produced  from  photolysis  of  C2H2  and  NO2 
with  laser  wavelengths  of  205  nm  and  226  nm  respectively  and  probed  by  TALIF  with 
the  same  laser  pulse[12].  Such  photodissociation/TALIF  tedmiques  have  been  used 
to  calibrate  the  TALIF  signal  collection  view  factor  in  a  reactor  cell  that  has  spatial 
regions  of  restricted  view.  The  reactor  cell  can  be  uniformly  filled  with  the  precursor 
molecule  which  results  in  a  spatially  uniform  atomic  density  and  corresponding  TALIF 
intensity  except  where  the  view  is  restricted.  The  data  from  a  spatial  scan  of  the 
photodissociation  induced  TALIF  signal  in  a  reactor  cell  can  be  used  to  determine 
the  cell  view  factor  for  use  with  other  atomic  sources,  such  as  a  discharge. 

Bengtsson  et  oZ.[13]  introduced  nitrous  oxide,  N2O,  as  a  suitable  precursor  for 
producing  ground  state  N  atoms  by  photodissociation  with  subsequent  TALIF  ab¬ 
sorption  to  the  (3p)  or  (3p)  state.  Their  data  indicated,  though,  that 

severe  quenching  of  the  upper  TALIF  state  occurs  in  less  than  0.5  Torr  of  N2O.  This 
quenching  would  limit  the  usefulness  of  the  technique,  as  the  benefits  of  raising  the 
N2O  pressure  to  increase  the  N  atom  yield  would  be  negated. 

In  this  work,  N  atom  collisional  partners  responsible  for  the  quenching  during  the 
nitrous  oxide  photolysis/N  atom  TALIF  event  were  identified  and  investigated.  The 
quenching  data  in  the  photolysis  region  was  compared  to  data  from  an  experiment 
in  which  candidate  quenching  species  were  isolated  and  mixed  with  atomic  nitro¬ 
gen  firom  a  flowing  microwave  discharge.  Fluorescence  decay  measurements  from  a 
TALIF  probe  determined  the  quenching  due  to  the  additive.  It  was  found  that  the 
N2O  molecule  itself  quenches  the  upper  TALIF  state  significantly  more  than  does  N2. 
The  quenching  rate  of  the  isolated  N2O  was  still  significantly  less  than  the  measured 
quenching  in  the  photol5^is  eiqieriment.  This  discrepancy  is  explained  by  accounting 
for  the  excess  photolysis  energy  that  is  partitioned  among  the  photofragments,  en¬ 
hancing  N  atom  kinetic  energy.  This  elevates  the  velocity  of  the  atomic  nitrogen  in 
the  photolysis  region  far  above  that  of  the  thermalized  atoms  in  the  flowing  afterglow 
of  the  microwave  discharge.  The  higher  velocity  during  photodissociation  results  in 
a  greater  collision  frequency  between  N  and  N2O  and  hence  the  higher  quenching 
rate.  Another  candidate  molecule  that  may  contribute  to  the  N  atom  upper  TALIF 
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state  quenchiiig  is  NO.  The  NO  density,  though,  was  estimated  to  be  small  in  the 
photolysis  region  and  its  contribution  to  the  quenching  negligible.  In  Chapter  5,  the 
photolysis  of  N2O  and  subsequent  N  atom  TALIF  detection  is  found  to  be  achiev¬ 
able  with  relatively  low  atomic  yield  at  low  pressures.  Substantial  quendiing  of  the 
TALIF  signal  by  N2O  eliminates  much  of  the  N  atom  yield  advantages  at  higher  N2O 
pressures. 

2.4  Behavior  of  atomic  nitrogen  and  excited  molecular  species 
produced  in  N2  and  N2/H2  high  frequency  discharges 

Much  research  has  been  performed  on  kinetic  modeling  of  discharges  in  both  pure 
N2  [14]  and  pure  H2  [15].  When  applying  these  pure  gas  discharges  as  atom  sources, 
the  fractional  dissociation  in  each  has  tended  to  be  low.  For  example,  a  standard  mi¬ 
crowave  discharge  with  pure  H2  has  been  reported  to  yield  <10%  dissociation[16, 17] 
while  microwave  nitriding  reactors  with  pure  N2  reportedly  yield  <  1%  dissociation[18, 
19].  Recent  studies  of  discharges  in  N2/H2  gas  mixtures,  though,  have  shown  that 
the  degree  of  dissociation  can  vary  significantly  from  the  pure  gas  cases.  In  DC, 
RF,  and  microwave  N2/H2  discharges,  small  percent  concentrations  of  H2  has  been 
found  to  result  in  near  complete  dissociation  of  the  H2  leading  to  an  efficient  source 
of  atomic  hydrogen[17].  It  has  been  theorized  that  the  heterogeneous  heavy  particle 
mllisinn  involving  metastable  N2(a^  additive  H2  is  responsible  for  the 

high  fractional  dissociation  of  hydrogen  in  this  gas  mixture  [20].  For  example,  in  a  DC 
discharge  with  [H2]/([H2]-I-[N2])  <  10%,  this  dissociative  path  involving  N2(a' 
has  been  credited  with  ~80-90  %  of  the  atomic  hydrogen  production  [21].  As  the 
percent  of  H2  increases,  the  relative  contribution  of  this  nitrogen  metastable  to  the 
H  atom  production  decreases  due  to  quenching  of  the  N2(a'  ^E~)  by  molecular  and 
atomic  hydrogen[20]. 

The  effects  of  N2/H2  gas  mixtures  on  atomic  nitrogen  production  have  been  stud¬ 
ied  to  some  extent,  but  have  received  much  less  attention  than  the  effects  on  atomic 
hydrogen  production.  It  was  noted  in  early  gas  mixture  work  that  small  additions  of 
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H2  to  an  N2  discharge  caused  afterglow  emission  from  N  atom  recombination  to  inten¬ 
sify,  indicating  an  increase  in  atomic  nitrogen  concentration  [22].  This  phenomenon 
was  attributed  to  a  “poisoning”  of  the  walls  by  hydrogen  that  prohibited  wall  recom¬ 
bination  of  N  atoms  and  therefore  an  increase  in  N  atom  concentration.  More  recently 
in  flowing  N2/H2  DC  discharges,  production  of  atomic  nitrogen  has  been  reported  to 
increase  up  to  a  factor  of  3,  as  the  %  H2  added  is  varied  from  0  to  0.5%[23].  This 
large  increase  in  N  atom  production  occurs  [20]  mainly  because  the  H2  additive  alters 
the  DC  discharge  impedance  and  results  in  more  power  deposited  in  the  discharge, 
operating  at  constant  current.  The  increase  in  the  DC  discharge  impedance  as  H2  is 
introduced  has  been  found  to  occur  due  to  the  depletion  of  the  dominant  associative 
ionization  process  involving  collisions  of  N2(a'  and  N2(A  metastables,  since 

both  metastables  are  readily  quenched  by  molecular  and  atomic  hydrogen.  The  DC 
power  supply  responds  to  the  increased  impedance  by  increasing  the  voltage  applied 
across  the  DC  discharge  to  sustain  the  same  discharge  current,  and  therefore  signif¬ 
icantly  more  power  is  deposited  in  the  discharge.  It  has  been  suggested  as  well  that 
the  addition  of  H2  in  the  discharge  has  the  coupled  effect  of  “heating  up”  the  electron 
energy  distribution  function  (EEDF)[20,  23,  24]  so  to  enhance  the  high  energy  tail  of 
the  EEDF  where  electrons  have  sufficient  energy  to  create  direct  impact  dissociation. 
The  mean  electron  energy  experiences  the  greatest  increase  as  the  %  H2  is  varied  from 
0-10  %,  since  the  electric  field  increases  significantly  in  this  range  as  discussed  above. 
As  the  N2/H2  disdiarge  becomes  further  diluted  with  H2,  the  mean  electron  energy 
continues  to  slowly  rise  due  to  the  declining  influence  of  the  large  inelastic  losses  of 
N2[20]. 

It  is  of  interest  to  determine  whether  a  more  applicable  nitrogen  atom  source, 
such  as  a  high  frequency  rf  or  microwave  discharge,  would  benefit  from  a  similar 
addition  of  H2  to  the  N2  feed  gas  while  operating  at  a  constant  power  deposition 
to  the  discharge.  Ricard  reported  variations  in  the  N  atom  density  in  the  flowing 
after^ow  of  a  microwave  nitriding  reactor  as  H2  was  introduced  into  an  Ar-N2  feed 
gas  mixture  at  44  Torr[5].  Ricard  measured  a  ~10%  peak  increase  in  N  atom  density 
with  the  addition  of  0.1%  H2  in  the  gas  mix.  With  further  H2  addition,  there  was  then 
a  sharp  decrease  in  atomic  nitrogen,  with  the  N  atom  density  at  1%  H2  being  80% 
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of  that  in  the  pure  Ar-N2  discharge.  Ricard  attributed  the  initial  increase  in  atomic 
nitrogen  to  the  depletion  of  the  N2(a')-  N2(A)  associative  ionization  processes  and 
subsequent  increased  E/N,  as  discussed  above,  which  enhanced  the  electron  impact 
dissociation  of  N2.  No  discussion  was  offered,  though,  on  the  sharp  decrease  in  N 
atom  density  at  higher  concentrations  of  H2. 

Chapter  6  presents  a  detailed  investigation  of  the  heavy  particle  interaction  of 
nitrogen  and  hydrogen  species  in  high  frequency  discharges  and  the  resulting  effects 
on  the  density  of  atomic  nitrogen.  Atomic  nitrogen  TALIF  was  utilized  along  with  LIE 
ofN2(A3E+)andN2(B  g)  species  as  well  as  emission  spectroscopy  of  electronically 
excited  molecular  stafes  and  to  determine  the  behavior  of  several  important  species 
produced  within  eadi  disdiarge.  A  low  power  parallel  plate  rf  discharge  cell  was 
analyzed  to  obtain  details  of  the  state  specific  heavy  particle  interaction  during  short 
periods  after  the  discharge  has  turned  off.  The  coupled  behavior  of  atomic  nitrogen, 
N2(X^E+,v>5),  N2(A  ®E+)  and  N2(B  ^Ilg)  in  the  pure  nitrogen  post-discharge  was 
clearly  observed.  The  addition  of  small  amounts  of  H2  to  the  rf  discharge  led  to 
measurements  of  vastly  increased  quenching  rates  of  N2(A  ),  whidi  was  attributed 
to  the  formation  of  atomic  hydrogen. 

A  higher  power  flowing  microwave  discharge  system  was  also  analyzed  with  N2/H2 
gas  mixtures.  Diagnostics  on  this  system  were  limited  to  downstream  N  atom  TALIF 
and  optical  emission  spectroscopy  of  the  discharge  volume  and  late  afterglow.  The 
data  from  the  microwave  system  demonstrated  enhancement  in  atomic  nitrogen  pro¬ 
duction  for  specific  N2/H2  gas  mixtures.  Analysis  of  this  data,  along  with  the  pulsed  rf 
discharge  data,  led  to  the  development  of  a  model  that  predicts  an  increase  in  atomic 
nitrogen  production  due  to  enhanced  mean  electron  energy  when  small  amounts  of 
H2  are  added  to  an  N2  discharge.  The  model  presented  in  Chapter  6  extends  the 
theory  to  attribute  the  observed  decline  in  N  atom  density  at  higher  %H2  to  the  de¬ 
population  of  N2(X^S+,v')  states,  where  v'  designates  a  highly  vibrationally  excited 
state.  This  depopulation  of  the  N2(X^S+,v')  states  is  theorized  to  result  in  a  decrease 
in  N2  dissociation  via  the  collisional  channels  and  electron  impact  channels  involving 
N2(X^E+,v').  The  combination  of  data  from  the  different  diagnostics  supports  this 
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theory  in  finding  this  N  atom  decline  at  higher  %B.2  more  prominent  at  higher  pres¬ 
sures,  where  the  dissociation  channels  involving  N2(X^S+,v')  would  be  contributing 
more  to  the  total  N2  dissociation. 

Finally,  a  gas  mixture  e3q)eriment  was  performed  in  the  flowing  microwave  dis¬ 
charge  where  ammonia,  NH3,  was  substituted  for  H2.  As  with  the  N2/H2  mixtures, 
small  amounts  up  to  10%  of  NH3  were  introduced  into  the  nitrogen  feed  gas.  The 
enhancement  in  N  atom  production  for  the  N2/NH3  gas  mixtures  were  remarkably 
similar  to  the  N2/H2  mixtures  for  all  conditions  tested.  This  similarity  is  an  indi¬ 
cation  that  both  H2  and  NH3,  when  dilute  in  N2,  experience  very  high  fractional 
dissociation  during  the  ms  that  the  flowing  gas  is  resident  within  the  discharge 
volume.  This  high  dissociation  rate  leaves  a  similar  chemistry  of  fundamental  species 
in  the  post  discharge  of  either  gas  mixture  and  hence  the  similar  effect  on  atomic 
nitrogen  concentration. 

2.5  Surface  and  volume  loss  of  atomic  nitrogen  in  a  discharge 
reactor 

The  net  loss  of  atomic  nitrogen  in  a  low-pressure  pure  N2  discharge  reactor  is  rela¬ 
tively  slow  when  compared  to  many  other  discharge  produced  radicals.  In  a  pure  N2 
discharge  void  of  impurities,  the  gas  phase  losses  axe  known  to  be  very  slow,  primarily 
through  three-body  recombination  with  N2[22].  The  surface  loss  rate  of  atomic  nitro¬ 
gen  in  pure  N2  on  pyrex  and  certain  metal  surfaces  is  also  considered  to  be  relatively 
lotv[25].  In  contrast,  the  introduction  of  small  amounts  of  certain  gas  phase  impurities 
can  cause  a  significant  increase  in  the  decay  rate  of  atomic  nitrogen.  The  occurrence 
of  unwanted  gas  impurities  is  most  often  a  result  of  a  virtual  leak  from  material  in  the 
reactor  system  or  an  air  leak,  resulting  in  O2  or  water  vapor  in  the  reactor.  A  minute 
amount  of  these  oxygen  containing  impurities  in  the  N2  discharge  can  lead  to  the 
formation  of  reactive  species,  such  as  NO,  or  increase  the  N  atom  surface  reactivity, 
both  of  which  can  be  destructive  to  the  atomic  nitrogen  concentration.  With  enough 
O2  present  in  the  system,  these  impurity  induced  N  atom  depletion  processes  could 
become  dominant  compared  to  other  loss  channels. 
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The  low  N  atom  volume  loss  rate  in  a  pure  N2  discharge  often  leaves  surface 
recombination  as  the  main  loss  channel[14],  depending  on  the  reactor  configuration. 
The  surface  reactivity  of  atomic  nitrogen,  therefore,  becomes  important,  not  only 
where  a  reaction  is  occurring  on  a  wafer  surface,  but  also  as  the  N  atoms  interact  with 
other  smfaces  within  the  discharge  cell.  This  collection  of  surface  losses  can  constitute 
the  primary  atomic  loss  channel  and  have  a  significant  influence  in  controlling  the  N 
atom  flux.  It  is  well  known  that  creating  a  high  flux  N  atom  source  is  complicated 
by  the  difficulty  in  obtaining  a  high  degree  of  dissociation  firom  the  N2  feed  gas.  This 
low  fractional  dissociation  is  due  to  the  large  bond  strength  of  N2  and  the  dbanneling 
of  large  amounts  of  input  energy  into  the  Vibrational  manifold  of  N2(X^S+,v).  This 
limitation  in  the  N  atom  flux  of  a  low  pressure  discharge  source  underscores  the  need 
to  maintain  reactor  surfaces  with  reduced  N  atom  re-association.  An  understanding  of 
atomic  nitrogen  recombination  and  reaction  rate  on  various  materials  would  therefore 
be  of  great  advantage  in  the  design  of  a  reactor  employing  an  N  atom  source. 

In  Chapter  7,  the  atomic  nitrogen  loss  due  to  both  volume  and  surface  reaction 
in  a  parallel  plate  rf  reactor  was  investigated  using  a  pulse  N2  discharge  and  TALIF 
techniques.  The  surface  loss  of  atomic  nitrogen  was  analyzed  by  determining  re¬ 
combination  rates  of  atomic  nitrogen  on  various  surfaces  within  the  parallel  plate  rf 
discharge  reactor.  Metal  electrode  surfaces  as  well  as  various  substrate  materials  and 
other  reactor  surface  materials  were  investigated  for  their  reactivity  with  atomic  ni¬ 
trogen  within  a  discharge  environment.  The  loss  rate  of  N  atoms  was  also  quantified 
due  to  minute  O2  impurities  introduced  into  an  parallel  plate  pulsed  rf  N2  discharge. 
The  introduction  of  O2  simulates  the  unwanted,  but  often  encountered,  effect  of  small 
amounts  of  air  or  water  vapor  infiltrating  the  reactor.  A  kinetic  model  was  developed 
for  the  N2/O2  discharge  and  solved  numerically  for  conditions  that  matdbed  the  ex¬ 
periment.  The  experimental  depletion  of  N  atom  density  with  the  addition  of  0.2% 
O2  to  the  N2  rf  discharge  was  accurately  depicted  by  the  kinetic  model. 
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CHAPTER  3 


Experimental 


3.1  Laser  light  generation  and  signal  detection 

Several  laser  configurations  were  utilized  in  this  study  for  detection  of  various  species 
produced  in  gas  discharges.  The  primary  laser  technique  was  atomic  nitrogen  TALIF 
with  laser  wavelength  at  207  nm,  which  was  applied  in  experiments  in  all  of  the 
following  chapters.  In  Chapter  4,  experiments  with  a  second,  but  more  well  known, 
atomic  nitrogen  TALIF  scheme  with  laser  absorption  at  211  nm  were  done  as  a 
comparison  to  the  207  nm  sdieme.  In  the  investigation  in  Chapter  5,  ultraviolet  laser 
photodissociation  of  N2O  was  coupled  with  atomic  nitrogen  TALIF.  Single  photon 
LIF  techniques  were  also  employed  in  Chapter  6  for  the  detection  of  NaCA  ®S+)  and 
N2(B  ^IIj)  produced  in  Na  and  N2/H2  discharges.  The  experimental  details  of  these 
laser  techniques  are  described  in  this  section. 

For  the  two  atomic  nitrogen  TALIF  techniques,  laser  excitation  was  either  at  207 
nm  or  211  nm  with  the  fluorescence  detected  at  747  nm  or  870  nm  respectively.  An 
energy  level  diagram  of  the  N  atom  TALIF  schemes  along  with  an  experimental  laser 
excitation  spectrum  for  each  scheme  is  shown  in  Fig.  1.  Tunable  radiation  to  provide 
the  207/211  nm  laser  pulse  energy  for  the  two  N  atom  TALIF  schemes  was  generated 
by  a  Nd-YAG  pumped  dye  laser  system  plus  two  non-linear  optic  crystal  frequency 
conversions.  The  second  harmonic  output  of  the  Nd-YAG  laser  pumped  LDS-698  dye 
to  produce  a  676/698  nm  dye  laser  output.  This  dye  output  was  then  mixed  with 
the  residual  1064  nm  Nd-YAG  fundamental  beam,  resulting  in  the  production  of  a 
414/422  nm  pulse  which  was  then  polarization  rotated  by  90°.  The  414/422  nm  beam 
was  then  passed  through  a  BBO  crystal  for  frequency  doubling  to  produce  the  desired 
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207/211  nm  beam.  A  diagram  of  the  experimental  laser  configuration  and  fiuorescence 
detection  for  the  two  TALIF  schemes  along  with  the  associated  electronics  is  shown 
in  Fig.  2. 

This  laser  configuration  was  chosen  above  a  more  popular  configuration  [7, 13]  for 
211  nm  generation  due  to  the  ease  and  repeatability  of  alternating  between  211  nm 
and  207  nm  output.  The  violet  414/422  nm  beam  in  this  configuration,  however,  has 
a  relatively  wide  linewidth  on  the  order  of  1  cm“^,  since  it  is  created  through  mixing 
with  the  Nd-YAG  fundamental.  This  results  in  a  relatively  low  frequency  doubling 
efficiency  (~5%)  in  the  BBO  and  a  '^^l  cm“^  linewidth  207/211  nm  beam  that  gives 
a  lower  two-photon  atomic  absorption  rate  per  laser  energy.  However,  the  resultant 
207/211  nm  output  of  150  /xJ  proved  to  be  ample  for  this  effort  and  the  alignment 
advantages  warranted  the  use  of  this  new  configuration.  Indeed  use  of  greater  powers 
would  raise  the  possibility  of  significant  ionization  from  the  excited  atomic  state.  A 
blue  dye  was  rejected  as  source  of  the  414/422  nm  beam  since  the  blue  dye  lifetimes 
are,  in  general,  much  shorter  than  LDS-698  and  a  standard  blue  dye  was  not  found 
that  efficiently  tuned  from  414  nm  to  422  nm. 

The  induced  fluorescence  from  the  discharge  cell  was  collected  with  a  single  7.5 
cm  focal  length  f/1  lens.  The  fluorescence  was  passed  through  a  10  nm  bandwidth 
interference  filter  centered  at  747/870  nm  and  imaged  onto  an  RCA-31034  or  a  Hama¬ 
matsu  R928  photomultiplier  tube.  In  cases  where  TALIF  measurements  were  made 
within  the  discharge  volume,  a  normally-off  gated  PMT  configuration  was  used  to 
Tninimi'zp  the  Saturating  effects  of  focusing  the  intense  discharge  emission  onto  the 
PMT.  The  PMT  would  be  gated  on  for  only  -^^l  /xs  during  the  period  of  fiuorescence 
and  remain  blind  for  the  remainder  of  the  10  Hz  duty  cycle.  A  Hamamatsu  gated 
socket  was  used  in  conjunction  with  the  R928  PMT  for  this  application.  The  time 
resolved  TALIF  signal  was  digitized  with  a  LeCroy  9450  oscilloscope  and  processed 
with  Windows  based  Lab  View  software. 

In  Chapter  6  the  laser  detection  efforts  turned  from  atomic  nitrogen  to  the  elec¬ 
tronically  excited  molecular  nitrogen  states  of  N2(A  and  N2(B  ^Hp).  In  these 
cases,  single  photon  laser  induced  fluorescence  techniques  were  employed.  For  prob¬ 
ing  the  N2(A  metastable  state,  an  LIF  technique  was  used  with  excitation  from 
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Figure  1:  (a)  Energy  level  diagram  of  N  atom  TALIF  schemes  with  excitation  at 
either  207  nm  or  211  nm  (b)  Experimental  laser  excitation  spectrum  for  each  TALIF 
scheme 
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414/422  nm  Discharge 


Figure  2:  Diagram  of  experimental  laser  configuration  and  fluorescence  detection  for 
atomic  nitrogen  TALIF  with  laser  excitation  at  either  207  nm  or  211  nm 
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the  N2(A  ®E+,v=0)  state  to  N2(B  ^np,v=3).  This  LIF  scheme  is  shown  as  the  lower 
transition  in  Fig.  3,  overlaid  with  the  potential  energy  curves  of  the  molecular  ni¬ 
trogen  states.  An  excitation  spectrum  for  this  scheme  is  shown  in  Fig.  4  (a),  where 
the  bandhead  with  peak  laser  absorption  is  seen  to  occur  at  687.5  nm.  The  second 
harmonic  output  of  the  Nd-YAG  laser  pumped  LDS-698  dye  to  produce  the  687.5  nm 
dye  laser  output.  The  fluorescence  from  this  scheme  for  A  v=2  peaks  at  763  nm  and 
is  observed  with  a  10  nm  bandwidth  interference  filter  centered  at  760  nm. 

The  LIF  probe  of  N2(B  ^Ilg)  involved  excitation  from  the  N2(B  ^n5,v=l)  state  to 
N2(C  ®n„v=0)  as  shown  in  the  upper  transition  in  Fig.  3.  The  laser  excitation  peak 
was  within  a  dominant  bandhead  at  357.6  nm  as  seen  in  the  excitation  spectrum  in 
Fig.  4  (b).  The  laser  pulse  train  involved  the  Nd-YAG  second  harmonic  output  pump- 
iug  LDS-698  dye  producing  a  715  nm  dye  laser  output  which  was  frequency  doubled 
to  generated  the  357.6  nm  excitation  pulse.  The  fluorescence  from  the  N2(C  *n„v=0) 
state  was  collected  through  a  10  nm  bandwidth  interference  filter  centered  at  337 
nm,  which  corresponds  to  the  Av=0  bandhead  of  the  N2(C-B)  transition.  Diagrams 
of  the  experimental  laser  configuration  and  fluorescence  detection  of  the  LIF  probes 
of  the  N2(A  ®S+,v=0)  and  N2(B  ®n5,v=l)  states  are  shown  in  Fig.  5.  It  should  be 
noted  that  the  laser  wavelengths  for  both  molecular  LIF  schemes  were  conveniently 
within  the  range  of  the  LDS-698  dye,  already  in  use  with  the  N  atom  TALIF  work. 

In  Chapters  4  and  6,  the  laser  diagnostic  experiments  were  augmented  with  ob¬ 
servations  of  emission  spectra  from  discharge  species.  The  discharge  emission  and 
afterglow  spectra  were  collected  with  a  Digikrom  240  monochrometer  combined  with 
the  Hamamatsu  R928  PMT.  The  spectral  response  of  the  system  was  calibrated  us¬ 
ing  a  tungsten  ribbon  standard  lamp.  In  several  cases  in  Chapter  6,  narrow  band 
interference  filters  were  used  with  the  R928  PMT  to  monitor  the  emission  intensity 
from  specific  molecular  bandheads  under  changing  discharge  conditions. 
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Figure  3:  Laser  induced  fluorescence  sdiemes  overlaid  on  N2  potential  energy  curves 
including  the  A  B  ^IIj,  and  C  states 
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Laser  Excitation  Wavelength  (nm) 
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Figure  4:  Laser  excitation  spectra  for  the  LIF  probes  of  (a)  N2(A  *E+,v=0)  and  (b) 
N2(B  ®np,v=l)  states 
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Figure  5:  Diagram  of  experimental  laser  configuration  and  fluorescence  detection  for 
LIF  of  (a)N2(A  v=0)  and  (b)  N2(B  ®n^,v=l)  states 
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3.2  Downstream  microwave  discharge  reactor 

In  Chapters  4,  5,  and  6,  experiments  were  conducted  with  a  continuous  microwave 
flowing  discharge  in  a  2.5  cm  diameter  quartz  tube  with  an  AsTex  S-1000  2.45  GHz 
microwave  excitation  source.  A  sketch  of  the  microwave  discharge  reactor  and  as¬ 
sociated  optical  diagnostics  are  shown  in  Fig.  6.  The  AsTex  microwave  source  was 
capable  of  applying  1  kW  of  power  to  the  discharge,  but  nominally  100-200  W  were 
used  in  the  experiments.  The  construction  of  the  Astex  unit  did  not  leave  access  for 
laser  diagnostics  in  the  discharge  volume,  but  species  could  be  probed  just  a  few  cen¬ 
timeters  downstream  from  the  discharge.  Suprasil  laser  entry  and  exit  windows  were 
mounted  to  the  quartz  tube  5  cm  downstream  from  the  discharge  region  along  with 
a  quartz  fluorescence  exit  window  mounted  perpendicular  to  the  others.  A  side-arm 
gas  inlet  was  included  so  that  gas  could  be  added  downstream  of  the  discharge  but 
before  the  laser  probe  region.  The  gases  used  through  the  side-arm  port  were  NO, 
which  was  added  to  titrate  the  atomic  nitrogen  from  the  discharge,  and  NjO,  which 
was  added  to  test  its  collisional  quenching  effect  on  the  excited  atomic  N  (3p)  ■*S|/2 
state.  Gas  flow  and  pressure  in  the  microwave  discharge  system  were  governed  by  a 
mass  flow  controller/  downstream  throttle  value  combination  which  were  controlled 
by  an  MKS-146  control  unit.  An  oil  free  molecular  drag  pump  and  diaphragm  pump 
in  series  evacuated  the  discharge  cell. 

3.3  Pulsed  parallel  plate  rf  reactor  cell 

The  rf  parallel  plate  discharge  system  used  in  Chapters  6  and  7  was  a  versatile  ex¬ 
perimental  tool  when  used  in  conjunction  with  laser  diagnostic  probes  described  in 
Section  3.1.  The  rf  pulsed  discharge  was  synchronized  with  the  laser  probe  to  study 
the  temporal  behavior  of  various  nitrogen  species  while  the  entire  discharge  cell  could 
be  translated  in  two  dimensions  to  study  the  spatial  distribution  of  the  species.  The  rf 
parallel  plate  plasma  cell,  shown  in  Fig.  7,  was  a  scaled  down  version  of  a  commercial 
plasma  processing  reactor  with  top/bottom  interchangeable  metal  disk  electrodes, 
5  cm  in  diameter  and  an  interelectrode  gap  of  2.2  cm.  The  bottom  electrode  was 
grounded  while  the  top  was  driven  by  a  10  MHz  rf  power  supply.  An  adjustable 
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N2^2  Gas  In 


Figure  6;  Sketch  of  AsTex  microwave  reactor 
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Figure  7:  Sketch  of  parallel  plate  pulsed  rf  reactor  cell 


matching  network  was  also  in-line  with  the  power  supply.  The  discharge  cell  was 
mounted  on  a  X-Z  translation  stage  that  allowed  laser  probing  of  the  axial  and  radial 
regions  of  the  interelectrode  gap.  The  cell  could  be  fitted  with  electrodes  made  from 
similar  or  different  metals.  The  discharge  could  be  run  with  both  electrode  surfaces 
exposed  to  the  discharge  or  a  thin  substrate  could  be  placed  on  the  bottom  electrode. 
The  gas  flow  to  the  cell  was  maintained  by  the  same  MKS-146  control  system  as 
described  in  Section  3.2. 

Temporal  analysis  of  the  evolution  and  decay  of  the  atomic  nitrogen  concentration 
was  accomplished  with  TALIF  measurements  synchronized  to  a  pulsed  rf  discharge 
with  a  10  Hz  repetition  rate.  The  Nd:YAG  laser  was  triggered  by  a  digital  delay 
generator,  whidi  sjmchronized  the  laser  with  the  gated  PMT  collecting  the  TALIF 
signal.  The  rf  discharge  on  and  off  times  could  be  set  anywhere  within  the  10  Hz 
cycle,  allowing  TALIF  measurements  to  be  done  during  or  after  the  discharge.  The 
discharge-off  and  laser-on  times  had  a  jitter  of  less  than  10  ns. 

3.4  TALIF  Signal  Analysis 

The  TALIF  signal,  once  recorded  with  the  digitizing  oscilloscope,  was  analyzed  by 
either  integrating  the  signal  intensity  over  time  or  observing  the  time  decay  of  the  up¬ 
per  atomic  state  of  the  TALIF  scheme.  The  amplitude  and  shape  of  the  TALIF  signal 
naturally  depend  on  the  atomic  ground  state  concentration  along  with  the  radiative 
lifetime  and  collisional  quenching  properties  of  the  upper  TALIF  state.  The  TALIF 
signal  waveform  is  influenced  as  well  by  experimental  factors  such  as  fluorescence  de¬ 
tection  efliciency  and  laser  pulse  shape  and  energy.  Since  information  about  atomic 
concentration  and  collisional  quenching  rates  was  extracted  from  the  TALIF  signal 
data,  it  was  important  to  account  for  the  dependence  of  the  signal  on  all  influential 
parameters.  A  theoretical  expression  for  the  temporally  resolved  TALIF  signal  is 
derived  here  which  includes  the  known  parameters  on  whidi  the  signal  may  depend. 
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3.4.1  TALIF  Upper-state  Population 

The  rate  equation  for  the  laser-excited  upper  state  population,  N*,  in  an  N2  discharge 
is  [10] 

dN*  1 

—  =  -Ar-(i,lJV,]  +  i),  (3.1) 

where  is  the  ground  state  atom  population,  is  the  quenching  rate  of  N*  by  N2, 
Tr  is  the  radiative  lifetime  of  N*,  and  W(t)  is  the  rate  of  two-photon  excitation  jfrom 
Ng  to  N*.  The  loss  terms  can  be  expressed  in  terms  of  the  fluorescence  quantum  yield, 
Q,  where 


(3.2) 


l  +  A:5[iV2]r/ 

The  normalized  temporal  profile  of  the  laser  pulse,  F(t),  determines  the  time  depen¬ 
dency  of  W(t)  as[26] 


(3.3) 


where  is  the  two-photon  absorption  rate  coefficient  in  cm^s,  E  is  the  laser  pulse 
energy  density  in  J/cm^,  and  hv  is  laser  photon  energy.  If  a  Gaussian  temporal  profile 
for  the  laser  pulse  is  assumed,  then 


F{t)  = 


,-tV2s 


(3.4) 


where  the  FWHM  of  the  laser  temporal  profile  is  2s^2ln{2).  Eq.  3.1  can  then  be 
rewritten  as 


dN*  J_  ^  1 

dt  Qtt  {hvY 


Eq.  3.5  is  in  the  form 


dN* 

dt 


-h  VN*  =  Z(t), 


whidi  has  the  general  solution  of 

=  J  Z{t)e!^<^dt  +  C 


or 


(hz/)2  27rs2 


(3.5) 

(3.6) 

(3.7) 

(3.8) 
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Now  using  the  general  formula  [27] 


where  a=  ^,b=  aiid  c  =  0  gives 


JV-  =  k  f-  -  5^)  +  c] 

Ay/w{huys  [  Vs  2QTrJ 


The  boundary  condition  of  lim  N*=0  is  met  when  C— 1,  giving 

•'  t-+-00 


N*  = 


40r(hi/)2s 


•*  _  jtll. 


^s^Iaq^t?  eif  I  — 


S  2QTr 


+  1 


(3.11) 


3.4.2  TALIF  integrated  intensity 


Assuming  the  detected  fluorescence  in  the  experiment  is  from  the  strongest  transition 
from  N*,  where  the  Einstein  coefficient  meets  the  condition,  A20  »  then  the  time 
integrated  fluorescence  signal,  I,  in  the  experiment  can  be  written  as 


D 

=  ^/  N*it)dt 

Tr  J—00 


(3.12) 


where  D  is  the  signal  generated  per  fluorescence  photon  and  includes  experimental 
S3rstem  factors.  Combining  Eq.  3.11  and  Eq.  3.12  gives 

j  =  e»V4QV2  s  \  ^  r  e-^/Qrr^^  (3.I3) 

4y/^{hl')^TrS  [./-e»  V^  ^QTt  )  •'-oo 

Substituting  x  =  t  —  gives 

/  =  q-AIaq^t?  [  r  erf  f-)  ,  (3.14) 

4\/^{hl/yTrS  U-oo  \S/  J-ao  J 

which  is  solved  using  the  general  formula  [27] 

I  e-*®  erf(&x)dx  =  [e-‘“erf(te)  -  erf  (bx  +  ^)]  (3.15) 

where  a  =  jf-  and  6  =  -.  Eq.  3.14  then  becomes 

VTr  S 

_  DNgO^ 

~  A^ihu^TrS 

-Qrr  erf  (3.16) 
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Evaluating  Eq.  3.16  gives 

^ I-'?’-'  I 

(3.17) 

or  more  simply 

,  DQNgC^^'iE^  ^  ^ 

2y/^{hufs  ' 


3.4.3  TALIF  detection  and  resulting  signal 

The  temporal  shape  of  the  TALIF  signal  collected  in  the  experiment  may  differ  slightly 
from  that  of  N*(t)  due  to  the  instrument  response  function  of  the  detection  system. 
The  temporally  resolved  TALIF  signal,  $(t),  can  be  determined  by  convolving  the 
time  dependent  fluorescence,  N*(t)/rr,  with  the  detection  system  sensitivity  and  re¬ 
sponse  function,  DQ.{t),  to  give 


«(*)  =  ?/’  JV*(C)ft(i-C)<iC. 

Tp  J—OO 


(3.19) 


The  time  response  of  detection  is  assumed  to  behave  as 


fi(f)  = 

0 


(3.20) 


where  S  is  the  instrument  response  time  constant.  Inserting  Eq.  3.11  and  3.20  into 
Eq.  3.19  gives 


\  J  A  /ZT/ L-AO  ^ 


Ay/i^{hvY5rrS 


f  f' 


erffJ-4- 

Vs  2Qr, 


-'j+l 

(3.21) 


f  [ 

J—OO  1^ 


erffS-^ 

Vs  2QTr 


-t-1  (3.22) 


Substituting  x  =  C  -  s^/2QTr  and  f  =  t  -  s^/2QTr  gives 


*(^)  = 

Ay/T^{hvY5TrS 


f_^  [erf  -Hi]  (3.23) 
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where  the  general  formula  of  Eq.  3.15  is  again  used  to  give 


^y/Ti{hvY5TrS 


1 _ 1_ 

5  Qtt 


erf  (I)  +  l]  -  e4«-*>’erf  g  -  ±  +  ^)  |_^  (3.24) 


and  evaluating  while  assuming  5  <  Qr^, 


m = 


A-^{hvf5TrS  J  ^ 


erf  ^  +  1  —  q\(s  qtt^^  erf 


If  s  s 

7”  2^ 


.  (3.25) 


Resubstituting  for  t'  gives 

~  Ay/it{hl/fTrS{l  -  ^ 


(3.26) 


which  is  the  temporally  resolved  TALIF  signal  accounting  for  the  instrument  response 
function. 

An  example  of  using  Eq.  3.26  to  fit  the  experimental  TALIF  signal  is  shown  in 
Fig.  8.  The  TALIF  upper  state  decay  time  was  determined  by  the  fit,  while  the  laser 
FWHM  was  measured  independently  with  a  photodiode  to  be  6.0  ns  and  the  temporal 
response  of  the  detection  system  was  measured,  as  discussed  in  Section  3.4.4,  to  be 
6.64  ns. 


3.4.4  Determination  of  instrument  time  constant 

The  temporal  response  of  the  detection  system  was  limited  by  the  electron  transit 
time  spread  within  the  photomultiplier  tube.  The  instrument  time  constant,  S,  was 
determined  experimentally  by  monitoring  pulse  laser  light  with  the  detection  system 
since  the  decay  time  of  the  laser  pulse  should  be  much  shorter  than  6.  Convolving 
the  laser  temporal  profile  with  the  instrument  response  function  gives 

A(t)  =  f‘  FiCMt  -  C)rfC  (3-27) 

J  — CX3 
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Figure  8:  Experimental  TALIF(207  nm)  data  with  theoretical  fit  where  the  upper 
state  decay  time  at  3  Torr  is  determined  to  be  22.8  ns 
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where  A(t)  is  the  laser  scatter  signal  as  measured  by  the  detection  system.  Inserting 
F{t)  and  Ct{t)  as  defined  earlier, 


A(t)  =  ^ 

Sy/^sJ-oo 

A(t)  =  f‘  (3.28) 

S\/27CS  J-oo 


which  is  solved  with  the  general  formula  of  Eq.  3.9  with  a  =  l/2s^  and  b  =  —1/26, 


giving 

A(t)  =  fe<'/-‘[erf(^-^)H.l]e-/^  (3.29) 

This  function  may  be  used  to  fit  the  data  from  a  laser  scatter  detection  and  solve  for 
the  instrument  time  constant. 


3.4.5  Calibration  of  TALIF  to  absolute  N  atom  concentration 

The  absolute  atomic  nitrogen  concentration  was  calibrated  by  NO  titration  [22],  with 
the  titrant  gas  introduced  in  the  fiowing  afterglow  of  the  microwave  discharge,  just 
prior  to  the  TALIF  viewing  region.  An  example  plot  of  the  TALIF  signal  intensity 
as  a  function  of  the  amount  of  NO  introduced  is  shown  in  Fig.  9  for  a  1  Torr  and  2 
Torr  N2  microwave  discharge.  The  fast  titration  reaction  proceeds  as 

NO  +  N  N2  +  O.  (3.30) 


As  the  controlled  flow  of  NO  is  slowly  increased,  the  vanishing  of  the  TALIF  signal 
denotes  the  exhaustion  of  atomic  nitrogen.  The  known  amount  of  NO  introduced  to 
the  system  when  the  signal  reaches  badcground  indicates  the  absolute  concentration 
of  atomic  nitrogen  that  corresponds  to  the  full  TALIF  signal  prior  to  titration.  In 
Fig.  9,  the  TALIF  signal  was  calibrated  so  that  0.073  arbitrary  imits  of  TALIF  signal 
represented  IxlO^^cm”®  of  atomic  nitrogen  for  that  particular  data  set. 
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Figure  9:  Examples  of  titration  of  atomic  nitrogen  by  the  addition  of  nitrous  oxide 
in  the  flowing  afterglow  of  a  microwave  discharge 


3.4.6  Fluctuation  of  N  atom  TALIF  signal 

The  fluctuation  of  TALIF  signal  strength  between  laser  pulses  is  considered  in  this 
section.  Depending  on  the  conditions  there  may  be  many  sources  that  contribute  to 
the  fluctuation  of  the  fluorescence  signal  measured.  The  two  major  sources  in  this 
experiment  were  determined  to  be  laser  pulse  energy  fluctuation  and  quantum  effects 
from  the  fluorescence  photon  collection.  Some  other  sources  considered,  but  found 
to  be  of  lesser  importance,  were  contributions  from  fluctuations  in  atomic  nitrogen 
concentration,  discharge  emission  background,  laser  scatter,  and  PMT  dark  noise. 

The  contribution  of  quantum  effects  to  the  TALIF  signal  fluctuation  can  be  an¬ 
alyzed  by  a.«=isnTning  that  the  number  of  fluorescence  photons,  /,  collected  for  each 
pulse  follows  a  Poisson  probability  distribution  with  a  mean  value  /.  The  variance, 
a^,  of  this  distribution  is  given  as  [28,  29] 

<^  =  ((/-7)')  =  7-  (3.31) 

Since  the  TALIF  signal  intensity,  I,  is  directly  proportional  to  /,  it  follows  that  the 

variance  in  I  due  to  quantum  effects,  related  as 

^{Quantum)  ^  ^  (3.32) 

The  contribution  of  laser  fluctuations  to  the  TALIF  signal  fluctuation  involves 
more  complexity  than  the  quantum  effects.  This  is  because  the  upper  TALIF  state  is 
populated  by  a  two-photon  process  and,  therefore,  the  resulting  fluorescence  photon 
count  will  not  follow  the  same  distribution  as  the  statistical  distribution  of  the  laser 
photons.  The  number  of  incident  laser  photons  within  one  laser  pulse  can  be  defined 
for  statistical  purposes  as  x,  with 

^  =  Aa  (3.33) 

(H 

where  E  is  the  laser  energy  density  per  pulse  in  J/cm^,  A  is  the  cross  sectional  area  of 
the  laser  focus  region,  and  hi^  is  the  laser  photon  energy.  If  we  assume  the  laser  photon 
energy  to  be  coherent  radiation,  the  photon  statistics  will  follow  that  of  a  coherent 
state,  which  is  a  Poisson  distribution[30].  The  Poisson  probability  distribution,  P, 
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for  the  laser  radiation  can  be  written 

00  .,x 

=  (3.34) 

x=0  *• 

where  ii  is  the  expectation  value,  <  a:  >,  or  the  mean  number  of  incident  laser  photons 
per  pulse.  By  combining  Eq.  3.18  and  Eq.  3.33  the  TALIF  signal  intensity,  J,  can  be 
expressed  as 

^  ^  {hu)^  ^  ^  (3.35) 

The  fluctuation  of  TALIF  signal  due  to  the  laser  fluctuation  would,  therefore,  produce 
a  variance,  o’/(ioser)>  proportional  to.  the  variance  in  the  quantity  which  is 

written 

-  {^)r)  =  (**)-  {^))  +  {x^y  =  .  (3.36) 

In  solving  this,  first  the  expectation  values  of  x^  and  x^  are  found  by 

x=0 

/  \  ^  up 

v)  =  |e~'*  =  /i^  +  6At^  +  7//2  +  ;i  (3.37) 

x=0 

SO  that 

^{Laser)  ^  +  //.  (3.38) 

Since  the  mean  number  of  laser  photons  per  pulse,  /i,  is  a  very  large  integer,  it  «^fln 
be  assumed  that  p?  ^  ^  fx.  The  following  approximations  for  the  mean  TALIF 

signal  intensity  and  its  variance  are  therefore  valid,  extending  from  Eq.  3.35,  Eq. 
3.37  and  Eq.  3.38  , 

I  (X 


^{Laser)  ^ 

(3.39) 

Combining  the  the  above  expressions  gives 

^^(Laser)  OC  I  ^ 

(3.40) 
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The  variance  in  the  experimentally  measured  TALIF  signal  intensity  will  be  a  sum  of 
the  variances  of  all  the  contributing  effects.  For  this  experiment,  the  variance  in  the 
measured  I  is 


2  2  2 

“  ^I{La$er)  “I"  ^/(Quantum) 


(3.41) 


where  o  and  6  are  constants. 

A  set  of  TALIF  signal  strengths  was  experimentally  generated  for  analysis  of  the 
signal  fluctuation.  The  integrated  TALIF  signal  intensity,  J,  is  expressed  in  units 
of  picoVolt-seconds  (pVs)  as  it  appears  on  the  digitizing  oscilloscope.  The  digitized 
signal  is  averaged  for  128  laser  pulses,  from  whidi  the  average  signal,  7,  and  the 
standard  deviation,  a/,  are  determined.  The  variance  of  I  is  plotted  as  a  function 
of  7  over  a  range  of  signal  strengths  in  Fig.  10.  The  TALIF  signals  were  varied 
over  a  range  of  several  orders  of  magnitude  by  adjusting  the  discharge  conditions  in 
a  microwave  cell  to  vary  the  downstream  atomic  nitrogen  density.  Also  shown  in 
Fig.  10  is  a  flt  of  Eq.  3.41  to  the  experimental  data  (solid  line)  and  each  major 
contribution  term  shown  independently  (dashed  lines).  The  results  show  the  lower 
level  signals  are  dominated  by  quantum  effect  fluctuation  and  that  the  effects  of  laser 
fluctuation  eventually  take  over  at  higher  signal  levels. 

These  different  regions  of  fluctuation  behavior  become  influential  when  considering 
the  effect  of  correcting  each  individual  TALIF  signal  pulse  for  its  corresponding  laser 
pulse  energy.  This  laser  correction  can  be  friiitful  when  applied  to  the  higher  level 
signals,  virtually  eliminating  the  variance  due  to  laser  fluctuation  and  leaving  a  small 
remaining  standard  deviation.  The  average  signal,  I,  can  then  be  obtained  to  a  small 
uncertainty  with  minimal  averaging  in  this  case.  On  the  other  hand,  the  lower  level 
signals  would  not  beneflt  from  eliminating  the  variance  due  to  laser  fluctuation  since 
^(Quantum)  ^  major  Contributor  in  this  range.  It  is  clear  that  averaging  over  a 
substantial  number  of  signal  pulses  will  be  necessary  to  obtain  an  accurate  measure 
of  I  for  small  signal  levels. 
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CHAPTER  4 


Alternative  excitation  scheme  for  atomic  nitrogen 

TALIF 


4.1  Introduction 

In  an  attempt  to  improve  the  performance  of  N  atom  TALIF  techniques,  this  chapter 
details  a  search  for  an  alternative  to  the  traditional  N  atom  TALIF  scheme  described 
in  Section  2.2.  The  well  known  N  atom  TALIF  scheme  that  employs  laser  absorption 
at  211  nm  has  several  limitations.  The  fluorescence  wavelength  of  868  nm  is  within 
a  spectral  range  of  poor  sensitivity  for  many  common  photomultiplier  detectors.  In 
addition,  868  nm  coincides  with  intense  molecular  N2  emission  that  is  often  associated 
with  a  discharge  source.  The  presence  of  this  emission  can  significantly  lower  the 
sensitivity  of  the  TALIF  diagnostic.  Finally,  the  upper  state  of  this  sdieme 
is  known  to  have  a  relatively  large  collisional  quenching  rate.  The  uncertainties 
associated  with  large  corrections  in  signal  strength  made  for  quenching  efiects  limit 
the  sensitivity  and  accuracy  of  N  atom  concentration  measurements. 

The  traditional  TALIF  technique  excites  the  atomic  nitrogen  from  the  ground 
(2p^)  state  to  the  (3p)  with  a  laser  wavelength  of  211  nm.  The  seardi 
for  an  alternate  TALIF  scheme  began  by  considering  upper  states  with  higher  energy 
than  (3p)  ^D°.  As  seen  in  the  N  atom  energy  level  diagram  in  Fig.  11,  (3p)  and 
(3p)  ^S°  states  exist  at  95532  cm“^  and  96750  cm  respectively.  The  (3p)  ^P°  state 
will  not  be  considered  since  the  (3p)  ^P°  4-  (2p®)  ^83/3  laser  excitation  is  not  allowed 
by  two-photon  transition  selection  rules,  as  will  be  proven  in  Section  4.4.4.  The  two- 
photon  absorption  from  the  ground  to  the  (3p)  state  is  allowed  and  requires  a  laser 
wavelength  of  207  nm,  which  is  no  more  difficult  to  generate  than  the  211  nm  laser. 
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The  fluorescence  from  this  new  TALIF  scheme  occurs  at  ~745  nm  for  the  (3p)  ^S|^2 
(3s)  transition.  Laser  excitation  of  the  207  nm  transition  has  only  been  mentioned 
once  in  previous  literature,  when  Bengtsson  et  a/.  [13]  utilized  the  scheme  to  measure 
the  radiative  lifetime  of  the  (3p)  ^S|^2  state.  No  further  experimentation  was  reported 
with  the  207  nm  technique  and  no  discussion  was  included  on  the  viability  of  this 
TALIF  scheme  compared  to  the  traditional  211  nm  scheme. 

La  this  chapter,  the  performance  of  the  alternative  TALIF  scheme  with  207  nm 
absorption  is  compared  to  the  well  known  TALIF  scheme  with  211  nm  absorption. 
The  new  TALIF  scheme  with  a  shorter  fluorescence  wavelength  clearly  eases  the  re¬ 
quirements  of  infrared  spectral  range  on  the  detector.  The  new  scheme’s  fluorescence 
at  745  nm  lies  in  a  spectral  region  of  plasma  emission  that  is  found  to  have  compa¬ 
rable  intensity  with  that  of  the  868  nm  region  where  the  traditional  TALIF  scheme 
fluoresces.  This  leaves  both  schemes  equally  unfortunate  in  being  nearly  coincident 
with  strong  molecular  nitrogen  emission. 

The  relative  efficiencies  of  the  two  TALIF  techniques  ultimately  depend  on  the 
TALIF  signal  intensities.  By  direct  experimental  comparison  of  the  TALIF  signals, 
the  new  scheme  with  excitation  at  207  nm  was  found  to  have  a  significantly  larger 
excitation  cross-section.  In  addition,  at  higher  N2  pressure,  the  advantage  of  the  new 
TALIF  scheme  increases  due  to  the  low  N2  quenching  rate  of  the  (3p)  ^S|/2  state, 
which  was  found  in  this  work  to  be  nearly  an  order  of  magnitude  smaller  than  that 
of  the  (3p)  ^D?/2  level. 

4.2  Theory 


Whenever  a  TALIF  signal  measurement  is  discussed  to  in  this  work,  it  refers  to  time 
integrated  fluorescence  intensity,  I,  which  was  collected  by  the  PMT  and  the  output 
digitized.  This  integrated  TALIF  intensity  and  its  relation  to  other  system  parameters 
was  analyzed  in  detail  in  Section  3.1  where  it  was  concluded  in  Eq.  3.18  that 

DQN^a^ 

2y/^{hvys  ’ 

where  D  is  the  signal  generated  per  fluorescence  photon,  Q  is  the  fluorescence  quan¬ 
tum  5deld,  Ng  is  the  ground-state  atom  population,  is  the  two-photon  absorption 
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associated  wavelengths. 
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rate  coefficient  in  cm^s,  E  is  the  laser  energy  density  per  pulse  in  J/cm^,  and  hv  is 
laser  photon  energy.  The  laser  temporal  width  parameter,  s,  is  defined  so  that  the 
Gaussian  FWHM  of  the  laser  temporal  profile  is  2s^2/n(2). 

When  experimentally  comparing  the  TALIF  schemes  with  207  nm  and  211  nm 
laser  excitation,  the  following  conditions  were  met.  Care  was  taken  to  maintain  the 
ground  state  atom  population,  Ng,  at  a  constant  level.  Temporal  profile  measurements 
of  the  207  nm  and  211  nm  laser  pulses  indicated  that  s  was  the  same  for  both  schemes. 
The  incident  laser  power,  as  measured  by  a  joulemeter,  was  adjusted  so  that  E  was 
also  the  same  for  both  schemes.  This  allows  the  fluorescence  time  integrated  signal 
intensity  ratio  to  be  simplified  to 

■^207  _  <^207  -^207  ^207  . 

■^211  (T211  D211  Q211  I'm 

where  the  factor  *^211/^207  nc&r  unity  with  a  value  of  0.06.  As  seen,  the  TALIF  signal 
strengths  may  differ  between  the  two  schemes  as  a  result  of  contrast  in  the  two-photon 
absorption  rates,  the  detection  efficiencies,  and/or  the  fluorescence  quantum  yields. 

4.3  Experimental 

The  generation  of  tunable  radiation  to  provide  the  207/211  nm  laser  pulse  energy 
for  the  two  TALIF  schemes  compared  in  this  section  is  described  in  detail  in  Section 
3.1.  Also  detailed  in  Chapter  3  are  the  microwave  discharge  cell  and  gas  hanHling 
configuration  and  the  fluorescence  detection  system.  When  experimentally  comparing 
the  TALIF  schemes,  a  constant  atomic  nitrogen  concentration  was  maintained  in 
the  flowing  afterglow  region  of  the  microwave  disdiarge  over  a  varied  range  of  N2 
pressure.  The  N2  flow  rate  and  microwave  input  power  were  adjusted  to  levels  that 
were  calibrated  by  titration,  discussed  in  Section  3.4.5,  to  provide  a  constant  N  atom 
concentration. 
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4.4  Results  and  Discussion 

4.4.1  Discharge  Emission  Background 

The  continuous  plasma  emission  badcground  demands  attention  whenever  an  N  atom 
TALIF  technique  is  conducted  within  an  N2  discharge  region.  A  typical  plasma  emis¬ 
sion  spectrum  for  a  pure  N2  100  Watt  microwave  discharge  is  shown  in  Fig.  12,  along 
with  the  spectral  positions  of  the  fluorescence  at  868  nm  and  745  nm  respectively  from 
the  TALIF(211)  excitation  to  the  (3p)  upper  state  and  the  TALIF(207)  excitation 
to  the  (3p)  *S°  upper  state.  As  Fig.  12  shows  the  TALIF(211)  fluorescence  directly 
coincides  with  the  peak  emission  from  the  B  ®n5(v^=2)  — >  A  (V^=l)  bandhead 
of  N2.  The  N2  first  positive  emission  is  coincident  as  well  with  the  TALIF(207)  sig¬ 
nal,  mth  the  fluorescence  multiplet  seen  to  fall  on  the  short  wavelength  shoulder  of 
the  B  *nj(v'=4)  ->  A  ^E+(v"=2)  bandhead.  This  demonstrates  that,  unfortunately, 
continuous  backgroimd  noise  from  strong  molecular  bandheads  is  encountered  when 
collecting  the  fluorescence  of  either  TALIF  technique. 

4.4.2  Experimental  Comparison  of  TALIF  Signals 

While  the  background  interference  is  comparable  for  the  two  schemes,  the  advan¬ 
tage  of  using  TALIF(207)  over  TALIF(211)  is  firmly  established  in  the  experimental 
comparison  of  the  TALIF  signal  strengths.  The  TALIF  signals  were  analyzed  in  the 
flowing  afterglow  region  of  the  microwave  discharge  with  the  N  atom  concentration 
maintained  at  approximately  6x10^^  cm'^.  Each  TALIF  signal  was  obtained  with 
incident  laser  pulse  energy  of  150  //J.  Typical  temporally  resolved  fluorescence  sig¬ 
nals  for  the  two  N  atom  schemes  at  a  pressure  of  3  Torr  of  N2  are  shown  in  Fig. 
13.  The  integrated  intensity  ratio,  I207/I211J  for  this  pair  of  signals  is  ~10.  In  Table 
1,  l2or/l2ii  integrated  intensity  ratios  are  listed  for  various  background  N2  pressures 
from  2  Torr  to  6  Torr.  The  TALIF(207)  signal  is  seen  to  be  considerably  greater  than 
the  TALIF(211)  signal  for  each  condition.  The  detector  sensitivity  ratio,  D207/D2U, 
is  also  given  in  Table  1,  and  includes  effects  of  interference  filter  transmission  and 
photomultiplier  tube  spectral  sensitivity.  The  GaAs  cathode  of  the  RCA-31034  pho¬ 
tomultiplier  tube  provides  the  best  available  infrared  spectral  response  curve  out  to 
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Figure  12:  Pure  N2  microwave  discharge  emission  spectrum  with  atomic  nitrogen 
fluorescence  lines  from  both  TALIF  schemes  shown  as  stick  figures. 
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Figure  13:  Temporally  resolved  N  atom  TALIF  signals  with  excitation  by  207  nm 
and  211  nm  with  equivalent  laser  power  and  beam  shape  in  3  Torr  of  N2. 
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the  TALIF(211)  fluorescence  at  868  nm.  If  one  elected  to  use  a  more  common  pho¬ 
tomultiplier  tube  that  is  less  sensitive  at  868  nm,  the  sensitivity  ratio  and  therefore 
the  signal  intensity  ratio  would  be  even  higher  in  favor  of  the  TALIF(207)  scheme. 


Pressure 

2  Torr 

3  Torr 

4  Torr 

5  Torr 

6  Torr 

Signal  Intensity 
I207/I211 

8.6±0.9 

10.4±1.0 

11.6±1.2 

12.8±1.3 

14.8±1.5 

Detection  Sensitivity 
D207/D2H 

1.20±0.14 

1.20±0.14 

1.20±0.14 

1.20±0.14 

1.20±0.14 

Quantum  Yield 
Q207/Q211 

2.1±0.2 

i6±0.3 

3.0±0.4 

3.4±0.5 

3.8±0.6 

Excitation  Rate 
(Exp) 

3.5±0.6 

3.5±0.6 

3.4±0.7 

3.3±0.7 

3.4±0.6 

Table  1:  TALIF(207):TALIF(211)  Comparison. 


4.4.3  Collisional  Quenching  of  TALIF  Upper-State 

In  any  case,  D207/D211  accounts  for  only  a  relatively  small  fraction  of  the  departures 
of  I207/I211  from  unity.  Furthermore,  the  I207/I211  signal  ratio  in  Table  1  is  seen  to 
increase  as  the  background  N2  pressure  increases.  This  phenomenon  is  due  to  difier- 
ent  collisional  quenching  rates  of  the  upper  (3p)  ^S|/2  and  (3p)  N  atom  states 
by  N2.  A  Stem-Volmer  plot,  shown  in  Fig.  14,  quantifies  the  collisional  quenching 
behavior  of  both  TALIF  signals  with  the  relationship 

—  =  kgW^]  +  i  (4.3) 

where  the  time  constant  for  fluorescence  decay,  r^,  is  determined  by  least  squares 
fit  of  the  signal  to  the  time  dependent  fluorescence  convolved  with  the  instrument 
response  function  of  the  detection  system.  The  collisionless  radiative  decay  time, 
Tp,  can  be  determined  from  the  intercept  of  the  Stem-Volmer  plot  at  zero  pressure. 
For  the  (3p)  ^D°  upper  state,  it  is  determined  that  Tr  =  41(5)  ns,  and  for  the  (3p) 
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^S°  upper  state  it  is  found  that  Tt  =  26(2)  ns.  Both  of  these  measurements  agree 
well  with  the  most  recent  studies[13,  31].  The  slope  of  each  plot  represents  the  rate 
of  the  upper  state  quenching  by  background  N2  gas.  The  quenching  rate  constant, 
kg,  of  the  respective  N  atom  upper  states  by  N2  are  determined  to  be  kg(^D7/2)  = 
1.63(16) xlO’’  (Torr-s)-^  =  5.1(5) xlO"^"  cmVs  and  =  2.1(3) xlO®  (Torr-s)-^ 

=  6.7(9)  X 10"^^  cm®/s.  Our  measurement  is  within  experimental  uncer¬ 
tainty  of  the  value  =  4.6(6)  x  10“^°  cm®/s  determined  by  Copeland  et  al. 

[31]  The  value  of  kg(^5|/2)  has  not  been  previously  reported. 

The  effect  of  collisional  quenching  on  the  TALIF  signal  intensity  is  best  quan¬ 
tified  by  the  quantum  yield  of  the  transition,  as  given  in  Eq.  3.2.  The  ratio  of 
the  quantum  yield  values  between  the  two  transitions,  Q207/Q2ii>  is  shown  in  Ta¬ 
ble  1  for  each  pressure.  The  quantum  yield  ratio  is  seen  to  increase  with  increasing 
pressure  due  to  relatively  efficient  collisional  quenching  of  the  (3p)  state.  The 
quenching  rate  of  the  (3p)  ^83^2  state,  on  the  other  hand,  gives  rise  to  only  minimal 
change  in  the  quantum  yield  over  a  large  pressure  range.  This  indicates  that  a  strong 
TALIF(207)  signal  can  be  maintained  over  a  broad  range  of  background  pressure 
without  significant  correction  terms  for  quenching,  whereas  the  TALIF(211)  signal 
strength  diminishes  significantly  with  increasing  pressure. 

The  ratio  of  the  two-photon  excitation  rate  coefficients,  <^207/0211,  for  the  two 
TALIF  schemes  can  be  foimd  from  Eq.  4.2  to  be 


7207  <?2ii  L)2ii  vloi 
7211  Q207  7)207  *^211 

Table  1  includes  the  ratio  of  two-photon  excitation  rate  coefficients  calculated  from 
Eq.  4.4  for  each  pressure  using  the  I207/I211J  D207/D2ii>  and  Q207/Q211  data  contained 
therein.  In  contrast  to  the  integrated  signal  intensities,  the  relative  cross  sections 
are  seen  to  be  independent  of  pressure,  as  is  expected,  with  a  o’^7/o^ii  value  of 
approximately  3.5. 


4o7/4ii 
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Pressure  (Torr) 


Figure  14:  Stem-Volmer  plot  of  atomic  nitrogen  (3p)  and  (3p) 

decay  with  quendiing  by  N2. 


4.4.4  Comparison  of  Theoretical  Two-Photon  Absorption  Cross 
Sections 


A  theoretical  calculation  of  the  two-photon  absorption  cross  sections  of  both  of  the 
atomic  nitrogen  TALIF  sdiemes  is  considered  here  as  a  further  comparison  of  the 
techniques.  The  validity  of  such  a  comparison,  though,  is  subject  to  the  soundness  of 
the  approximations  made  in  order  to  keep  the  calculation  feasible.  The  two-photon 
absorption  cross  section  to  be  calculated,  in  units  of  cm'*,  depends  only  on  the 
properties  of  the  absorbing  atom  and  not  the  properties  of  the  incident  laser[32].  This 
cross  section  is  related  to  the  two-photon  excitation  coefficient,  discussed  early 


as 


(7(2)  = 


(4.5) 


where  g(u)  is  the  normalized  laser  line-shape  function  and  is  the  photon  statistical 

factor,  conamonly  assumed  to  be  equal  to  2  for  simultaneous  two-photon  absorption 
from  a  multimode  laser[26,  32].  The  values  of  g{u)  and  G^^^  can  be  assumed  to 
remain  constant  when  alternating  between  207  nm  and  211  nm  laser  radiation  in 
this  experiment.  This  allows  that  the  ratio  of  the  calculated  theoretical  two-photon 
absorption  cross  sections,  to  directly  compared  to  the  experimentally 

determined  ratio  of  two-photon  excitation  coefficients,  cr^r/cr^i,  in  Table  1.  The 
two  photon  absorption  cross  section  from  the  ground  state  may  be  obtained  from 
perturbation  theory  as  [33,  34] 


M,M' 


{J^M^\z\i){i\z\JM) 


1/4  —  1/ 


(4.6) 


where  1/  is  the  laser  frequency  and  e  is  the  fundamental  charge.  The  \JM)  eigenstates 
are  labeled  so  that  the  prime  designates  the  excited  state  and  unprimed  symbols  relate 
to  the  ground  state  of  the  atom.  A  virtual  intermediate  state  is  represented  by  a 
summation  over  a  complete  set  of  intermediate  states,  ji),  and  hz/j  is  the  energy  of  the 
ith  intermediate  state  above  ground.  For  this  type  of  calculation  to  be  feasible,  the 
usual  convention  is  to  choose  a  single  intermediate  eigenstate  of  the  atom  that  provides 
a  large  contribution  to  the  sum  of  the  dipole  moment  products  over  all  intermediate 
states.  Previous  theoretical  calculations  of  approximate  two-photon  absorption  cross 
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sections  of  light  atoms  have  been  mostly  directed  toward  atomic  oxygen  [34].  These  0- 
atom  TALIF  calculations  have  shown  that  the  single  intermediate  state  approximation 
produces  results  in  good  agreement  with  atomic  oxygen  TALIF  experiments. 

The  two-photon  cross  section  in  Eq.  4.6  may  be  expanded  with 


\JM)=  {LMlSMs\JM)\LMlSMs)  (4.7) 


and  the  sum  over  intermediate  states  approximated  by  a  single  eigenstate  to  give 


(ftc)2  2J-fl 


S  J2{LMlS{M-Ml)\J'M) 


‘  Ml 


{LMlS{M  -  Ml)\JM) {UMl \z\JM) 


Vi  — V 


(4.8) 


noting  that  AMi=0  for  linear  polarized  laser  absorption.  Assigning  ground  state, 
(2p^)  and  intermediate  state,  (3s)  ^P,  quantum  values,  the  cross  section  simpli¬ 
fies  to 


(2)  _ 


I(£'ON|10)P  |(10|^|00>p.  (4.9) 


The  dipole  moment  matrices  can  be  separated  into  radial  and  angular  components 
by 

\{L'Ml\z\LM^)\^  =  \{L'\z\L)\^  (LMaa\L'Mz)\  (4.10) 


The  radial  dipole  matrix  term  can  be  evaluated  as[35,  36] 


~  3^  (2£  +  1) 


(4.11) 


where  A  is  the  wavelength  of  the  transition  and  Al',l  is  the  Einstein  A  coefficient 
between  the  two  multiplet  terms,  V  and  L.  The  two-photon  cross  section  then  becomes 


<r?)  = 


10247r^c2  {ui  - 1/)2  ^  {L'0^M\J'M)^  |(1010|L'0)p 


(4.12) 

For  the  TALIF  schemes  with  ^D°  and  upper  states,  having  L'  =  2  and  L'  = 
0  respectively,  all  of  the  terms  in  Eq.  4.12  are  non-zero.  This  indicates  that  the 
(3p)^D7y2  {2p®)^S3y2  (3p)  ^83/2  (2p^)^S3^2  absorptions  are  both  allowed 
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by  two-photon  transition  selection  rules.  When  considering  a  ^P°  upper  state,  though, 
the  quantum  value  L'  =  1  is  applied  which  makes  the  term  |(1010|L^0)p  in  Eq.  4.12 
equal  to  zero  and  therefore  the  two-photon  cross  section  is  zero.  This  confirms  the 
statement  in  Section  4.1  that  the  transition  (3p)  -4—  (2p®)  ^83^2  allowed 

by  two-photon  transition  selection  rules.  By  utilizing  the  transition  probability  data 
from  Wiese  et  al.,  Eq.  4.12  produces  the  numerical  results  shown  in  Table  2. 


Two-Photon  Transition 

Excitation  Wavelength 

Two-Photon  Cross  Section 

(3p)  (2p»)  ‘S|„ 

211  nm 

2.86  X  10-3®cm^ 

(3p)  <-  (2p»)  ‘SI,; 

- 

0 

(3p)  ^83/2  (2p*)  ^83/2 

207  nm 

3.91  X  10-^“cm^ 

Table  2:  Theoretical  Cross  Section  Comparison. 


The  theoretical  derivation  produces  a  ratio  of  This  differs  be¬ 

yond  the  uncertainty  tolerance  with  the  experimentally  determined  ratio  of  ct^qj/ o'lii  = 
3.5.  This  difference  in  two-photon  absorption  rate  comparison  may  be  attributed  to 
the  approximation  made  in  the  theory  that  the  single  (3p)^P  eigenstate  of  the  N 
atom  can  replace  the  sum  over  all  intermediate  states  in  the  general  expression.  For 
atomic  oxygen,  this  approximation  has  been  shown  through  detailed  calculation  and 
experiment  to  work  well,  but  no  data  exists  for  atomic  nitrogen.  It  is  quite  possible 
that  a  vast  number  of  ^P  states  with  energy  levels  significantly  greater  than  that  of 
(3p)^P  combine  to  contribute  significantly  to  the  sum  of  intermediates,  which  would 
affect  the  accuracy  of  this  theoretical  calculation. 

4.5  Conclusion 

An  alternative  two-photon  laser  induce  fiuorescence  scheme  for  ground  state  atomic 
nitrogen,  labeled  TALIF(207),  has  been  compared  with  the  well  known  TALIF(211) 
scheme.  The  TALIF(207)  technique  has  been  found  experimentally  to  be  superior  to 
TALIF(211)  in  two-photon  absorption  rate  by  a  factor  of  3.5.  A  theoretical  derivation 
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of  the  same  two-photon  absorption  rate  ratio  was  also  done,  with  the  caveat  that  cer¬ 
tain  approximations  were  applied  with  respect  to  the  intermediate  state  of  the  TALIF 
schemes.  This  theoretical  calculation  also  demonstrated  the  TALIF(207)  technique 
to  be  favorable,  but  suggests  a  factor  of  1.4.  It  was  shown  that  the  TALIF(207) 
scheme  possesses  an  upper  state  quenching  rate  by  N2  almost  an  order  of  magni¬ 
tude  lower  than  TALIF(211).  The  TALIF(207)  scheme  was  also  found  to  be  slightly 
more  attractive  for  detection  with  near  visible  fluorescence,  while  being  equivalent 
to  TALIF(211)  in  laser  generation  feasibility  and  susceptibility  to  discharge  emission 
background  noise.  In  general,  the  TALIF(207)  scheme  should  provide  a  much  im¬ 
proved  signal  strength  under  most  N2  discharge  conditions  with  decreased  concern 
for  corrections  for  collisional  quenching  effects,  a  process  that  is  always  challenging. 
With  these  beneflts  and  no  observed  disadvantages,  it  is  strongly  recommended  that 
the  TALIF(207)  scheme  be  utilized  in  place  of  TALIF  (211)  for  remote  detection  and 
analysis  of  ground  state  atomic  nitrogen. 
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CHAPTER  5 


Photodissociation  of  N2O  and  subsequent  N  atom 

TALIF  detection 


5.1  Introduction 

A  technique  of  molecular  photodissociation  that  produces  atomic  nitrogen  with  sub¬ 
sequent  detection  of  the  N  atom  with  the  same  laser  by  TALIF  would  be  of  interest  in 
the  field  of  photodiemistry.  This  technique  would  also  have  a  general  application  in 
the  calibration  of  N  atom  TALIF  experiments.  With  a  uniform  distribution  of  precur¬ 
sor  molecules  in  a  cell  without  a  discharge,  the  laser  source  at  the  TALIF  wavelength 
would  generate  a  constant  density  of  N  atoms  at  the  focus  and  then  immediately 
detect  them  with  TALIF.  This  would  provide  a  spatially  uniform  TALIF  signal  that 
could  be  Tised  to  generate  a  view-factor  correction  for  any  cell  positions  where  the 
fiuorescence  was  partially  blocked. 

Nitrous  oxide,  N2O,  has  been  previously  reported  to  be  a  suitable  precursor  for 
producing  ground  state  N  atoms  by  photodissociation  with  subsequent  TALIF  ab¬ 
sorption  to  the  (3p)  ^Ty^i2  or  (3p)  '*83/2  state[13].  An  energy  level  diagram  for  N2O 
is  presented  in  Fig.  15(a)  along  with  an  adiabatic  correlation  diagram  in  Fig.  15(b) 
which  shows  the  dissociation  pathways  firom  N2O  excited  states  below  6  eV[37]. 

Only  the  N2O  excited  states  that  have  C5,  or  bent,  symmetry  are  energetically 
accessible  with  207  nm  (6  eV)  radiation  absorption  from  the  ground  state.  The  N2O 
ground  state  S5unmetry  is  in  the  Coo»  point  group  which  corresponds  to  ^  A'  in  the 
Cs  point  group.  Therefore,  transitions  to  the  2^A'(B)  and  1^A"(A)  excited  states  are 
allowed  from  the  ground  state,  where  as  transitions  to  the  l^A'(a)  and  l^A"(b)  states 
are  not  spin  allowed.  The  four  combinations  of  dissociative  products  shown  on  the 
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Figure  15:  (a)  Energy  level  diagram  for  N2O  states  of  Coo»  symmetry  (linear)  and 
symmetry  (bent)  (b)  Adiabatic  correlation  diagram  showing  dissociation  pathways 
from  N2O  states  below  6  eV,  which  all  involve  bent  N2O  excited  states 
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adiabatic  correlation  diagram  are  the  possible  photofragments  from  absorption  of  207 
Tim  radiation.  Since  the  photodissociation  channel  resulting  in  N2  and  0(^D)  involves 
allowed  transitions,  this  combination  of  photofragments  would  be  produced  much 
more  readily  than  NO  and  N(^S)  photofragments.  The  distribution  of  photofragment 
combinations  produced  from  N2O  with  incident  207  nm  laser  radiation  are  expected 
to  be  effectively  the  same  as  published  results  using  205  nm  dissociation  [38].  The 
primary  photodissociation  processes  and  the  percent  of  photolysis  products  produced 
by  each  are[38,  39] 

N2O  ^  iV2(X^E)  +  0{^D)  (~  90%) 

N2O  H  N2{X^E)  +  OCS)  (~  5%) 

N2O  ^  N2{X^^)  +  OfF)  (-  2%).  (5.1) 

A  secondary  dissociative  branch  that  has  been  shown  to  exist,  but  with  an  upper 
bound  of  2%  of  the  primary  process,  is 

N2O  ^  NOiX^n)  +  iV(^5).  (5.2) 

Despite  the  seemingly  low  yield  of  N  atoms,  the  photolysis  of  N2O  is  a  more  efficient 
source  of  atomic  nitrogen  than  the  photolysis  of  NO2  or  NO  with  A  >  200  nm[39]. 
Bengtsson  and  coworkers  reported  N2O  photolysis  and  two-photon  absorption  [13]  as 
the  first  step  in  a  process  that  further  excited  the  N  atoms  from  the  (3p)  state  to 
a  series  of  Rydberg  states  by  single  photon  absorption  using  a  secondary  dye  laser 
beam.  The  paper  established  success  in  producing  ground  state  N  atoms  in  a  few 
himdred  mTorr  of  N2O  and  detecting  a  TALIF  signal.  The  data  in  the  paper,  though, 
indicated  that  the  upper  (3p)  state  of  the  TALIF  scheme  was  severely  quenched  as  the 
nitrous  oxide  pressure  was  raised  toward  0.5  Torr.  The  quenching  rate  was  not  stated 
but,  by  inspection,  is  obviously  much  greater  than  that  of  N2,  which  was  reported  in 
Chapter  4.  The  authors  claimed  that  a  lack  of  a  well  defined  collision  partner  in  the 
interaction  region  prevented  any  analysis  of  the  quenching.  This  severe  quenching 
could  limit  the  usefulness  of  this  technique,  since  a  higher  pressure  nitrous  oxide 
environment  would  produce  more  N  atoms,  but  the  resulting  TALIF  signal  could  be 
quenched  to  a  degree  that  renders  it  unusable. 
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In  this  work,  nitrous  oxide  photodissociation  and  atoniic  nitrogen  TALIF  detection 
were  both  accomplished  with  207  nm  laser  radiation.  The  quenching  rate  of  the  (3p) 
^83/2  upper  state  in  the  region  of  N2O  photolysis  was  experimentally  determined.  The 
quenching  rate  of  the  (3p)  ^S^/2  state  by  N2O  was  then  independently  investigated 
using  a  microwave  discharge  N  atom  source  and  mixing  N2O  in  the  flowing  afterglow 
region.  The  quenching  of  the  (3p)  ^83^3  state  in  the  photolysis  environment  was 
confirmed  by  experiment  to  be  very  high.  The  (3p)  ^S^/2  quenching  by  N2O  alone  was 
found  to  be  quite  high,  but  less  than  the  effective  quenching  rate  during  photolysis. 
It  is  proposed  here  that  the  N  atom  collision  frequency,  and  thus  the  quenching  rate, 
is  much  higher  in  the  photolysis  experiment  than  in  the  flowing  afterglow  experiment 
due  to  the  excess  kinetic  energy  given  to  the  N  atom  photofragment.  The  NO  density 
produced  by  N2O  photolysis  is  estimated  to  be  a  very  small  fraction  of  the  total  gas 
density  in  this  experiment.  Upon  considering  the  upper  limit  of  the  NO  quenching 
rate,  N2O  is  determined  to  be  the  dominant  quenching  partner  with  NO  only  a 
minor  contributor.  This  is  contrary  to  the  assumptions  by  Bengtsson  et  c/.[13],  where 
analysis  of  quenching  rates  were  avoided  since  the  collision  partners  were  not  well 
defined. 

5.2  Experiment 

The  AsTeX  microwave  discharge  cell,  detailed  in  Section  3.2,  was  utilized  for  this 
set  of  experiments.  The  TALIF(207)  technique  was  used  as  described  in  Chapter  3. 
The  experimental  conditions  included  either  a  cell  with  a  5  seem  slow  flow  of  pure 
nitrous  oxide  with  no  discharge  or  a  flowing  N2  microwave  discharge  with  a  quendiing 
gas  mixed  in  downstream  from  the  discharge.  When  a  quenching  gas  was  added  to 
the  flowing  afterglow  region  of  the  microwave  discharge,  the  N2  partial  pressure  was 
maintained  at  2  Torr.  With  the  TALIF  signal  quenching  effects  of  N2  well  established 
from  data  in  Section  4.4,  the  quenching  effects  of  the  added  gas  could  be  easily 
extracted.  The  analysis  method  to  determine  the  fluorescence  decay  lifetime  from 
the  raw  TALIF  signal  is  derived  in  Section  3.4. 
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5.3  Results  and  Discussion 

5.3.1  Fraction  of  Photolyzed  N2O  in  TALIF  Region 

An  estimate  of  the  number  density  of  N2O  that  experience  photodissociation,  [N2O]*, 
in  the  TALIF  region  can  be  calculated  for  this  experiment.  The  number  of  incident 
photons  per  pulse,  /jl,  relating  to  the  0.2  mJ  of  laser  energy  is  fi=2xl0^^  photons. 
The  number  density,  [N2O]*,  within  the  laser  focal  region  can  then  be  given  as 

(5-3) 

where  la/Io  is  the  fraction  of  laser  light  absorbed  within  the  focal  region,  A  is  the 
cross  sectional  area  of  the  beam  focus  and  i  is  the  length  of  the  focal  region.  The 
Beer-Lambert  law  will  be  assumed  valid  when  applied  to  the  laser  focal  region.  Laser 
light  absorption  may  then  be  expressed  as[39] 

=  (5.4) 

where  k  is  the  absorption  coefficient  of  N2O  in  (Torr  •  cm)“^  and  p  is  the  N2O  pressure. 
With  k  =  10~®(Torr  •  cm)~^  for  207  nm  light[40],  p  <1  Torr,  and  i  <1  mm,  it  is  seen 
that  in  this  experiment  kpl  «C  1.  This  means  that  Eq.  5.4  can  be  approximated  as 
la/h  ^  kpl  which,  combined  with  Eq.  5.3,  gives 

[iV,0]*  «  (5.5) 

For  a  nitrous  oxide  pressure  of  0.5  Torr  and  a  beam  diameter  of  200  pm,  the  result 
is  [N2O]*=3xl0^^cm“^  which  corresponds  to  <2%  of  the  [N2O]  in  the  region.  This 
small  [N20]*/[N20]  estimate  indicates  that  the  N2,  0,  and  NO  photofragments,  as 
well  as  any  molecules  produced  in  fast  secondary  reactions,  are  not  likely  to  be  major 
contributors  to  the  N  atom  (3p)  ^83/2  state  quenching  in  this  experiment. 

5.3.2  Collisional  Quenching  Results 

Analysis  of  the  fluorescence  decay  from  the  N  atom  (3p)  ^83^2  ^  ^1^®  presence 

of  photolyzed  N2O  was  done  for  a  cell  with  a  slow  flow  of  N2O  at  a  pressure  of  50 
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mToiT  to  750  mTorr.  The  decay  constant  determined  in  the  analysis  is  equal  to  l/r^, 
which  is  defined  as 

r  =  +  (5.6) 

where  is  the  radiative  decay  lifetime  and  the  sum  is  over  all  the  species  present  after 
the  N2O  photolysis,  with  [M]  being  the  number  density  and  kq  the  quenching  rate 
of  each  species.  Since  the  density  of  each  of  the  collision  partners  in  the  photolysis 
interaction  region  is  not  known,  an  effective  quenching  rate,  fc' ,  due  to  the  collective 
quenching  of  all  species  is  introduced  as 

(5.7) 

i 

where  fc'  is  coupled  with  the  number  density  of  nitrous  oxide  in  the  cell  prior  to 
photolysis.  The  results  in  Fig.  16  show  that  the  species  in  the  N2O  photolysis 
region  quench  the  N  atom  upper  state  very  rapidly  with  an  effective  quendiing  rate 
of  A:'=5.8xlO’^(Torr-s)“^  or  1.8xlO~®cmVs 

In  order  to  isolate  the  quenching  effects  of  N2O  as  a  collision  partner  with  the 
N  atom  (3p)  ^83/2  state,  N2O  was  introduced  downstream  of  a  2  Torr  N2  microwave 
disdiarge.  The  decay  constant  of  the  TALIF  upper  state  in  this  case  behaves  as 

—  =  kg(^ff20)[^20]  +  ^j(Ar2)[Ar2]  H - .  (5.8) 

*d  Ty. 

The  quenching  rate,  kg^ff^o),  is  found  by  varying  only  the  partial  pressure  of  nitrous 
oxide  in  the  post-discharge  mix.  The  result,  determined  from  the  plot  in  Fig.  16,  is 
that  A:5(Ar20)=l -3x10^ (Torr  •  s)“^  or  4.2xl0“^®cm®/s. 

Since  N2  and  0(^D)  are  the  primary  products  of  N2O  photodissociation  at  207  nm, 
at  least  a  small  density  of  N2  can  be  expected  in  the  N2O  photolysis  cell.  CoUisional 
quenching  of  the  TALIF  upper  state  by  N2(X^S)  is,  therefore,  also  considered  here. 
The  quenching  of  N  atom  (3p)  ^83^2  ^2  was  studied  in  detail  in  Section  4.4  with 

the  results  repeated  in  Fig.  16  for  comparison.  The  quenching  rate  due  to  molecular 
nitrogen  is  A:g(jv2)=2.1xl0®(Torr  •  s)~^  or  6.7xl0~^^cm®/s.  Since  this  quenching  rate 
is  smaller  than  kq(^N20)  aiid  tke  density  ratio  [N2]/[N20]  is  very  small,  the  collisional 
quenching  due  to  N2  is  calculated  to  be  negligible  compared  to  that  due  to  N2O. 
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Partial  Pressure  (Torr) 


Figure  16:  Stem-Volmer  plot  of  atomic  nitrogen  (3p)  ^S|y2  state  decay  with  quenching 
by  photolyzed  N2O  products,  N2O,  and  N2. 


The  quenching  by  nitric  oxide  is  also  of  interest,  since  NO  is  formed  by  the  sec¬ 
ondary  process  in  Eq.  5.2  and  is  the  product  of  a  secondary  reaction  involving  N2O 
and  0(^D)[41].  The  quantum  yield  of  NO  during  N2O  photolysis  at  147  nm  has  been 
found  to  be  as  much  as  0.5.  This  value,  though,  includes  contributions  from  the 
secondary  reactions  that  may  not  have  time  to  develop  before  the  <27  ns  lifetime 
TALIF  signal  is  quenched.  With  the  expected  quantum  yield  of  NO  being  <0.5  and 
the  fraction  of  [NjO]  photolyzed  being  <2%,  it  can  be  assumed  that  [N0]/[N20]<1% 
in  the  TALIF  region.  In  an  attempt  to  quantify  the  N  atom  (3p)  ^83/3  quenching  by 
NO,  a  mixture  of  1%  NO  in  N2  was  introduced  downstream  of  a  2  Torr  N2  microwave 
discharge.  Once  introduced  into  the  flowing  afterglow,  there  was  a  high  degree  of 
uncertainty  in  the  actual  density  of  NO  due  to  the  fast  rate  of  reaction  with  ground 
state  N  atoms.  This  is  the  same  reaction  used  for  N  atom  titration,  which  is  given  in 
Eq.  3.30.  It  was  only  possible  to  introduce  a  10  mTorr  partial  pressure  of  NO  before 
the  ground  state  N  atom  density  was  depleted.  No  change  in  the  decay  constant  was 
recorded  beyond  experimental  uncertainty  over  this  TniniTnal  pressure  range. 

Just  as  the  experimental  measurement  of  the  N  (3p)  ^8^/3  quenching  rate  by  NO 
has  proven  to  be  elusive,  an  estimate  of  the  rate  by  theoretical  means  is  likewise  diffi¬ 
cult.  A  quenching  rate  upper  limit  can  be  estimated  by  a  gas  kinetic  rate  calculation 
if  the  effective  collision  diameters  of  each  particle  is  known.  Unfortunately,  this  gas 
kinetic  rate  is  not  well  defined  for  collisions  involving  the  N  (3p)  ^83/3  state  due  to 
uncertainty  in  the  collision  diameter  of  this  atomic  state. 

Although  kfio  has  not  been  determined,  the  measured  value  is  relatively  large 
and  it  would  be  unlikely  to  find  that  k^o  exceeds  kn^o  by  more  than  an  order  of  mag¬ 
nitude.  Therefore,  when  considering  that  [N0]/[N20]<1%  in  the  photolysis/TALIF 
region  it  follows  that  A:iyr2o[N20];^>  /satoINO].  This  indicates  that  NO  contributes  only 
sli^tly  to  N  (3p)  ^S|^3  quenching  in  the  photodissociation  experiment., 

5.3.3  Discussion  of  TALIF  Quenching  by  N2O 

Quenching  of  the  N  atom  (3p)  ^8|/3  state  due  only  to  nitrous  oxide  was  found  to 
have  a  relatively  high  rate  of  A:5(jv2O)=4.2xl0-^°cm®/s,  eclipsing  the  rate  of  molecular 
nitrogen  by  more  than  a  factor  of  six.  The  measured  rate,  though,  is  still  a  factor 
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of  4  smaller  than  the  effective  quendiing  rate  measured  during  photolysis.  Despite 
the  different  quenching  rates  measured  in  the  photolysis  and  post-discharge  environ¬ 
ments,  it  is  proposed  that  N2O  is  the  primary  quencher  of  N  (3p)  ^83/3  in  both  cases. 
The  differences  in  the  quench  rates  between  the  two  environments  results  from  the 
quenching  rate  during  photolysis  being  enhanced  by  excess  kinetic  energy  imparted 
to  the  atomic  nitrogen  during  photodissociation.  If  one  assumes  the  quenching  cross 
section  to  be  constant  as  a  function  of  kinetic  energy,  then  the  factor  of  4  difference  in 
the  N2O  quenching  rates  results  from  a  proportional  difference  in  collision  frequency 
caused  by  a  factor  of  4  difference  in  N  atom  velocity.  With  the  average  translational 
kinetic  energy  of  the  atomic  nitrogen  expressed  as 

<  Ej.  >=  imS*  =|a:T  (5.9) 

it  is  seen  that  a  factor  of  4  difference  in  velocity  corresponds  to  a  factor  of  16  difference 
in  kinetic  energy  between  N  atoms  in  the  photolysis  region  and  the  post-disdiarge. 
The  N  atoms  in  the  flowing  after^ow  of  the  microwave  discharge  are  assumed  to 
be  at  T»300K  with  a  kinetic  energy  of  ~0.039  eV.  A  factor  of  16 x  applied  to  this 
energy  would  give  the  N  atoms  in  the  photolysis  region  a  kinetic  energy  of  0.6  eV. 
This  represents  only  a  portion  of  the  total  translational  energy  disposed  in  the  re¬ 
coil  of  NO(X^n)  and  N(^S)  photofragments.  The  total  translational  energy  can  be 
determined  from  conservation  of  energy  and  momentum  as 

<  >tot=  <  e„  >»,  (5.10) 

rriffo 

which  results  in  <  Etr  >tot=0.9  eV. 

The  total  excess  energy  available  for  distribution  between  the  photofragments 
can  be  determined  from  known  properties.  The  bond  dissociation  energy  that  is 
needed  to  produce  the  N  atom  is  Do(N-NO)=4.99  eV.  The  energy  delivered  by  an 
absorbed  207  nm  photon  is  5.99  eV,  leaving  1.00  eV  to  be  distributed  among  the 
photofragments.  This  allows  the  prospect  of  the  N  atom  fragment  receiving  the 
necessary  <  Etr  >n=0-6  eV  when  90%  of  the  excess  energy  is  channeled  to  fragment 
translational  energy  and  the  remaining  10%  is  partitioned  into  internal  energy  within 
the  NO  fragment.  Although  the  dynamics  of  the  primary  photodissociation  process 
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producing  N2  and  0(^D)  have  received  considerable  attention[42, 43, 44],  no  studies  to 
date  have  measured  the  partitioning  of  energy  in  the  NO  and  N(^S)  photofragments. 

Two  possible  methods  of  directly  measuring  the  partition  of  energy  of  the  NO  and 
N(^S)  photofragments  would  be  to  perform  rotationally  resolved  LIF  on  the  NO  or  to 


measure  the  doppler  width  of  the  TALIF  absorption  transition.  The  LIF  spectrum 
of  NO  would  indicate  the  amount  of  energy  transferred  to  internal  rotational  states 
of  the  molecular  and  the  resulting  kinetic  energy  to  the  NO  and  N  photofragments 
could  be  calculated.  A  more  direct  method  of  determining  the  velocity  of  the  N 
atom  photofragments  would  be  to  measure  the  doppler  broadening  of  the  two-photon 
absorption  line.  The  FWHM  of  the  doppler  broadened  line,  Auj^,  is  given  by  the 
formula[39] 

2^21n(2) 

Ai/d  =  — - Uo\  -nr  5.11) 


=  (5.11) 

where  is  the  peak  absorption  frequency,  M  is  the  atomic  weight,  R  is  the  gas 
constant,  T  is  the  temperature  of  the  absorbing  species,  and  c  is  the  speed  of  light. 
The  doppler  width  can  be  related  to  the  normalized  absorption  lineshape  function, 
for  laser  two-photon  absorption  as 


_  2y/ln{2) 

^/5F(2Al/2  + 


41n(2)(2i/  —  i/oY 
(2Aul  -h  Ai4  -t-  AvS) 


(5.12) 


where  A^t  is  the  FWHM  frequency  spread  of  the  laser  radiation  and  Ai/^  is  the  natu¬ 
ral  linewidth  of  the  two-photon  transition  which  is  usually  comparatively  small.  Thus, 
the  FWHM  of  the  two-photon  absorption  lineshape  that  includes  all  the  linewidth 
parameters  is  (2Ai/^  -I-  Ai/^  +  Ai/^)^/^. 

The  experimental  absorption  lineshape  spectrum  taken  in  the  afterglow  of  the 
microwave  discharge,  shown  in  Fig  1,  Section  3.1,  was  measured  to  have  a  FWHM 
of  1.6  cm~^.  In  the  afterglow,  the  atomic  nitrogen  has  an  approximate  temperature 
of  300  K,  which  would  give  a  doppler  width  of  Ai/d=0.16  cm"^.  This  doppler  width 
is  negligible  in  terms  of  its  contribution  to  the  1.6  cm"^  width  of  the  two-photon 
absorption  lineshape.  Assuming  a  very  small  natural  linewidth,  the  laser  radiation 
spread  dominates  the  line  broadening  effects  and  is  calculated  to  be  Ai>t=1.13  cm“^. 
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The  Ime  broadening  effect  can  be  estimated  for  the  case  where  the  recoiling  atomic 
nitrogen  photofragments  have  a  factor  of  16  x  greater  kinetic  energy  than  the  dis¬ 
charge  produced  atoms.  The  corresponding  velocity  of  a  recoiling  N  atoms  would 
produce  a  maximum  shift  in  the  absorption  frequency  of  '^^0.5  cm“^.  The  doppler 
shift  in  the  photolysis  case  would,  therefore,  cause  lineshape  broadening  of  a  few 
tenths  of  a  wavenumber  greater  than  in  the  discharge  case. 

The  small  predicted  difference  between  the  lineshapes  for  the  different  tempera¬ 
tures  indicates  that  it  would  be  difficult  to  accurately  determine  the  velocity  of  the 
photofragments  from  a  linewidth  measurement  using  this  laser  configuration.  A  much 
better  laser  configuration  for  this  measurement  would  be  one  that  provided  207  nm  ra¬ 
diation  with  a  radiation  spread  of  Ai/i  «0.2  cm~^  Neither  the  doppler  measurement 
nor  the  rotationally  resolved  NO  LIF  experiment  was  attempted  since  considerable 
alterations  in  the  experimental  set-up  would  have  been  necessary.  With  a  further 
investment  of  equipment  and  time,  the  energy  partitioning  between  the  NO(X2n) 
and  N  (3p)  ^S|/2  photofragments  could  be  quantified. 

5.4  Conclusion 

Photodissociation  of  nitrous  oxide  with  subsequent  N  atom  TALIF,  both  using  207 
nm  laser  radiation,  was  accomplished  in  a  low  pressure  N2O  cell.  A  calculated  esti¬ 
mate  showed  the  number  density  of  the  photolyzed  N2O  was  <2%  of  the  total  N2O 
density  and  that  [N0]/[N20]<1%  in  the  photolysis  region.  The  N  atom  (3p)^S3/2 
state  was  determined  to  be  quenched  by  N2O  in  the  photolysis  region  at  a  rate  of 
1.8xl0“®cm®/s.  A  separate  measurement  of  N  atom  (3p)  ^83/2  state  quenching  by 
isolated  N2O  in  the  flowing  afterglow  of  a  microwave  discharge  produced  a  quenching 
rate  of  A:5(Ar2O)=4.2xl0”^®cm^/s.  The  difference  in  the  quenching  rates  is  consistent 
with  '^90%  of  the  excess  photolysis  energy  being  channeled  into  translational  energy 
of  the  NO(X^n)  and  N(^S)  photofragments.  The  high  quenching  rate  of  N  (3p)  ^S|/2 
in  the  photolysis  region  is  consistent,  as  well,  with  quenching  entirely  from  collisions 
with  N2O. 
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CHAPTER  6 


Laser  Studies  of  Atomic  and  Molecular  Nitrogen 
Species  in  High  Frequency  N2  and  N2/H2 

Discharges 


6.1  Introduction 

Discharges  involving  N2  and  H2  gas  mixtures  have  recently  re-emerged  as  a  topic 
of  interest  in  gaseous  electronics  researdi.  The  interest  was  motivated  by  recent 
experimental  and  theoretical  evidence  that  certain  N2/H2  gas  mixtures  in  a  dis¬ 
charge  provide  an  increased  efficiency  in  atomic  hydrogen  or  nitrogen  production 
beyond  that  of  pure  H2  or  N2  discharges[20,  21].  The  potential,  therefore,  exists 
for  the  development  of  more  efficient  hydrogen  and  nitrogen  atomic  sources  by  tai¬ 
loring  the  gas  discharge  chemistry  through  mixtures  in  the  feed  gas.  The  effect  of 
high  fractional  dissociation  of  H2  for  small  amounts  of  H2  in  N2  has  been  clearly 
demonstrated  experimentally [17,  21,  23]  and  strong  theoretical  arguments  have  been 
presented  which  claim  that  heavy  particle  collisions  are  responsible  for  the  high  rate 
of  H2  dissociation  under  these  conditions[45].  Evidence  also  exists  that  dissociation 
of  N2  may  be  enhanced  in  N2/H2  mixtures.  It  has  been  predicted  that  small  per¬ 
centages  of  H2  in  N2  alter  the  heavy  particle  interaction  within  the  discharge  in  such 
a  way  that  there  is  an  enhancement  in  the  high  energy  tail  of  the  electron  energy 
distribution  function  (EEDF),  where  electrons  have  sufficient  energy  to  create  direct 
impact  dissociation [20,  21].  Some  experimental  evidence  exists  for  a  DC  discharge 
that  shows  an  increase  in  N  atom  density  for  small  additions  of  H2  in  N2  in  agreement 
with  theoretical  predictions[23].  The  mean  electron  energy  in  a  DC  discharge  shows 
the  greatest  increase  as  the  H2  concentration  is  varied  from  0-10%[45].  As  the  N2/H2 
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discharge  becomes  further  diluted  with  H2,  the  mean  electron  energy  continues  to 
slowly  rise  due  to  the  declining  influence  of  the  large  inelastic  losses  of  N2[20]. 

It  is  of  practical  interest  to  analyze  the  N  atom  production  within  high  frequency 
discharges  through  N2/H2  mixtures,  since  rf  or  microwave  discharges  would  most 
likely  be  employed  as  atomic  sources  in  industrial  applications.  In  this  chapter,  N2 
and  N2/H2  high  frequency  discharges  were  investigated  with  emphasis  on  the  volume 
processes  resulting  from  heavy  particle  collisions  and  their  relation  to  atomic  nitrogen 
production.  A  low  power  pulsed  rf  parallel  plate  discharge  and  a  flowing  microwave 
discharge  were  analyzed  separately  with  laser  diagnostics  and  optical  emission  spec¬ 
troscopy.  An  advantage  of  the  pulsed  rf  disdiarge  system  was  the  temporal  control 
of  the  electronic  excitation  whidi  provided  the  opportunity  to  observe  species  during 
the  evolution  or  decay  of  the  discharge.  The  parallel  plate  rf  cell  also  allowed  full 
optical  access  of  the  discharge  volume  and  proved  to  be  a  useful  tool  in  monitoring 
certain  effects  of  heavy  particle  interaction.  The  rf  cell,  though,  was  of  a  relatively 
jgTTiflll  scale  and  could  not  produce  the  high  fractional  dissociation  of  N2  that  is  de¬ 
sired  from  an  atomic  nitrogen  source  for  plasma  processing.  E3q)eriments  were  also 
performed,  therefore,  with  a  higher  powered  flowing  gas  microwave  excited  system. 
By  introducing  small  amounts  of  hydrogen  into  the  nitrogen  microwave  discharge 
and  observing  the  resulting  downstream  N  atom  density,  specific  gas  mixtures  were 
determined  to  be  favorable  for  low  level  enhancement  of  atomic  nitrogen  production. 
Emission  spectroscopy  of  the  microwave  discharge  along  with  the  heavy  particle  in¬ 
teraction  data  from  the  rf  disdiarge  cell  provided  input  for  a  kinetics  model  which 
explains  the  observed  behavior  of  the  N  atom  density  resulting  from  the  addition  of 
H2  into  the  N2  microwave  disdiarge. 

The  TALIF(207)  technique  was  used  in  the  pulsed  rf  discharge  to  observe  atomic 
nitrogen,  while  single-photon  LIF  was  applied  to  the  detection  of  N2(A®Sj{’)  and 
N2(B^n5)  molecular  states.  Emission  spectroscopy  was  also  used  to  gain  information 
on  the  molecular  states.  Each  of  the  experimental  techniques  were  described  in  Chap¬ 
ter  3.  These  optical  techniques  were  utilized  to  monitor  temporal  and  spatial  changes 
in  the  concentrations  of  N2(A^E;[^),  N2(B^np),  N2(C®nu),  and  atomic  N(2p®)^S3/2. 
Observations  of  N2(A®E+)  and  N2(B*np)  in  the  post-discharge  along  with  a  measured 
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N  atom  concentration  provided  information  on  the  corresponding  changes  within 
the  vibrational  manifold  of  the  electronic  ground  state  for  N2(X^E+,5<v<13).  Fur¬ 
thermore,  the  absolute  N  atom  concentration,  determined  by  titration,  allowed  an 
absolute  calibration  of  the  density  of  each  of  the  molecular  species  by  coupling  the 
temporal  LIF  and  emission  measurements  with  known  rate  constants  of  several  heavy 
particle  reactions  occurring  within  the  time  period  after  the  discharge  has  been  turned 
off. 

In  the  flowing  microwave  discharge,  the  TALIF(207)  technique  was  employed  to 
monitor  the  nitrogen  concentration  as  up  to  10%  H2  was  added  to  the  N2  feed  gas. 
Emission  spectroscopy  of  the  discharge  volume  and  the  late  afterglow  were  utilized 
to  observe  the  behavior  of  molecular  nitrogen  species.  The  optical  diagnostic  data 
for  this  N2/H2  microwave  discharge  agree  well  with  the  proposed  kinetic  model  which 
predicts  an  elevated  N  atom  density  occurring  at  small  %H2  due  to  an  enhancement 
in  mean  electron  energy.  The  microwave  data  also  indicate  an  eventual  decline  in  N 
atom  density  as  the  %H2  is  increased  above  1%.  It  is  suggested  in  this  work,  that  dis¬ 
sociation  channels  involving  vibrationally  excited  ground  state  N2(X^S+,v')  molecules 
are  influential  within  the  microwave  experiment,  especially  at  the  higher  pressures. 
This  assertion  allows  the  eventual  decline  in  N  atom  density  to  be  attributed  to  the 
depletion  of  N2(X^S+)  high  v'  states  and  the  resulting  decrease  in  N  atom  production 
from  the  V-V,  V-T  and  e-V  N2  dissociation  channels. 

6.2  Theory 

6.2.1  Rate  equations  for  N2/H2  Discharge  Kinetics 

Several  sets  of  rate  equations  are  presented  in  Tables  3  through  6  in  order  to  follow 
the  behavior  of  the  discharge  species  of  interest.  The  rates  given  in  the  tables  have 
the  unit  of  cm^/s  where  the  per  molecule  dimension  is  implied.  This  is,  by  no  means, 
a  comprehensive  list  of  reactions  within  the  discharge,  but  a  list  that  is  representative 
of  the  discharge  kinetics  that  are  most  influential  to  the  changing  atomic  flux  as  a 
function  of  gas  mixture.  Kinetic  equations  for  the  relatively  short  lived  N^iB^Ug) 
and  N2(C®nu)  states  are  also  given  in  Table  6.  This  is  not  because  of  the  influence  of 
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these  moleculax  states  on  the  discharge  kinetics,  but  because  these  states  have  been 
observed  by  LIF  and  emission  in  this  experiment  and  are  good  indicators  of  changes  in 
certain  discharge  properties,  such  as  mean  electron  energy  and  N2(X^S+,v)  behavior. 


Process 

Rate  (cm®/s) 

Ref 

(a) 

e  +N2(X1S+,v=0)  N  +  N  +  e 

3x10-^®  * 

[46] 

(b) 

e  +N2(XiS+,v=1...45)  N  +  N  +  e 

1.9xl0-“* 

[46] 

(cl) 

N2(X,v=45)+N2(X,v^)  — >  N  +  N  +  N2(X,v^-1) 

[47] 

(c2) 

N2(X,v=45)+N2  ->  N  +  N  +  N2 

4x10-12* 

[46] 

(d) 

N  +  N  +  N2(XiE+)  ^  N2+  N2(B3na,v=10,ll,12) 

_ .02  • j.-. 1 r m /XT  Cv>in— 16 

4.4x  10“®® 

\r  _ _ 2  rom 

[48] 

♦Rate  coefiB-dents  csJculated  for  E/N  =  5x10  V  cm^  [20] 


Table  3:  CollisionaJ  N2  Dissociation  and  N  atom  Destruction  Processes. 


Process 

Rate  (cm®/s) 

Ref 

(a) 

N(^5)  +  N2(A®E+)  ^  N(2P)  +  N2(X'S+,v) 

5.0x10"^^ 

[49] 

(b) 

H2  +  N2(a'  1E-)  ->  H  +  H  +  N2(X1E+) 

2.6x10-11 

[50] 

(c) 

H2  +  N2(A3E+)  H  +  H  +  N2(X1E+) 

2.4x10-1® 

[51] 

(d) 

H  +  N2(Aas+)  -»  H  +  N2(X‘E+) 

2.1x10-1® 

[52] 

(e) 

H  +  N2(a'  ‘E;)  -»  H  +  N2(X'E+) 

1.5x10-11 

[21] 

Table  4:  Processes  for  Quenching  of  N2  Metastables  and  a'  . 
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(a)  N2(a^  S^^)+N2(a^  N2+  N2+  6 

(b)  N2(A3E+)+N2(a'  ^E")  N++  N2+  e 

(c)  e  +N2(X1E+ v=0)  ^  Nt+  e 

(d)  e  +N2(A"E+)  ^  N++  e 


2.0x10"^“ 
5.0x10-“ 
2.2x10-1®  t 
2.5xl0-“f 


fRate  coejB5.cients  calculated  for  E/N  =  6xl0-i®  V  cm^  [46] 
Table  5:  Collisional  Ionization  Processes. 


Process 


(a)  N2(A®E+)  +  N2(X1E+ v>5) 

-4  N2(B®n5)+  N2(X1E+) 

(b)  e  +N2(X1E+ v=0)  -4  N2(B®ng)+  e' 

(c)  N2(A*E+)  +  N,(A’S+)  ->  N2(B^n,)+ N2(X1E+) 

(d)  NzjB’ns)  +  N2(X‘E+)  ^  N2(A*E+)+  N2(X‘E+) 

(e)  N2(BSn,)  +  N2(X‘E+)  ->  2N2(X'E+) 

(f)  N2{B>n,)  +  Hj  N2(A=E+)+  Hj 

(g)  NjCB’n,)  N2(A’e+)  +  h./ 

(h)  e  +N2(X'E+,v=0)  N2(C=II„)+  e 

(i)  N2«?n.)  ->  Ni(B°n.)  +  hv _ 

fRate  coefficients  calculated  for  E/N  =  6xl0-i® 


Rate  (cni®/s) 


3.0x10-“ 
1.2xl0-i2f 
2x10-1" 
3.0x10-11 
2x10-12 
5.0x10-11 
tb=5.0x10‘ 
3.3xl0-“t 
rc=4.0xl0' 
V  cm^  [46] 


[56] 
[54] 

[57] 

[58] 
[21] 
[58] 

•«s  [59] 
[54] 
®s  [59] 


Table  6:  Collisional  Processes  for  Excitation  and  Decay  of  N2,  B^IIo  and  C®II«. 
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It  should  be  noted  that  the  rate  of  the  process  in  Table  6(c)  includes  the  contri¬ 
bution  from  the  intermediate  reaction  step  [60,  61] 

->  N2(X)  +  N2(X)  +  +  lu'  (6-1) 

where  the  N2(C®Ilu)  that  is  formed  by  the  pooling  reaction  of  N2(A^S;|')  will  quickly 
radiate  to  the  N2(B*nj)  state.  Processes  marked  with  *  or  f  have  rates  that  are 
dependent  on  the  electron  energy  distribution  function  and  may  vary  significantly  with 
discharge  condition.  In  these  cases,  rate  coefficients  have  been  taken  from  fLoureiro 
et  al.  [46]  or  *Garscadden  et  o/.[20],  where  the  rate  coefficients  were  calculated  for 
an  E/N  of  6x10“^®  V  cm^  and  5x10“^®  V  cm^  respectively,  which  are  of  the  same 
magnitude  as  the  effective  E/N  in  this  work. 

6.2.2  Kinetic  theory  in  the  pulsed  rf  post-discharge 

At  the  termination  of  each  rf  pulse,  the  electrons  and  ions  recombine  within  the  first 
few  microseconds.  The  short-lived  neutral  species  that  were  produced  within  the 
discharge  will,  in  general,  quickly  decay  within  their  respective  coUisional  or  radiative 
time  periods.  The  metastable  species  formed  in  the  nitrogen  discharge,  such  as  atomic 
nitrogen,  N2(X^S+,v),  and  N2(A®S+),  decay  much  slower.  The  collisional  behavior 
of  these  metastables  after  the  discharge  has  turned  off,  especially  at  higher  pressures, 
is  also  known  to  regenerate  some  of  the  short-lived  species,  sudi  as  N2(B®ng).  The 
“post-discharge”  period  of  interest  for  this  work  will  be  defined  as  a  time  period 
starting  at  <^0.2  ms  after  the  rf  discharge  pulse  has  been  turned  off.  This  ensures 
that  the  electrons  and  ions  have  had  ample  time  to  recombine,  thus  simplifying  the 
kinetics  within  the  cell.  Measurements  and  calculations  have  been  made  during  the 
post-discharge  period  out  to  30  ms  after  the  disdiarge  was  turned  off. 

As  an  example  of  the  production  of  N2(B®nj)  in  the  post-discharge  by  collisional 
interaction  of  metastables.  Fig.  17  shows  the  emission  from  the  N2(B®ng,v=4)  state 
after  a  20  ms  rf  pulse  in  5  Torr  of  pure  N2.  The  emission  as  a  function  of  time  indicates 
a  rapid  drop  in  N2(B)  population  with  the  discharge  turn-off,  but  the  N2(B)  decay 
rate  diminishes  within  tens  of  microseconds  after  the  rf  discharge  turn-off  due  to  ex¬ 
citation  from  metastable  heavy  particle  collisions.  The  decay  of  the  remaining  N2(B) 
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population  in  the  post-discharge  is  then  dependent  on  the  decay  of  the  metastables. 
The  emission  from  the  N2(B)  state,  which  is  readily  observed,  can  therefore  be  a 
valuable  indicator  of  the  population  and  behavior  of  the  metastables.  Information  on 
the  behavior  of  the  vibrationally  excited  N2(X^E^,v)  species  is  of  particular  interest, 
since  direct  optical  detection  of  this  species  is  especially  difficult. 

A  kinetic  model  describing  the  temporal  behavior  of  the  metastable  species  during 
the  post-discharge  was  constructed  from  the  heavy  particle  reactions  in  Tables  3 
through  6.  Differential  rate  equations  involving  six  species  are  expressed  in  Eqs.  6.5 
through  6.10  which  include  the  primary  reactions  during  the  post  discharge  period 
in  the  absence  of  electron  impact  processes.  The  six  species  of  discharge  products 
include  atomic  N,  N2(A3E+),  N2(B3np),  N2(B3n5,v=10,ll,12),  N2(XiS+,v>5),  and 
N2(X^S+,v>12). 

The  N2(B,v=10,11,12)  population  is  a  subset  of  the  N2(B)  state  population  in 
the  model.  The  N2(B,v=10,ll,12)  states  are  analyzed  specifically  since  they  are 
selectively  poplulated  by  N  atom  recombination  and  have  been  observed  in  emission. 

The  species  N2(X^S^,v>5)  and  N2(X^S^,v>12)  are  subsets  of  the  full  vibrational 
manifold  of  N2(X^E+).  The  population  of  N2(X^E+,v>5)  is  significant  for  its  role  in 
exciting  N2(B)  states  in  process  6(a).  The  population  of  N2(X^E+,v>12)  was  included 
in  the  model  since  the  reaction 

N2(A)  -h  N2(X,v  >  12)  ^  N2(B,  V  =  10, 11, 12)  +  N2(X)  (6.2) 

is  another  significant  source  of  N2(B,v=10,ll,12)  production  other  than  the  N  atom 
recombination  process  of  Table  3(d). 

The  nonequilibrium  vibrational  kinetics  of  the  N2(X)  state  are  important  in  de¬ 
termining  the  behavior  of  the  N2(X^E+,v>5)  and  N2(X^E+,v>12)  states.  A  general 
characteristic  of  a  system  of  diatomic  molecules,  such  as  N2,  is  the  dominant  V-V 
exchange  in  low-lying  vibrational  levels  due  to  anharmonicity  of  the  molecule  as  an 
oscillator,  which  leads  to  a  V-V  up-pumping  mechanism[62].  For  example  in  molec¬ 
ular  nitrogen,  the  V-V  energy  exchange  with  v>v' 

N2(X,  v'  =  1)  +  N2(X,  v)  #  N2(X,  v'  =  0)  +  N2(X, y  +  1)  (6.3) 
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N2(B  -  A)  Emission  Intensity  (arb) 


Discharge  On 


Figure  17:  N2(B®n5->A^E;J')  emission  intensity  during  turn-off  and  early  post¬ 
discharge  period  of  a  pulsed  rf  nitrogen  discharge  with  a  20  ms  pulse  duration  at 
5  Torr 


displays  a  forward  reaction  rate  that  is  much  greater  than  the  reverse  process  rate. 
The  vibrational  energy  that  has  been  deposited  in  the  lower  N2(X,v)  levels  by  elec¬ 
tron  impact  excitation  is  then  transferred  to  the  higher  v  levels  through  V-V  energy 
exchange.  The  result  is  that  instead  of  a  Boltzmann  vibrational  distribution,  the 
V-V  up-pumping  creates  a  quasi-stationary  Treanor’s  distribution  with  a  plateau  in 
the  mid  vibrational  levels[62],  which  includes  the  N2(X^E+,v>5)  and  N2(X^E+,v>12) 
states.  An  example  of  such  a  distribution  can  be  found  in  Section  6.4.5  in  Fig.  34. 
The  N2(X)  V-V  energy  exchange  rates  at  5  Torr,  with  a  characteristic  time  of  less  than 
a  millisecond[62],  are  much  faster  than  the  N2(X,v)  collisional  reaction  rates  and  wall 
loss  rate  used  in  this  model.  Therefore,  as  vibrational  energy  is  lost  from  the  N2(X,v) 
manifold  through  heavy  particle  interactions  or  de-activation  by  wall  collisions,  the 
vibrational  distribution  is  adjusted  by  V-V  exchange  and  maintains  the  Treanor’s 
distribution  and  plateau.  The  magnitude  of  the  population  in  the  plateau  region 
decreases  according  to  the  loss  of  vibrational  energy.  Although  the  model  does  not 
account  for  V-V  energy  exchange,  the  differential  rate  equations  for  [N2(X^E+,v>5)] 
and  [N2(X^Z1+,v>12)]  include  loss  terms  that  represent  the  resulting  decrease  in  the 
hi^er  level  vibrational  states  due  to  the  reactions  and  wall  de-activation. 

Each  of  the  differential  rate  equations  for  the  metastable  species  include  a  term 
that  represents  a  rate  of  change  in  density  of  that  species  due  to  spatial  diffusion. 
The  diffusion  loss  in  the  rf  cell  is  highly  dependent  on  the  resulting  loss  rate  of 
the  diffusing  species  at  the  electrode  surface.  In  each  case  the  term  is  expressed 
in  the  form  -DjV(f/A|^[M],  where  M  is  the  species,  Djv/  is  the  diffusion  coefficient  of 
the  species  through  N2,  which  is  pressure  dependent,  and  Aat  is  the  effective  diffusion 
length,  which  depends  on  the  reactor  dimensions  and  the  surface  recombination  rates. 
The  specific  values  of  this  parameter  for  each  metastable  species  in  the  model  have 
been  determined  by  various  methods.  For  atomic  N,  a  value  of  £>i\r/A^=24  s~^  was 
determined  empirically  from  N  atom  temporal  decay  data  in  the  5  Torr  discharge.  A 
value  of  Dn^{x,v)/An2(x,v)='>^'^  for  N2(X,v>5)  and  N2(X,v>12)  in  the  rf  cell  at  5 
Torr  was  determined  from  a  theoretical  derivation  presented  in  Section  7.4.2.  A  value 
for  N2(A)  diffusion  loss  at  5  Torr  of  Dn^^a)/ was  found  by  modeling 
spatial  N2(A)  LIF  profiles  in  Section  6.3.4.  Diffusion  loss  for  the  N2(B)  state  was 
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not  includsd  in  the  model  because  of  the  very  short  radiative  lifetime  of  the  state 
compared  to  the  possible  diffusion  times. 

The  first  four  differential  rate  equations  presented  in  Eqs.  6.5  through  6.8  corre¬ 
spond  to  the  species  N,  NaCA^SJ),  and  Na(B3ng,v=10, 11,12)  which  have 

all  been  observed  by  optical  diagnostics  in  the  rf  discharge  cell.  The  final  two  rate 
equations  represent  N3(X^E^,v>5)  and  N3(X^S^,v>12)  in  Eqs.  6.9  and  6.10,  which 
have  not  been  monitored  in  the  discharge.  Absolute  densities  for  these  two  species 
may  be  determined  from  a  fit  of  all  the  equations  to  the  experimental  observations  of 
the  other  species. 

The  processes  involving  creation  or  destruction  of  atomic  nitrogen  are  listed  in 
Table  3(a)-(d).  The  direct  electron  impact  dissociation  processes  3(a)  and  3(b)  are 
not  included  in  the  rate  equation  for  the  post-discharge  since  the  electron  density 
becomes  negligible  within  the  first  few  microseconds  after  the  discharge  is  turned  oflP. 
The  dissociation  processes  3  (cl)  and  3(c2)  that  involve  N3(X,v=45)  are  excluded 
as  well  firom  the  post-discharge  rate  equation  since  these  mechanisms  depend  heav¬ 
ily  on  electron  impact  excitation  of  the  lower  v  states.  The  rate  of  change  in  the 
N3(X,v=45)  population  in  the  post-discharge  can  be  determined  from  V-V  and  V-T 
energy  exchange  rates  calculated  firom  SSH  theory[46].  The  calculations  show  that 
the  V-T  reaction 


N2{X,  v=A5)  +  N2^  N2{X,  v-U)  +  N2  (6.4) 


has  a  rate  of  2xl0“^^cm®/s  which  greatly  exceeds  the  reverse  V-T  process  and  the  V- 
V  energy  exchange  rates  involving  N3(X,v=45).  The  high  V-T  rate  of  de-activation 
is  a  general  property  for  Na(X,v>35)[20,  46].  In  the  absence  of  electron  impact 
excitation  of  the  lower  vibrational  levels,  the  Na(X,v=45)  population  will  experience 
a  decay  from  reaction  6.4  with  a  time  constant  on  the  order  of  microseconds.  This 
will  result  in  a  negligible  contribution  firom  processes  3 (cl)  and  3(c2)  during  the 
post-discharge  period  of  interest.  The  post-discharge  evaluation  of  atomic  nitrogen, 
therefore,  depends  only  on  N  atom  diffusion  and  process  3(d),  leading  to  the  rate 
equation 


(6.5) 
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The  rate  equation  for  the  N2(A®E+)  state  in  the  post-discharge,  which  includes 
all  the  production  and  destruction  processes  in  Tables  4  and  6  is 

^N2(A)  Ta 

+i[W2(B)]  +  (6.6) 

*  B 

Process  5(d)  involving  N2(A)  reaction  with  N2(a')  has  a  very  low  rate  in  the  post- 
disdiarge  and  is  not  included  in  Eq.  6.6. 

Similarly  including  all  the  processes  involving  N2(B®ng)  in  Tables  3  and  6  leads 
to  the  following  differential  rate  equation  for  the  post-discharge 
dfiVofB)! 

•  =  hia)[N2iX,v  >  5)][iV2(A)]  -F  A:6(c)[iV2(A)][Ar2(A)] 

-(*:6<4) +fc6(.))[JV2m][JV2(B)]  -  i[iV,(B)| 

+*3M[JVlliV][W2(Jf)]  *  (6.7) 

where  the  last  term  represents  selective  production  of  N2(B,v=10,ll,12)  from  three- 
body  N  atom  recombination. 

The  rate  equation  for  N2(B®nj,v=10,ll,12)  is  expressed  here  using  the  designation 
N2(B*)=N2(B®ng,v=10,ll,12).  The  rate  of  the  pooling  reaction  of  N2(A)  in  process 
6(c)  that  leads  specifically  to  the  N2(B*)  states  has  not  been  measured  but  is  assumed 
to  be  small  and  is  not  included  in  the  rate  equation.  The  rate  equation  for  N2(B*)  is 
then 

dfJV2(S*)l 

— =  ^6(a)[iV2(A,U>12)][iV2(A)]-KA:3(d)[iV][iV][iV2] 

-(isM  +  *6w)|iVj(X)][iV2(S-)]  -  -(JV2(B')].  (6.8) 

Tb 

The  post-disdiarge  rate  equations  for  the  N2(X,v)  states  each  include  two  loss 
terms.  The  kinetic  model  will  show  that  the  N2(X,v)  states  were  most  infiuential 
during  the  first  2  ms  after  the  discharge  was  turned  off  and  the  [N2(X,v)]  populations 
decayed  very  little  during  this  time  period.  At  later  times,  the  N2(X,v)  state  had  a 
negligible  effect  in  the  modeling  of  the  observed  species.  Therefore,  the  accuracy  of 
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the  decay  rate  of  the  N2(X,v)  state  in  the  model  with  its  lack  of  accounting  for  V-V 
interaction  was  not  a  concern.  The  post-discharge  concentration  of  N2(X^S+,v>5) 
is  governed  by  diffusion  loss  and  the  process  in  Table  6(a),  which  results  in  a  rate 
equation  of 

d[N^(X  v  >  5)1  ^  >  5)1  -  >  5)].  (6.9) 

^N2iX,v) 

The  post-discharge  rate  equation  for  N2(X^S+,v>12)  is 

djNijX  v  >12)]  ^  „  >  12)]  -  t6(.)[«2(^)l[Ar2(jr,,j  >  12)].  (6.10) 

^N2iv) 

The  six  differential  rate  equations  have  been  modeled  numerically  using  Euler’s 
method[63]  with  the  initial  densities  of  each  species  used  as  the  boundary  conditions. 
The  initial  densities  of  each  species  were  adjusted  so  that  the  resulting  temporal 
behavior  of  the  densities  fit  the  LIF  and  emission  diagnostic  data  collected  for  several 
of  the  species.  The  results  of  the  numerical  modeling  are  reported  in  Section  6.4. 

6.3  Experimental 

Two  high  frequency  discharge  sources,  a  pulsed  rf  parallel  plate  disdbarge  and  a  higher 
powered  flowing  gas  microwave  discharge,  both  shown  in  Fig.  18,  were  investigated 
separately  in  this  chapter. 

6.3.1  Pulsed  rf  discharge  cell 

The  rf  disdiarge  pulse  was  synchronized  with  the  pulsed  laser  probe  to  study  the 
temporal  bdiavior  of  various  nitrogen  species  within  the  discharge-on  period  and 
shortly  after  the  discharge  pulse  was  turned  off.  The  experimental  details  of  the 
rf  parallel  plate  cell  can  be  found  in  Section  3.3.  When  observing  the  rf  parallel 
plate  discharge  in  operation,  two  bright  disk  shaped  glow  regions  were  evident  within 
the  discharge  volume,  each  parallel  to  an  electrode.  These  regions  corresponded  to 
the  edge  of  the  modulated  discharge  sheath  near  each  electrode  where  there  was  a 
high  rate  of  excitation  by  electron  impact.  The  glow  regions  moved  further  from  the 
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Pulsed  rf 


(b) 


Figure  18:  Sketches  of  (a)  rf  parallel  plate  discharge  cell  and  (b)  flowing  gas  microwave 
discharge  cell 
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electrodes  and  were  more  expansive  for  lower  pressures.  At  1  Torr,  the  two  regions 
effectively  merged  as  one  glow  region. 

The  rf  signal  at  10  MHz  was  amplified  for  pulse  widths  varying  firom  50  /xs  to 
20  ms  with  a  repetition  rate  of  10  Hz.  The  parameters  of  the  pulsed  rf  discharge 
for  gas  pressmres  of  1,3,  and  5  Torr  were  set  and  maintained  so  that  breakdown  was 
not  observed  outside  of  the  interelectrode  space.  The  resulting  power  absorbed  in 
the  N2  discharge  at  5  Torr  was  8  W  for  continuous  rf  excitation.  The  power  was 
measured  by  capturing  voltage  and  current  waveforms  to  the  cell  and  analyzing  them 
to  account  for  the  current-voltage  phase  shift.  The  average  input  power  scaled  down 
as  the  10  Hz  pulse  width  and  duty  cycle  were  lowered.  The  resulting  input  power 
maintained  for  3  and  1  Torr  without  excess  breakdown  was  6  W  and  4  W  respectively. 
For  experiments  involving  a  changing  composition  of  N2/H2  gas  mixture  in  the  rf  cell, 
a  rbanging  discharge  impedance  was  known  to  somewhat  alter  the  rf  power  input, 
but  the  task  of  ensuring  a  constant  rf  input  power  for  each  gas  mixture  was  not 
attempted.  The  data  displayed  in  Section  6.4  shows  effects  from  the  changing  N2/H2 
gas  mixture  that  was  clearly  attributable  to  significant  changes  in  the  heavy  particle 
kinetics  and  not  due  to  slight  rf  power  variations.  A  slow  gas  flow  of  10  seem  per 
Torr  was  used  to  ensure  that  the  measured  loss  rate  of  the  discharge  products  was 
not  effected  by  bulk  gas  transport  from  the  cell. 

The  density  of  atomic  nitrogen  was  measured  in  the  rf  cell  by  the  TALIF(207) 
technique  described  in  Chapter  3.  The  TALIF  signal  intensity  was  calibrated  to 
the  absolute  N  atom  density  by  titration  of  the  atomic  nitrogen  by  NO  added  in  a 
separate  microwave  system.  The  spatial  resolution  of  the  TALIF  measurement  was 
determined  by  the  focal  volume  of  the  laser.  The  200  /xm  resolution  in  the  z  direction, 
which  designates  direction  parallel  to  the  interelectrode  axis,  was  most  important  in 
the  parallel  plate  cell  since  the  gradients  in  the  density  of  discharge  products  was 
typically  much  greater  in  the  z  direction  as  compared  to  the  radial  direction. 

The  N2(A®E+,  v=0)  and  N2(B®n5,v=l)  molecular  states  were  probed  by  separate 
single  photon  LIF  techniques,  also  detailed  in  Chapter  3.  In  eadb  LIF  technique,  the 
laser  was  collimated  to  approximately  1  mm.  This  gave  an  adequate  1  mm  spatial 
resolution  for  the  LIF  probes  in  the  z  direction,  with  the  resolution  in  the  less  critical 
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radial  direction  determined  by  the  field  of  view  of  the  fluorescence  collection  optics. 

6.3.2  Corrections  for  collisional  quenching  of  laser  excited 
states 

When  comparing  relative  LIF  signals  from  the  same  state  but  at  diflferent  background 
pressures,  the  quenching  of  the  ambient  gas  on  the  laser  excited  state  must  be  consid¬ 
ered.  The  quenching  must  be  accounted  for  to  accurately  relate  the  signal  intensity 
to  the  density  of  the  species  being  probed.  The  quenching  rate  by  N2  was  found 
in  Chapter  4  to  be  very  small  on  the  upper  N  (3p)  ^83/3  state  in  the  TALIF(207j 
technique.  For  the  N2(B^n5,v=l)  LIF  probe,  the  measured  fluorescence  decay  times 
of  the  N2(C®nu,v=0)  upper  state  did  not  change  noticeably  over  the  pressure  range 
of  interest.  Signal  corrections  due  to  quenching  were  therefore  not  needed  for  the 
TALIF(207)  or  the  N2(B®ng,v=l)  LIF  probe.  The  LIF  measurement  of  the  N2(A®S+, 
v=0)  state,  on  the  other  hand,  displayed  quite  noticeable  changes  in  fluorescence  de¬ 
cay  tunes  as  a  function  of  pressure  as  shown  in  Fig.  19  (a).  This  was  due  to  quenching 
of  the  upper  N2(B®n5,v=3)  state  by  N2,  which  was  determined  from  the  Stem-Volmer 
plot  in  Fig.  19(b)  to  have  a  rate  of  (4  ±  l)xlO-ii  cmVs,  which  agrees  with,  within 
experimental  uncertainty,  the  rate  measured  by  Piper[58].  This  quenching  rate  has 
been  used  to  correct  the  relative  signal  intensities  of  the  data  taken  at  each  of  the 
three  gas  pressures. 

6.3.3  Comparison  of  the  theoretical  sensitivity  of  the  of  N2(A) 
and  N2(B)  LIF  probes 

Since  LIF  techniques  which  probe  the  N2(A)  and  N2(B)  states  were  be  used  together 
in  this  study,  a  theoretical  analysis  is  presented  here  which  compares  the  fluorescence 
output  and  sensitivity  of  each  technique.  The  rate  of  change  of  population  of  an 
upper  N„  excited  state  due  to  laser  absorption  from  a  N„i  state  is  given  by[64]  to  be 

_  NmAmnP{v)  /■«  1 1  \ 

dt  87r/ii/3 

where  A„i„  is  the  Einstein  probability  of  spontaneous  emission,  u  is  the  frequency 
of  laser  radiation,  and  p{u)  is  the  spectral  energy  density  of  the  laser  radiation. 
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Figure  19:  (a)  Time  resolved  signal  for  LIF  probe  of  N2(A®E+,  v=0)  state  showing 
fluorescence  quenching  efiect  (b)  Stem-Volmer  plot  of  laser  excited  N2(B®nj,v=3) 
state  quenching  by  N2 


Assuming  pulsed  laser  radiation  and  integrating  over  the  time  of  a  single  LIF  pulse, 
the  number  of  fluorescence  photons  collected,  S^,  from  this  absorption  rate  is  derived 
as 


Sffi  oc 


mn 


(6.12) 


where  !„  is  the  radiation  energy  of  a  laser  pulse  and  is  the  fluorescence  quantum 
yield,  where 


Qmn  — 


■^mn 


(6.13) 


Amn  +  Km' 

The  fluorescence  quantum  yield  accounts  for  the  quenching  of  the  laser  excited  state 
by  collisions  with  N2  before  spontaneous  emission  can  occur. 


In  comparing  the  speciflc  LIF  probes,  the  m  and  n  subscripts  in  Eq.  6.12  denoting 
the  lower  and  upper  states  are  replaced  by  A,  B,  and  C,  which  represent  the  N2(A), 
N2(B)  and  N2(C)  states  respectively.  Recall  that  the  N2(A)  LIF  probe  involves  the 
N2(A)— >N2(B)  transition  while  the  N2(B)  LIF  probe  uses  the  N2(B)^N2(C)  transi¬ 
tion.  A  ratio  of  the  N2(B)  LIF  signal  to  the  N2(A)  LIF  signal  can  be  expressed  from 
Eq.  6.12  as 


Sb  _  Nb  Aqb  •fyg  QcB  ..  .V 

Sa~  Na  Aba  Iu^v%Qba  ^  ^  ^ 

A  good  approximation  is  that  Acb/Aba^tbI'^c  so  by  using  the  values  in  Table  6,  a 
ratio  of  AcB/AjByi=125  is  derived.  Also  from  the  information  in  Table  6  and  Eq.  6.13, 
it  is  found  for  5  Torr  of  N2  that  Qcn=0.92  while  Qb^=0.04.  The  laser  wavelengths  of 
688  nm  for  the  N2(A)  probe  and  358  nm  for  the  N2(B)  probe  results  in  v\/u^={^.lAl 
and  for  this  laser  source  it  is  common  that  Ii/b/I».^=10%.  Using  these  values,  a 
comparison  in  LIF  signal  sensitivity  gives 


Sa 


40x 


IVb 

Na 


(6.15) 


This  indicates  a  more  sensitive  LIF  probe  for  N2(B)  compared  to  N2(A)  detection. 
Under  most  discharge  conditions,  though,  the  N2(A)  population  will  be  much  greater 
than  the  N2(B).  In  this  study  it  was  common  that  the  N2(A)  population  was  100 
to  500  times  greater  than  N2(B)  which  would  result  in  a  Ss/S^  of  ~10  to  40%.  In 
summary,  the  sensitivities  of  the  N2(A)  and  N2(B)  LIF  probes  have  been  calculated 
to  be  relatively  similar  so  that  one  technique  will  not  severely  limit  the  range  of 
conditions  that  tandem  measurements  can  be  made. 
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6.3.4  Spatial  dependence  of  diffusion  in  rf  cell 

The  N2(A)  production  in  the  pulsed  rf  cell  was  not  homogeneous  but  was  found  to  be 
peaked  in  the  sheath  regions  as  shown,  for  example  at  5  Torr,  in  the  interelectrode 
spatial  distribution  plot  of  N2(A)  LIF  in  Fig.  20.  This  non-uniform  distribution 
caused  the  diffusion  loss  rate  of  the  N2(A®i:+,  v=0)  state  to  vary  with  position  and, 
in  many  cases,  to  vary  with  time  after  the  rf  pulse.  With  several  other  significant 
loss  terms  in  the  N2(A)  differential  rate  equation,  Eq.  6.6,  it  was  not  possible  to 
experimentally  isolate  and  study  the  N2(A)  loss  by  diffusion.  A  numerical  calculation 
was,  therefore,  done  to  find  a  spatial  coordinate  within  the  rf  cell  where  a  consistent 
rate  of  diffusion  loss  to  the  walls  occurred  over  the  time  period  of  interest.  The 
calculation  also  determined  the  value  of  DAr2(A)/A^2(^j  at  this  spatial  coordinate  for 
several  conditions  to  be  assigned  to  Eq.  6.6  in  the  kinetic  model.  The  decay  of  N2(A) 
was  modeled  using  the  two  dimensional  numerical  model  for  this  parallel  plate  reactor 
that  is  discussed  in  detail  in  Appendix  A.  The  numerical  model  used  the  experimental 
N2(A)  LIF  spatial  distributions  as  inputs  and  modeled  the  N2(A)  decay  considering 
only  the  diffusion  wall  loss  term  in  Eq.  6.6.  The  N2(A)  was  assumed  to  be  lost  at 
the  wall  surface  with  unit  probability  [65].  It  was  found  that  the  spatial  region  3  mm 
from  either  electrode,  as  shown  in  Fig.  20,  produced  a  diffusion  loss  rate  that  was 
most  consistent  over  the  time  period  0.2  ms  to  2  ms  after  the  rf  pulse  for  each  of 
the  initial  spatial  distributions.  The  LIF  probe  measurements  for  N2(A)  decay  were, 
therefore,  made  at  this  position,  3  mm  from  the  electrode.  The  measured  loss  rate 
from  the  model  was  then  assigned  as  the  value  of  ^N2{A)/^%2{a)  condition  in 

the  kinetic  model  that  was  used  to  fit  the  experimental  decay  data.  The  loss  rates 
due  to  diffusion  determined  by  the  model  are  shown  in  Table  7. 

6.3.5  Afterglow  emission  produced  by  N  atom  recombination 

Afterglow  emission  measurements  which  indicate  the  density  of  the  N2(B^nj,v=10- 
12)  states  were  also  made  since  these  states  are  greatly  influenced  by  N  atom  re¬ 
combination  in  the  post-discharge.  The  three-body  atomic  recombination  process. 
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Figure  20:  Spatial  distribution  of  N2(A)  relative  density  between  electrodes  in  the 
pulsed  rf  cell  at  1  and  5  Torr  with  a  5  ms  pulse  duration 


Pressure 

Diffusion  Coefficient 

Diffusion  Loss  Rate,  D/A^ 
Determined  from  Model 

1  Ton- 

170  cm^/s 

600  s"^ 

3  Torr 

57  cm^/s 

310  s-i 

5Torr 

34  cm^/s 

220  s-i 

Table  7:  N2(A)  Diffusion  Loss  Rates  at  3  mm  from  Electrode. 


listed  in  Table  3(d),  produces  N2(B®ng)  selectively  in  the  v=10,ll,  and  12  vibra¬ 
tional  states  [48].  A  spectrum  of  the  emission  20  ms  after  the  discharge  was  turned 
off  in  the  rf  pulsed  cell  is  shown  in  Fig.  21. 

The  emission  was  very  weak  from  these  bands  due  to  the  relatively  low  N  atom 
density  and  the  [N]^  dependency  of  the  emission.  For  each  measured  wavelength  in 
Fig.  21,  the  signal  through  the  monochrometer  was  recorded  by  photon-counting  over 
several  minutes  to  reduce  the  noise  level.  A  much  greater  signal  was  produced  when  an 
interference  filter  was  placed  directly  in-line  with  the  PMT.  Figure  21  also  shows  the 
afterglow  emission  through  a  5  nm  bandwidth  interference  filter  centered  at  585  nm. 
The  bandpass  of  the  filter  coincides  closely  with  the  N2(B,v=10  ->  A,v=6)  emission. 
This  filter  was,  therefore,  used  to  measure  the  temporal  decay  of  the  emission  from 
the  N2(B,v=10)  state  in  the  post-discharge. 

The  total  steady  state  population  of  N2(B,v=10,ll,12)  in  the  post-discharge  is  rep¬ 
resented  by  the  differential  rate  equation  in  Eq.  6.8,  so  information  on  the  fractional 
population  of  the  measured  v=10  state  compared  to  the  total  [N2(B,v=10,ll,12)]  is 
required.  The  state  specific  emission  intensity,  I(v',v"),  can  be  given  by[66] 

1(^,0  =  [N2(B,v')]Ay,v,^  (6.16) 

where  Ayy/  is  the  Einstein  A  coefficient  and  Avy'  is  the  wavelength  of  the  N2(B,v'->^A,v") 
transition.  The  wavelengths  of  the  transitions,  A^.v",  are  roughly  the  same  so  that 
[N2(B,v’)]  can  be  approximated  as 

lN2(B,v')](xi^  (6.17) 

Ay>y» 
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N2(B-A)  Late  Afterglow  Emission  (Normalized) 


Figure  21:  Av=4  emission  spectrum  and  filtered  spectrum  of 

afterglow  produced  by  N  atom  3  body  recombination  in  pulse  rf  cell 
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and  thus  for  N2(B)  states  in  the  post-discharge,  where  ^=10, 11,  or  12  and  v"— v'-4, 
as 

lN.(B,v01 _  _ fv',v"/-^v',v" _  /0 

[N2(B,  V  =  10,  11, 12)]  IiO,6/Aio,6  +  Ill,?/ All, 7  +  Il2,8/ Al2,8 

Table  8  shows  the  resulting  fractional  distribution  of  [N2(B,v)]  using  Eq.  6.16  and 
the  normalized  afterglow  emission  intensities  in  Fig.  21.  The  results  show  that 


(V,v") 

Av',v"(s  ^) 

Ref  [67] 

Ill, 7 

[N2(B,v')J 

1N2(B,v=10,11,12)] 

(12,8) 

1.06x10^ 

0.282 

(11.7) 

9.58x10^ 

1.D0 

0.512 

(10,6) 

8.13x10^ 

0.205 

Table  8:  N2(B)  vibrational  state  population  distribution  20  ms  after  rf  discharge  pulse 


[N2(B,v=10)]  accounts  for  ~20%  of  the  total  N2(B)  state  density  excited  by  atomic 
recombination. 

6.3.6  Microwave  discharge  cell 

The  second  discharge  source  for  this  experiment  was  a  flowing  gas  microwave  system 
with  an  AsTex  S-IOOO  power  source  shown  in  Fig.  18(b)  and  described  in  Section 
3.2.  Microwave  power  was  maintained  at  100  W  transmitted  with  near  zero  reflected 
power  for  all  gas  flow  and  mixture  conditions.  A  total  gas  flow  condition  of  50  seem 
per  Torr  was  maintained  to  keep  the  gas  flow  velocity  and  residence  time  constant 
within  the  discharge.  Because  of  limited  optical  access,  laser  diagnostics  could  not 
be  performed  directly  in  the  discharge  volume  of  this  system,  but  were  done  a  few 
centimeters  downstream  from  the  discharge.  The  atomic  nitrogen  concentration  in 
the  post  discharge  region  was  monitored  by  the  TALIF(207)  scheme  described  in 
Chapter  3.  A  relatively  short  gas  transport  time  between  the  discharge  and  TALIF 
region  was  desired  so  that  changes  in  N  atom  wall  loss  coefficients  as  a  function  of  gas 
mixture  would  have  a  minimum  effect  on  the  resulting  N  atom  concentration.  Since 
the  excited  discharge  products  had  tens  of  milliseconds  of  collisional  decay  time  before 
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they  reached  the  laser  probe  region,  the  N2(A3S+,v=0)  and  N2(B®n3,v=l)  states  were 
depleted  below  the  detectable  limits  of  the  LIF  techniques.  A  strong  amber  colored 
N2(B®nj)— ^N2(A^E+)  emission  was  observable  in  the  laser  probe  region.  The  notable 
intensity  of  this  downstream  afterglow  emission  was  an  indication  of  the  relatively 
high  atomic  production  within  the  microwave  discharge,  since  the  N2(B,v=10,ll,12) 
emitting  states  were  populated  by  three-body  recombination  of  atomic  nitrogen. 

Without  the  luxury  of  the  LIF  probes,  the  behavior  of  N2  excited  states  were  mon¬ 
itored  through  the  emission  intensity  from  the  first  and  second  positive  systems  and 
the  first  negative  system  within  the  disdiarge  cavity  region.  A  single  planar  mirror 
reflected  the  discharge  cavity  emission  which  was  collected  by  a  sin^e  7.5  cm  focal 
length  lens  and  focused  onto  the  slits  of  a  Digikrom  240  monochrometer.  A  Hama¬ 
matsu  R928  PMT  was  mounted  at  the  output  slits  of  the  monochrometer.  For  the  first 
positive  emission,  the  N2(B,v=4— >A,v=2)  transition  was  monitored  at  750  nm,  the 
second  positive  system  was  monitored  at  337  nm,  representing  the  N2(C,v=0->B,v=0) 
transition.  For  the  first  negative  system,  the  (B,v=0-)'X,v=0)  transition  at  391 
nm  was  used.  Emission  from  the  NH  radical  was  also  monitored  within  the  discharge 
at  336  nm  representing  the  NH(A-X)  transition. 

6.4  Results  and  Discussion 

The  temporal  density  measurements  of  the  discharge  products  in  the  pulsed  rf  cell 
with  pure  N2  are  presented  first  with  a  calibration  of  absolute  concentration  made 
for  several  of  the  species.  The  behavior  of  the  measured  species  in  mixtures  of  N2/H2 
is  then  presented  along  with  a  discussion  of  the  evident  changes  in  heavy  particle 
kinetics  as  a  function  of  gas  mixture.  The  effects  of  N2/H2  mixtures  in  the  flowing 
microwave  discharge  are  then  analyzed  with  the  optimal  mixture  for  enhanced  N 
atom  production  identified.  Finally  a  mixture  of  N2/NH3  is  studied  in  the  flowing 
microwave  discharge  and  the  results  compared  to  that  of  the  N2/H2  mixtures. 
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6.4.1  Calibration  of  N  atom,  N2(A),  N2(B),  and  N2(X,v)  ab¬ 
solute  concentrations  in  the  pulsed  rf  parallel  plate  dis¬ 
charge 


Differential  rate  equations  for  the  post-discharge  period  involving  N,  N2(X^E+,v>5), 
N2(A®E+),  N2(B®np),  N2(X^E+,v>12),  and  N2(B®nj„v=10-12)  were  expressed  in  Eqs. 
6.5  through  6.10.  A  numerical  model  using  Euler’s  method  has  been  applied  to  this  set 
of  six  equations.  The  differential  rate  equations  for  eadi  species,  Mj,  where  i=l,...,6, 
are  of  the  form 


dMj 

dt 


/i(t,  Ml,  ...jMe) 


(6.19) 


with  the  formula  for  the  Euler  method  being 


Mi, jfc+i  =  Mi^  -f  hfi{tk,  Mi,fc, ...,  Me,*;)  (6.20) 

which  advances  a  solution  by  a  time  interval  of  h  from  time  t*;  to  t*;+i=t*;-t-h.  A  time 
stepsize  interval  value  of  h=10“^  seconds  was  used  here  which  provided  the  necessary 
stability  and  accuracy  for  the  numerical  analysis.  The  initial  density  values  of  each 
of  the  species  were  adjusted  to  produce  the  best  fit  to  a  series  of  LIE  and  emission 
data  taken  after  a  20  ms  rf  discharge  pulse  at  5  Torr  N2.  The  data  that  was  used 
included  (i)  TALIF  measurements  of  absolute  [N]  from  2  to  20  ms  after  the  discharge 
was  turned  off,  (m)  emission  measurements  of  relative  [N2(B,v=10-12)]  from  0.2  to  30 
ms  after  discharge-off,  (iii)  LIF  measurements  of  the  relative  decay  of  [N2(A)]  from 
0.2  to  1  ms  after  discharge-off,  and  {iv)  LIF  measurements  of  the  relative  decay  of 
[N2(B)]  from  0.2  to  1  ms  after  discharge-off. 

A  unique  set  of  initial  density  values  was  determined  that  fit  the  [N]  temporal 
data  as  well  as  the  temporal  data  of  relative  measurements  of  the  excited  N2  states. 
This  was  accomplished  by  first  establishing  a  relationship  between  the  specific  state 
population  [N2(B®nj,v=10-12)]  and  the  total  population  [N2(B®np)].  Information 
of  relative  N2(B,v)  population  was  obtained  by  analyzing  the  first  positive  system 
spectrum  during  a  time  period  when  an  overlap  in  N2(B®np)  LIF  and  N2(B®np,v=10- 
12)  emission  data  occurred.  An  extended  spectrum  is  shown  in  Fig.  22(a)  of  the 
N2(B®np,v'-^A®S+,v")  emission  taken  0.2  ms  after  a  20  ms  rf  pulse  in  5  Torr  of  N2. 
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Figure  22:  (a)  NzCB^npy-^  ,N2(A3E+,v")  emission  spectrum  0.2  ms  after  rf  discharge 
pulse  of  20  ms  duration  at  5  Torr  N2  (b)  N2(B)  vibrational  population  distribution 
relative  to  [N2(B,v=10)] 
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The  spectrum  includes  emission  from  N2(B)  vibrational  states,  v — 2  through  v— 12 
where  the  spectral  response  of  the  system  was  calibrated  using  a  tungsten  ribbon 
standard  lamp. 

Measurements  of  emission  from  v=0  and  1  were  beyond  the  spectral  range  of  the 
detector  being  used.  Using  these  ejqperimental  emission  intensities  with  Eq.  6.16,  the 
relative  distribution  of  vibrational  states,  v=2  through  12,  normalized  to  the  v=10 
state  is  plotted  in  Fig.  22(b).  The  experimental  points  agree  very  well  with  relative 
N2(B,v)  population  data  generated  by  Plain  et  oJ.[68]  and  Nagpal  et  ol.[69],  which 
both  display  a  sharp  drop  in  population  between  v=2  and  3  and  small  population  rises 
at  v=6  and  v=10.  Considering  the  similarity  of  the  aforementioned  data  with  the 
v=2  through  12  relative  populations  in  Fig.  22(b),  the  final  two  populations  at  v=0 
and  v=l  were  fit  according  to  the  data  of  Nagpal  et  al.[69].  A  vibrational  temperature 
determined  from  v=0  through  2  appears  to  be  much  less  than  one  determined  using 
higher  vibrational  levels.  The  v=0  relative  population  of  Nagpal  rises  only  ~66% 
above  v=2  and  v=l  rises  -^50%  above  that  of  v=2. 

A  ratio  of  [N2(B)]/[N2(B,v=10-12)]=88  was  determined  from  the  N2(B,v)  distri¬ 
bution  in  Fig.  22(b)  at  0.2  ms  after  an  rf  pulse  of  20  ms  duration  in  5  Torr  of  N2. 
This  ratio  couples  the  initial  values  of  [N2(B®ng,v=10, 11,12)]  and  [N2(B*n5)]  for  the 
numerical  model  of  the  post-discharge  rate  equations  for  the  given  conditions. 

The  numerical  model  solution  of  Eqs.  6.5  through  6.10  with  the  best  fit  to  the 
data  resulted  in  the  temporal  plot  of  the  densities  of  the  species  as  shown  in  Fig.  23. 
The  initial  density  values  that  were  used  in  the  fit  are  listed  in  Table  9.  Of  these 
values,  only  [N2(X,v  >  5)],  [N2(A)],  [N2(B)]  and  [N2(X,v  >12)]  were  variables  in  the 
fitting  routine,  while  [N]  had  been  measured  absolutely  and  [N2(B,v=10,ll,12)]  was 
dependent  on  [N2(B)] 


[N] 

[N2(X,v  >  5)] 

[N2(A)] 

INzCB)] 

[N2(X,v  >12)] 

[N2(B,v=10,11,12)] 

7.9x10^=* 

1x10^^ 

8x10^^ 

8x10** 

9x10^=" 

5x10^ 

Table  9:  Initial  density  values  in  cm“®  for  numerical  model  resulting  in  best  fit  of 
data  for  20  ms  rf  pulse  at  5  Torr 
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By  identifying  some  general  trends  in  the  numerical  analysis  of  the  differential 
rate  equations,  approximations  can  be  made  which  simplify  the  kinetic  relationships 
between  the  species  under  certain  post-discharge  conditions.  The  most  useful  trend 
in  the  rate  equations  involves  the  value  of  d[N2{B)]/dt  in  the  post-discharge,  which 
equates  to  the  difference  in  the  production  and  loss  rate  terms  of  [N2(B)]  in  Eq.  6.7. 
The  numerical  model  for  the  conditions  of  a  20  ms  rf  pulse  at  5  Torr  showed  that  the 
loss  rate  for  [N2(B)]  was  only  slightly  greater  than  the  production  rate  in  the  post¬ 
discharge.  Specifically,  during  the  time  period  of  0.1  to  <~2  ms,  when  -d[N2{B)]ldt 
is  seen  in  Fig.  23  to  be  the  steepest,  the  sum  of  the  production  rate  terms  is  still 
calculated  to  be  within  0.1%  of  the  loss  rate  terms.  During  the  time  period  from  3 
to  30  ms,  the  difference  in  the  production  and  loss  rates  drops  to  within  ~0.001%. 
This  trend  justifies  the  approximation  that  the  production  and  loss  terms  of  Eq.  6.7 
be  set  equal  from  the  time  0.1  ms  after  the  rf  pulse  out  to  at  least  30  ms.  This  allows 
Eq.  6.7  to  be  simplified  to  the  steady  state  approximation  of 


— ^(e)  [fV2  (X)]  [iV2  (B)]  -f  kz{d)  [-AT]  [AT]  [JV2]  =  —kz{a)  [^^2  {X,  v  >  5)]  [iV2  (A)] 

-keic)[N2{A)][N2iA)] 


+kei<D[N2{X)][N2{B)] 

+hMB)]  (6,21) 


which  is  a  useful  form  since  all  of  the  terms  on  the  RHS  of  Eq.  6.21  also  reside  in  the 
rate  equation  for  N2(A),  Eq.  6.6.  It  should  be  noted  that,  in  terms  of  kinetics,  the 
first  term  on  the  LHS  of  Eq.  6.21  represents  a  loss  mechanism  of  N2(A)  states  that 
are  excited  by  collision  to  the  N2(B)  state,  but  do  noi  relax  back  to  the  N2(A)  state. 
Substituting  Eq.  6.21  into  Eq.  6.6,  the  rate  equation  for  N2(A3S+)  then  becomes 


d[N2iA)] 

dt 


-*4(.)[iV][Ar2(X)]  -  ^(JV,(^)1 
-ke^4N2{X)][N2{B)]  +  k^4N][N][N2] 


(6.22) 


which  shows  the  major  losses  of  N2(A)  in  the  post-discharge  to  be  from  (i)  colli- 
sional  quenching  of  N2(A)  to  N2(X)  by  atomic  nitrogen,  {ii)  diffusion  loss  of  N2(A) 
to  the  walls,  and  {Hi)  reduction  of  N2(A)  by  excitation  to  N2(B)  followed  by  colli- 
sional  quenching  of  N2(B)  to  N2(X)  by  N2(X).  The  only  production  term  of  N2(A) 
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Time  after  discharge  (ms) 


Figure  23:  Temporal  behavior  of  species  in  post-discharge  where  points  are  e3q)eri 
mental  data  and  solid  lines  represent  the  fit  generated  by  the  numerical  model 


represented  in  this  post-discharge  kinetic  model  is  from  N  atom  recombination,  where 
most  of  the  N2(B,v=10,ll,12)  states  created  in  the  three-body  recombination  process 
will  cascade  down  to  N2(A)  by  collisional  quenching  or  radiation.  Note  the  sharp 
knee  in  decay  of  N2(A)  at  ms.  The  initial  rapid  decay  of  [N2(A)]  occurs  be¬ 
cause  the  quendiing  loss  term,  -k4(o)[N][N2(A)],  is  initially  much  greater  than  the 
three-body  recombination  production  term,  k3(<i)  [N][N][N2],  which  is  independent  of 
[N2(A)].  The  sharp  knee  in  the  plot  occurs  when  [N2(A)]  has  decreased  to  a  value 
where  the  quenching  and  production  processes  equilibrate.  Thereafter,  the  slow  decay 
in  [N2(A)]  depends  on  the  decay  of  [N].  The  same  general  shape  is  seen  in  the  [N2(B)] 
model  due  to  its  coupling  with  N2(A).  Since  the  number  densities  of  the  species  dur¬ 
ing  this  time  period  did  not  allow  for  accurate  LIF  measurements,  the  sharpness  of 
the  transition  in  decay  rates  were  not  confirmed  by  experiment. 

Equation  6.21  can  also  be  rearranged  to  determine  a  relationship  for  [N2(B)]  during 
the  post-discharge  as 

fAofS')!  =  >  5)][Ar2(A)]  -i-  A:6(c)[A2(A)][Ar2(A)]  +  A:3([<i)[iV][Ar][iV2] 

(A:6(rf)  +  A:6(e))[A2(A)]-M/rB 

(6.23) 

In  the  numerical  analysis  with  the  20  ms  rf  pulse  at  5  Torr,  the  production  term  from 
the  associative  reaction  of  N2(A)  was  found  to  be  <1%  of  total  production  at  0.1  ms 
after  the  rf  pulse  was  tumed-off,  which  dropped  to  <0.001%  for  the  time  period  1  ms 
to  30  ms  after  the  rf  pulse.  The  assumption  can  be  made  that  k6(c)[N2(A)][N2(A)]  is 
negligible  compared  to  other  production  terms  under  these  discharge  conditions  and, 
in  general,  can  be  neglected  whenever  [N2(A)]  <  [N2(X,v  >  5)].  It  was  also  found  in 
the  numerical  analysis  for  the  time  period  of  0.1  to  1  ms  after  the  rf  pulse  that  the 
production  term  involving  N  atom  recombination,  k3(d)[N][N][N2],  was  <1%  compared 
to  k6(o)[N2(X,v  >  5)][N2(A)].  The  approximation  that  the  k3(d)[N][N][N2]  production 
term  may  also  be  neglected  from  0.1  ms  to  1  ms  is  valid  as  well.  Thus,  under  the 
given  discharge  conditions  for  the  time  period  0.1  to  1  ms  after  the  rf  pulse,  Eq.  6.23 
can  be  reduced  to  an  expression  for  [N2(B)]  of 


\No(B)]  =  %)[^2(A,u>5)][Ar2(A)] 

^  ^  (A:6(d)+A:6(e))[iV2(A)]-|-l/rB' 


(6.24) 
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The  specialized  conditions  which  make  Eq.  6.24  a  valid  approximation  and  create 
a  simple  relationship  between  N2(X,v>5),  N2(A),  and  N2(B),  have  been  found  to 
exist  under  a  broad  range  of  discharge  conditions.  In  many  esses,  a  window  of  time 
exists  after  a  discharge  pulse  when  this  approximation  can  be  applied.  The  short¬ 
lived  N2(B)  produced  during  the  rf  discharge  has,  in  general,  been  fully  depleted  by 
0.1  ms  after  the  rf  pulse,  but  the  N2(A)  from  the  discharge  is  still  measurable  and  is 
regenerating  measurable  quantities  of  N2(B)  from  collisions  with  N2(X,v  >  5).  This 
represents  the  time  period  that  Eq.  6.24  is  valid  imtil  [N2(A)]  becomes  so  depleted 
that  the  production  of  N2(A)  by  N  atom  recombination  becomes  a  factor,  which 
occurred  around  1  ms  for  the  5  Torr/20  ms  rf  pulse  case. 

An  experimental  test  of  the  validity  of  Eq.  6.24  under  various  discharge  conditions 
was  performed  by  comparing  the  exponential  decay  rates  of  the  N2(A)  and  N2(B) 
LIE  signals  between  0.1  and  1  ms  after  the  pulse.  This  test  was  performed  under  the 
assertion  that  [N2(X,v  >  5)]  experienced  negligible  change  over  this  time  period.  For 
Eq.  6.24  to  be  considered  valid,  the  experimental  data  must  show  that  the  N2(A) 
and  N2(B)  decay  rates  are  similar,  or 

1  rf[iV2(^)] _ 1  d[N2{B)]  (.25) 

[Ar2(A)]  dt  ~  [N2{B)]  dt  '  ^  ^ 

Examples  of  the  LIF  data  iu  Fig.  24  were  taken  at  5  Torr  for  pulse  durations  of  5 
ms,  10  ms  and  20  ms.  This  particular  data  shows  very  similar  decay  behavior  for  the 
N2(B)  and  N2(A)  states.  In  general,  it  has  been  found  that  the  higher  pressure  and 
longer  pulse  data  resulted  in  N2(B)  and  N2(A)  decay  rates  that  were  most  similar. 
The  higher  pressure  and  longer  pulse  discharge  provides  an  environment  beneficial  to 
population  of  the  N2(X,v  >  5)  state  by  vibrational  pumping  and  thereby  promoting 
the  condition  [N2(A)]  <  [N2(X,v  >5)]  which  is  necessary  for  Eq.  6.24  to  be  valid. 
The  data  in  Fig.  24,  therefore,  verifies  that  Eq.  6.24  is  valid  at  5  Torr  for  pulse 
durations  >  1  ms.  The  validity  of  Eq.  6.24  has  also  been  found  to  be  reasonable 
at  3  Torr,  but  very  questionable  under  the  situations  examined  at  1  Torr.  The  rate 
of  collisional  excitation  at  1  Torr  is  proportionally  lower  than  with  the  3  and  5  Torr 
discharge  conditions,  which  results  in  a  lower  [N2(X,v  >5)]. 
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Figure  24:  Comparison  of  N2(A)  and  N2(B)  temporal  decay  rates  at  pulse  durations 
of  (a)  20  ms,  (b)  10  ms,  and  (c)  5  ms  where  the  similarity  in  N2(A)  and  N2(B)  decays 
indicate  that  [N2(A)]  <  [N2(X,v  >  5)]  in  each  case 
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An  expression  to  determine  [N2(X,v  >  5)]  can  be  derived  &om  Eq.  6.24  to  give 


[JV2(A-.o  >  5)] 


[Ar2(B)]  {{hid)  +  A;6(e))[^2(^)]  +  I/tb) 


(6.26) 


^  [N2{A)]  h^a)  '  ' 

This  approximate  expression  is  very  useful  under  the  appropriate  pulsed  disdiarge 
conditions,  since  [N2(B)]/[N2(A)]  can  be  measured  experimentally  during  the  specific 
window  of  time  after  the  rf  discharge  pulse  when  Eq.  6.24  holds  true.  The  time 
of  0.2  ms  after  the  rf  discharge  pulse  was  chosen  to  make  a  variety  of  calibrated 
LIF  and  emission  measurements  under  various  discharge  conditions.  The  results  of 
the  measurements  at  0.2  ms  after  the  rf  pulse  are  given  in  Section  6.4.2  and  include 
calculated  determinations  of  [N2(X,v  >5)]. 

The  above  derivation  can  be  extended  to  use  the  calibrated  emission  representing 
[N2(B,v=10,11,12)]  to  give  an  approximate  ©qpression  for  the  density  of  the  N2(X,v 
>12)  states  which  possess  enough  vibrational  energy  to  excite  N2(A,v=0)  to  the 
N2(B,v=10)  state.  The  expression  is  then 

(JV2(X,  c  >  12)1  =  ^  ^ 

[A2(AjJ  «6(o) 

A  simplified  expression  for  [N2(B,v=10,ll,12)]  can  be  also  be  found  by  assuming, 
as  for  [N2(B)],  that  the  production  and  decay  rate  terms  of  Eq.  6.8  are  approximately 
equal  between  0.1  and  30  ms  after  the  rf  discharge  pulse.  This  expression  can  be 
written  as 

[1V2(B,p  ^  10, 11, 12)1  =  W  Wf  ^  (6.28) 

^  ^  ’  /j  {k6(d)  +  ke(e))[N2{X)]  +  l/TB 

In  the  numerical  analysis  with  the  20  ms  rf  pulse  at  5  Torr,  the  k3(rf)[N][N][N2(X)]  pro¬ 
duction  term  proved  to  be  a  factor  of  20x  greater  than  the  k6(a)[N2(X,v  >12)][N2(A)] 
production  term  for  these  specific  states.  For  the  conditions 


(6.27) 


*3MCT[A1[Ar2(A-)]  >  *6(.)[Ar2{X,t>>12)l[Af2(A)] 
+  *6W)[JV2(X)1  »  1/tb 
a  simplification  can  be  made  to  Eq.  6.28  giving 

[iV2(B.v  =  10,ll,12)l  =  ^^^fg^. 


(6.29) 


(6.30) 
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This  expression  is  the  approximation  behind  the  common  assertion[23]  that  the  N 
atom  density  is  directly  proportional  to  the  square  root  of  the  nitrogen  afterglow 
emission  from  the  first  positive  system. 

6.4.2  Pulsed  rf  parallel  plate  discharge  results 

The  absolute  atomic  nitrogen  concentration  measured  at  the  center  of  the  rf  cell  is 
shown  in  Fig.  25  for  various  discharge  pulse  durations  up  to  20  ms.  The  N  atom 
measurement  for  each  pulse  duration  was  done  early  in  the  post-discharge,  0.2  ms 
after  the  discharge  was  turned  off.  This  plot  represents  the  rate  of  build  up  of  the 
atomic  species  as  a  function  of  the  “discharge  on”  time.  The  N  atom  density  at  5 
Torr  is  seen  to  continue  to  rise  beyond  the  20  ms  time  period.  This  is  a  result  of 
the  diffusion  rate  of  the  atomic  nitrogen  through  N2  slowing  at  higher  pressure  and 
resulting  in  a  longer  time  period  for  the  atomic  production  to  reach  equilibrium  with 
the  surface  losses. 

The  trend  of  increasing  [N]  with  longer  rf  disdiarge  pulse  length  explains  the 
measured  increasing  rate  of  N2(A)  and  N2(B)  decay  rates  for  longer  pulses  shown  in 
previously  in  Fig.  24.  The  greater  N  atom  density  accelerates  the  depletion  of  N2(A) 
by  the  primary  coUisional  quenching  reaction  of  N  -f  N2(A). 

A  plot  is  shown  in  Fig.  26  of  the  evolution  of  each  of  the  species  in  the  center  of 
the  rf  cell  at  5  Torr  of  N2  as  a  function  of  rf  pulse  duration.  The  temporal  behavior 
of  the  N2(A)  state  is  quite  interesting  as  it  builds  to  a  substantial  concentration  for 
pulse  durations  less  than  1  ms,  but  is  depleted  for  longer  pulses.  The  atomic  nitrogen 
density  has  had  time  to  build  during  the  longer  pulses  and  is  responsible  for  increasing 
the  quenching  rate  of  N2(A)  and  lowering  its  concentration.  The  smoothness  of  the 
[N2(X,v  >  5)]  evolution,  which  was  calculated  from  the  vastly  different  temporal 
functions  of  [N2(A)]  and  [N2(B)],  is  an  indication  that  the  kinetic  derivations  produced 
an  accurate  result.  The  N2(X,v  >  5)  density  approaches  equilibrium  at  a  slightly 
slower  rate  than  the  measured  N  atom  density,  most  likely  due  to  the  smaller  diffusion 
rate  of  N2(X,v)  through  N2.  Neither  species  is  readily  quenched  by  volume  reactions 
and  both  have  a  relatively  low  wall  loss  rate. 
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Figure  25:  Plot  of  atomic  nitrogen  density  measured  by  TALIF  in  the  center  of  the 
pulsed  rf  pure  N2  discharge  cell  for  rf  pulse  durations  of  0.05  ms  to  20  ms 
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Figure  26:  Temporal  evolution  of  [N2(X,v  >  5)],  [N],  [N2(A)],  and  [N2(B)]  at  5  Torr 
N2  in  rf  cell  with  varying  pulse  durations 
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The  plot  in  Fig.  27  shows  the  measured  [N2(B,v=10)]  along  with  the  calculated 
evolution  of  [N2(X,v  >  12)]  under  the  same  conditions  as  Fig.  26.  On  close  inspection, 
it  is  seen  that  the  higher  v-states  of  N2(X)  do  not  increase  as  quickly  with  increasing 
pulse  length  as  did  [N2(X,v  >  5)].  A  result  is  that  the  [N2(B,v=10)]  actually  decreases 
since  the  N2(X,v  >  12)  population  gain  can  not  counter  the  depletion  rate  of  its  co¬ 
reactant,  N2(A),  for  pulse  durations  greater  than  1  ms. 

The  plot  in  Fig.  28  again  investigates  the  behavior  of  [N2(X,v  >  5)],  similar  to 
Fig.  26,  except  at  a  pressure  is  3  Torr.  A  decrease  in  collisional  effects  is  noticeable 
in  this  plot.  The  [N2(A)]  evolution  does  not  turn  over  quite  as  quickly  or  deplete  as 
much  as  at  5  Torr  due  to  a  lower  N  atom  concentration.  The  atomic  nitrogen  and 
N2(X,v  >  5)  densities  also  tend  toward  an  equilibrium  concentration  faster  due  to 
increased  diffusion  to  the  walls. 

With  the  densities  and  behavior  of  the  measured  species  well  characterized  at 
3  and  5  Torr  of  N2,  hydrogen  was  then  added  in  small  amoimts  into  the  pulsed 
rf  discharge  cell.  Since  the  quenching  of  N2  electronic  metastables  by  hydrogen  has 
been  an  important  aspect  of  recent  theories  involving  N2/H2  discharges[20,  21],  direct 
observation  of  this  effect  in  the  pulsed  rf  cell  was  a  goal  of  this  experiment.  Figure 
29  shows  the  effect  of  the  first  4%  of  H2  introduced  into  the  N2/H2  discharge  on 
the  population  of  the  N2(A^S^)  metastable  and  the  N2(B®nj)  excited  state,  both 
measured  by  LIF.  As  can  be  seen,  the  densities  of  the  coupled  N2(A)  and  N2(B)  states 
are  reduced  by  more  than  an  order  of  magnitude  by  the  slight  addition  of  H2.  Since  the 
[N2(A)]  <C  [N2(X,v  >  5)]  condition  holds  in  the  gas  mixture,  the  N2(X,v  >  5)  density 
can  again  be  calculated  with  Eq.  6.26  using  the  LIF  data.  The  resulting  [N2(X,v 
>  5)]  trend  is  interesting  in  that  it  increases  slightly  as  the  electronic  metastable 
is  quenched  by  hydrogen.  This  result  is  consistent  with  theoretical  work[20]  which 
predicts  an  enhancement  of  the  middle  vibrational  states  of  N2(X)  in  a  2  Torr,  80  mA 
DC  disdiaxge.  This  enhancement  is  said  to  result  from  a  two-fold  effect  of  increased 
electron  impact  excitation  of  vibrational  states  and  a  cascading  of  higher  v-states, 
removed  through  vibrational-translational  energy  exchanges  with  atomic  hydrogen 
and  molecular  hydrogen. 

The  effect  of  N2/H2  mixtures  in  the  rf  pulsed  disdiaxge  has  also  been  studied 
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Figure  27:  Temporal  evolution  of  [N2(X,v  >  12)],  [N],  [N2(A)],  and  [N2(B,v=10)]  in 
5  Torr  N2  rf  cell  with  varying  pulse  durations 
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Figure  28:  Temporal  evolution  of  [N2(X,v  >  5)],  [N],  [N2(A)],  and  [N2(B)]  at  3  Torr 
N2  in  rf  cell  with  varying  pulse  durations 
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Figure  29:  [N2(A)]  and  [N2(B)]  determined  by  LIF  and  [N2(X,v  >  5)]  calculated  from 
[N2(B)]/[N2(A)]  for  H2  added  to  N2  in  a  pulsed  rf  discharge  at  3  Torr 


by  observing  both  atomic  nitrogen  and  N2(A)  over  the  cycle  of  a  20  ms  rf  pulse  at 
3  Torr,  shown  in  Fig.  30.  Measurements  were  taken  for  both  a  pure  N2  disdiarge 
and  then  a  1%  mixture  of  H2  in  N2.  The  behavior  of  the  atomic  nitrogen  was  very 
siTnilar  under  both  gas  conditions.  The  atomic  nitrogen  density  was  elevated  slightly 
in  the  hydrogen  mixture,  but  both  cases  were  within  experimental  uncertainty  of  one 
another.  The  hydrogen  seems  to  have  no  effect  on  the  post-discharge  loss  rate  of  the 
N  atoms  either.  The  N2(A)  state,  on  the  other  hand,  ejqperienced  a  much  greater  loss 
rate  with  the  hydrogen  mixture.  During  the  time  period  between  1  ms  to  5  ms  into 
the  discharge  pulse  and  then  immediately  after  the  discharge  turn-off,  it  was  clear 
that  the  N2(A)  was  quenched  at  a  much  higher  rate  in  the  N2/H2  mixture.  This  loss 
rate  can  be  accounted  for  by  the  atomic  hydrogen  quenching  process  in  Table  4(d) 
which  is  higher  than  the  associated  quenching  rate  due  to  N  atoms  in  Table  4(a). 
A  hydrogen  fractional  dissociation  within  the  pulsed  rf  disdiarge  of  ~1  to  5%  would 
easily  account  for  the  observed  loss  rate  of  N2(A). 

The  rf  pulsed  discharge  cell  has  allowed  for  some  unique  determinations  of  absolute 
concentrations  of  nitrogen  disdiarge  species  and  their  behavior  within  N2/H2  gas 
mixtures.  The  information  from  the  fundamental  observations  made  within  the  pulsed 
rf  cell  will  be  extrapolated  to  the  kinetic  model  for  the  higher  power  flowing  microwave 
system. 


6.4.3  Flowing  microwave  discharge  results 

The  behavior  of  the  atomic  nitrogen  density  in  the  flowdng  N2/H2  microwave  discharge 
as  a  function  of  %H2  has  been  investigated  by  TALIF  measurements  for  several  pres¬ 
sures  typical  of  downstream  microwave  plasma  processing,  where 

Figure  31  shows  an  increase  in  N  atom  density  in  response  to  very  small  additions  of 
H2  to  the  N2  feed  gas  for  all  pressures  studied.  In  each  case  the  N  atom  concentration 
peaked  and  then  began  to  decline  as  more  H2  was  added.  The  peak  concentration  at 
0.5  Torr  was  a  ~50%  increase  over  pure  N2  while  at  3  Torr,  the  peak  was  limited  to  a 
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Figure  30:  Behavior  of  [N]  and  [N2(A)]  during  and  after  a  3  Torr,  20  ms  rf  pulse 
discharge  with  N2  and  with  1%  H2  in  N2 
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~15%  increase.  For  the  higher  pressure  data  in  Fig.  31,  the  decreasing  trend  in  the 
N  atom  density  was  more  prominent  and  occurred  at  a  lower  %H2. 

A  kinetic  model  of  atomic  nitrogen  production  in  this  N2/H2  microwave  discharge 
system  is  established  here  using  the  rate  equations  of  Tables  3  through  6  and  the 
experimental  data  from  optical  diagnostics.  The  primary  channels  of  N2  dissociation 
and  N  atom  destruction  in  the  N2/H2  microwave  discharge  are  listed  in  Table  3.  The 
electron  impact  processes,  3(a)  and  3(b),  have  rates  that  are  highly  dependent  on  the 
electron  energy  distribution  function  and  may  vary  extensively  with  discharge  condi¬ 
tion.  Under  certain  discharge  conditions,  electron  impact  dissociation  of  molecules  in 
the  N2(X^S+)  upper  vibrational  manifold,  or  e-V  dissociation  listed  as  process  3(b), 
can  be  a  significant  contributor  to  the  overall  N  atom  production  rate[46].  Loureiro 
calculated  under  the  conditions  of  5  Torr  of  N2  at  Tp=500  K  and  an  E/N  of  6x10"^® 
V  cm^  that  e-V  dissociation  proceeds  at  a  rate  more  than  an  order  of  magnitude 
greater  than  process  3(a),  which  involves  direct  electron  impact  from  N2(X^E+,v=0). 
The  vibration-vibration  (V-V)  and  vibration-translation  (V-T)  processes,  3(cl)  and 
3(c2)  respectively,  become  important  when  V-V  energy  exchange  propagates  within 
the  N2(X^E+)  manifold  and  produces  a  non-equilibrium  vibrational  distribution  that 
includes  a  significant  population  in  the  last  bound,  v=45  state.  The  example  rate 
coefficient  for  the  combined  V-V  and  V-T  dissociation  processes  under  the  above 
conditions  was  extracted  from  Loureiro  et  cl.  [46]. 

The  model  of  N  atom  production  presented  here  for  the  N2/H2  system  involves 
two  separate  mechanisms  that  are  independently  responsible  for  (i)  the  enhancement 
of  N  atom  density  at  small  %H2  and  (zi)  the  decline  in  N  atom  density  with  further 
increase  of  %H2. 

6.4.4  Kinetics  of  observed  enhancement  in  N  atom  density 
for  small  %H.2 

A  kinetic  mechanism  has  been  established  for  the  increase  in  N  atom  density  that  was 
observed  for  all  pressures  tested  when  small  fractions  of  H2  were  added  into  the  Nj 
microwave  discharge.  First,  it  is  assumed  that  the  small  quantities  of  H2  introduced 
in  the  microwave  discharge  were  dissociated  at  very  near  100%  efficiency,  as  has  been 
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measured  in  previous  experiments [21, 17].  This  impressive  efficienqr  is  mainly  due  the 
collisional  dissociative  process  of  H2  with  the  N2(a'  metastable[20,  21]  as  shown 
in  Table  4(b).  A  similar  process  involving  N2(A  in  Table  4(c)  is  seen  to  have 
a  lower  rate.  This  dissociation  produces  considerable  amounts  of  atomic  hydrogen 
which  proceed  to  quench  the  N2(a'  ^E")  and  N2(A  3S+)  metastables  at  a  very  high 
rate,  as  given  in  processes  4(d)  and  4(e).  This  volume  loss  rate  of  the  nitrogen 
metastables  by  the  abundant  H  atoms  exceeds  the  metastable  volume  loss  rate  in 
pure  N2  where  the  net  volume  loss  of  metastables  was  primarily  due  to  collisions  with 
N  atoms. 

The  depletion  of  the  N2(a'  ^E~)  and  N2(A  ^E+)  metastables  is  of  consequence  to 
various  collisional  processes  which  are  important  in  a  pure  N2  discharge  in  defining  the 
chemical  composition.  Most  notable  is  the  depleting  effects  on  the  associative  ioniza¬ 
tion  f;ha.Tmfils  involving  metastable  collisions.  The  rates  of  the  associative  ionization 
processes,  shown  in  Table  5(a)&(b),  indicate  that  these  processes  are  the  primary 
production  mechanisms  of  Nj  in  a  typical  pure  N2  discharge,  in  excess  of  the  electron 
impact  processes  5(c)&(d).  The  decrease  of  N2(a'  ^E“)  and  N2(A  ^EJ)  metastable 
density  as  H2  is  added  leaves  both  reactants  of  the  associative  ionization  processes 
depleted,  which  should  result  in  a  sharp  decline  in  the  Nj  population  within  a  small 
%H2.  Since  the  associative  ionization  process  relies  on  heavy  body  collisions,  the  de¬ 
cline  in  Nj  population  will  be  most  noticeable  at  higher  pressures.  Evidence  of  this 
decline  in  Nj  can  be  found  by  observing  the  Nj (B^S+  -^X^S+)  emission  from  the 
discharge.  Figure  32  shows  the  decline  of  Nj  emission  intensity  as  H2  is  added  for 
the  pressures  of  interest.  The  magnitude  of  the  drop  in  Nj  emission  intensity  over 
the  first  few  percent  of  H2  added  at  3  Torr  is  seen  to  be  quite  similar  to  the  measured 
decline  in  Fig.  29  of  the  N2(A  ^S+)  metastable  in  a  3  Torr  pulsed  rf  discharge  as  a 
few  percent  of  H2  was  added. 

The  loss  of  this  N2  ionization  flux  has  been  found  to  result  in  the  increase  of 
the  discharge  impedance  in  a  DC  discharge[20],  and  a  similar  increase  is  assumed 
to  occur  in  the  microwave  discharge.  The  applied  microwave  field,  which  is  held  at 
a  constant  forward  power,  responds  to  the  impedance  increase  with  an  increase  in 
the  effective  E/N  in  the  discharge.  This  change  “heats  up”  the  EEDF  within  the 


103 


Emission  Intensity  (arb  units) 


disdiarge,  enhancing  the  population  of  the  high  energy  tail  of  the  EEDF.  The  high 
energy  tail  of  the  EEDF  includes  the  range  of  electron  energies  that  are  responsible 
for  the  electron  impact  dissociation  of  N2  as  in  Table  3(a)&:(b).  Since  the  heating 
of  the  EEDF  tail  accompanies  the  depletion  of  associative  ionization,  the  EEDF 
enhancement  will  follow  the  decline  in  Nj  population.  The  Nj  emission  is  seen  in 
Fig.  32  to  have  a  greater  decline  at  higher  pressure,  which  indicates  that  the  heating 
of  the  EEDF  will  occur  in  the  low  %H2  range  and  will  be  greatest  at  higher  pressures. 

The  observation  of  the  N2(C®ntt  emission  from  within  the  microwave 

discharge  region  can  provide  information  on  the  changes  in  the  EEDF.  The  excitation 
of  the  N2(C*n„)  state,  whidi  lies  11  eV  above  the  ground  state,  can  be  achieved  by  a 
pooling  reaction  between  two  N2(A  states[60],  but  is  also  quite  readily  populated 
by  the  direct  electron  impact  process 

e  +  iV2(X^S+)  ^  iV2(C®n«)  +  e.  (6.32) 

An  actinometry  method  was  used  by  Loureiro  and  co-workers  [24]  within  a  DC  dis¬ 
charge  which  determined  that  excitation  to  the  N2(C^na)  state  in  the  active  disdiarge 
occurred  mainly  by  the  direct  electron  impact  process  in  Eq.  6.32.  The  short  radia¬ 
tive  lifetime  of  N2(C®n„),  coupled  with  its  proximity  in  energy  to  the  N2  dissociation 
threshold  of  9.75  eV,  makes  the  N2(C^IItt->B®n3)  emission  a  good  indicator  of  the 
rdative  mean  electron  energy  and  the  availability  of  electrons  with  sufficient  energy 
for  dissociation  by  direct  electron  impact.  Figure  33  shows  the  N2(C®n« 
emission  at  337  nm  as  H2  is  introduced  under  the  discharge  conditions  of  interest. 
The  data  in  Fig.  33  is  consistent  with  the  proposed  kinetics  as  the  increase  in  emission 
is  within  the  first  few  %H2  and  is  sharpest  at  the  higher  pressures,  which  corresponds 
to  the  conditions  where  the  greatest  depletion  of  associative  ionization  occurs.  The 
slight  decline  in  N2(C*n«->B^np)  emission  at  higher  %H2  is  attributed  simply  to  the 
decreasing  partial  pressure  of  N2  in  the  gas  mixture. 

The  kinetic  model  predicts  that  this  enhancement  in  the  tail  of  the  EEDF  for  small 
%H2  will  lead  to  an  increase  in  atomic  nitrogen  production  by  the  direct  electron 
impact  processes  of  Table  3(a)&(b).  The  data  in  Fig.  31  confirms  this  increased  N 
atom  production  at  all  pressures  for  the  initial  range  of  %H2.  For  the  higher  pressure 
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data,  though,  a  lesser  rise  in  N  atom  production  indicates  that  a  parallel  medianism 
is  causing  a  decline  in  the  atomic  nitrogen  density  as  the  %H2  is  increased  further. 
The  kinetics  of  this  parallel  mechanism  causing  a  decline  in  N  atom  density  will  be 
discussed  in  the  next  section. 

6.4.5  Kinetics  of  observed  decline  in  N  atom  density  at  higher 
%H2 

The  atomic  nitrogen  concentration  reaches  a  peak  at  a  small  %H2  for  each  disdiarge 
pressure,  as  seen  in  Fig.  31,  and  then  decreases  as  further  H2  is  added.  For  N2/H2 
pressures  of  1  Torr  and  above,  the  decline  in  N  atom  density  is  too  rapid  to  be 
eq)lained  by  decreasing  N2  partial  pressure  in  the  mixture.  Neither  can  the  decreasing 
N  atom  concentration  be  attributed  to  a  decline  in  dissociation  by  direct  electron 
impact,  since  the  N2(C®n„->B^np)  emission  in  Fig.  33  does  not  indicate  a  “cooling” 
of  the  high  energy  tail  of  the  EEDF.  Another  theory  that  does  not  adequately  explain 
the  experimental  results  was  presented  by  Gordiets  et  al,  who  proposed  that  a  similar 
decline  in  [N]  observed  in  their  DC  discharge  data  was  due  to  the  reaction 

N(2D)  +  NH->N2  +  H.  (6.33) 

The  result  of  this  reaction  within  a  N2/H2  discharge  is  that  a  metastable  state  of 
atomic  nitrogen  is  removed,  whereas  in  a  pure  N2  discharge  the  metastable  would 
have  returned  to  the  groirnd  N(^S)  state.  A  calculation  by  Gordiets[21]  resulted  in 
a  theoretical  decline  that  was  not  steep  enough  to  match  their  DC  discharge  data 
and  does  not  produce  an  accurate  model  of  the  rapid  decline  in  this  work  either.  A 
further  analysis  of  the  kinetics  of  the  discharge  species  is  needed  to  understand  the 
behavior  of  atomic  nitrogen  at  the  higher  %H2. 

Along  with  direct  electron  impact  dissociation,  N2  may  be  dissociated  as  well 
by  heavy  particle  interactions  involving  vibrationally  excited  N2(X^E+)  states.  To 
give  an  idea  of  the  general  trend  of  the  N2(X)  vibrational  distributions  for  N2/H2 
mixtures,  theoretical  derivations  by  Garscadden  et  al.  for  a  2  Torr  DC  discharge  at 
80  mA  are  displayed  in  Fig.  34  for  0  to  10%  H2.  Dissociation  involving  the  higher 
vibrational  states  of  N2(X)  may  occur  by  V-V  or  V-T  energy  transfer  processes,  as 
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given  in  Table  3(cl)&(c2)  respectively.  These  dissociation  channels  are  depleted  as 
H2  is  added  due  to  the  increasing  rate  of  depopulation  of  high  v  N2(X^E+)  states  by 
V-T  energy  exchange  between  N2  and  H2.  The  depletion  of  N2(X,v  >  40)  states  as 
the  %H2  rises  toward  10%  is  quite  obvious  from  the  distributions  in  Fig.  34.  The 
V-T  rate  coefficient  for  N2/H2  interaction  in  1%  H2  has  been  calculated  to  be  more 
than  an  order  of  magnitude  greater  than  the  N2-N2  V-T  rate  coefficient[24].  This 
theoretical  calculation  suggests  that  the  rate  of  dissociation  by  the  V-V  and  V-T 
processes  drops  by  3  orders  of  magnitude  as  the  %H2  increases  from  0%  to  0.5%  in 
the  N2/H2  mixture.  The  actual  depletion  of  the  high  V  states  of  N2(X^S^)  by  H2  in 
this  experiment  may  be  somewhat  less  than  in  theory  due  to  the  loss  of  H2  by  the 
high  fractional  dissociation  at  low  %H2.  In  any  case,  if  the  V-V  and  V-T  processes  are 
contributing  substantially  to  the  dissociation  of  N2,  the  atomic  nitrogen  production 
will  drop  accordingly  as  hydrogen  is  introduced. 

The  loss  of  electron  impact  dissociation  from  the  mid  to  upper  vibrational  levels 
of  N2(X^S+),  or  e-V  dissociation,  could  be  a  contributor  to  the  observed  [N]  decline 
at  higher  %H2.  For  the  initial  H2  addition  to  N2,  Fig.  34  shows  a  population  rise  of 
[N2(X,5<v<35)]  within  the  vibrational  plateau  region.  This  effect  was  observed  in 
the  [N2(X,v>5)]  behavior  in  the  pulsed  rf  cell  in  Fig.  29  as  well.  The  upper  v^  states 
(v^>25)  within  this  plateau  region  were  then  depleted  as  the  %H2  was  increased 
to  several  percent.  Dissociation  may  occur  from  these  high  v'  states  by  impact  of 
electrons  with  energies  less  than  4  eV,  of  which  there  are  several  orders  of  magnitude 
more  than  the  10  eV  electrons  needed  for  dissociation  from  v=0[20].  The  process  of 
e-V  dissociation,  therefore,  can  be  modeled  to  experience  a  more  gradual  decline  over 
a  %H2  of  1-10%  as  the  density  of  v^>25  states  of  N2(X)  is  gradually  depleted. 

Dissociation  by  V-V,  V-T,  or  e-V  energy  transfer  is  typically  more  prominent  for 
the  higher  pressure  regimes  where  a  high  frequency  of  heavy  particle  collisions  allow 
the  transfer  of  vibrational  energy  throughout  the  N2(X,v)  manifold.  Therefore,  the 
depletion  of  these  N2  dissociation  channels  will  have  a  more  obvious  effect  on  the  N 
atom  density  at  higher  pressures.  The  data  in  Fig.  31  clearly  support  this  aspect  of 
the  kinetic  model  as  the  atomic  nitrogen  concentration  at  3  Torr  declines  by  almost 
one-third  of  its  peak  value  as  further  H2  is  added,  while  at  0.5  Torr  there  is  essentially 
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no  decline  beyond  that  due  to  the  decreasing  N2  partial  pressure. 

6.4.6  N2/NH3  gas  mixtures  in  the  flowing  microwave  dis¬ 
charge 

As  an  alternative  to  the  addition  of  molecular  hydrogen  to  the  N2  microwave  dis¬ 
charge,  NH3  was  added  in  small  percentages.  The  effects  of  the  ammonia  on  the 
downstream  density  of  atomic  nitrogen  was  monitored  by  TALIF  and  the  results  are 
shown  in  Fig.  35.  The  results  of  the  N2/NH3  mixtures  are  strikingly  similar  to  the 
N2/H2  data.  For  each  pressure,  the  N  atom  density  displays  a  modest  increase  for 
%NH3  <  1%  and  then  declines  as  the  %NH3  is  increased  up  to  10%.  The  heavy  par¬ 
ticle  kinetics  described  for  the  N2/H2  mixtures  may  also  be  applied  to  the  N2/NH3 
discharge,  since  in  both  cases  the  additive  precursor  will  be  dissociated  to  a  high  de¬ 
gree,  leaving  similar  discharge  products  containing  N  and  H.  In  comparing  the  mixes, 
it  appears  that  each  trend  in  the  N  atom  density  occurs  with  a  slightly  higher  per¬ 
centage  of  NH3  than  with  H2.  This  could  be  an  indication  that  the  heavy  particle 
dissociative  collisions  that  occur  with  the  initial  introduction  of  H2,  may  occur  at  a 
slightly  lower  rate  for  NH3  than  for  H2. 

6.5  Conclusion 

Two  separate  high  frequency  N2/H2  discharge  systems  were  studied  by  optical  diag¬ 
nostic  techniques  with  a  quantitative  analysis  performed  on  several  of  the  discharge 
products.  Experimental  measurements  in  a  low  power  rf  parallel  plate  pulsed  dis¬ 
charge  and  a  flowing  microwave  discharge  each  provided  unique  information  on  the 
kinetics  within  a  N2/H2  discharge.  The  pulsed  rf  discharge  cell  allowed  for  a  combi¬ 
nation  of  optical  diagnostics  for  the  detection  of  N('‘S),  N2(A®E+),  N2(B^n,)  states. 
The  diagnostic  data  were  fitted  by  a  numerical  kinetic  model  which  calibrated  the  rel¬ 
ative  diagnostic  signals  to  the  absolute  density  of  these  states  and  provided  a  method 
to  derive  [N2(X^E+,v)]  from  the  measured  quantities. 
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The  evolution  of  the  absolute  concentrations  of  these  discharge  products  was 
studied  for  N2  and  N2/H2  discharges.  At  3  and  5  Torr  N2,  the  densities  [N]  and 
[N2(X^E^,v>5)]  measured  0.2  ms  after  the  rf  pulse  were  found  to  increase  with  in¬ 
creasing  rf  pulse  duration.  The  rate  of  [N]  and  [N2(XiS+,v>5)]  build  up  tapered  off 
at  the  longer  pulses  due  to  diffusion  losses.  The  measured  [N2(A)]  at  5  and  3  Ton- 
increased  for  the  pulse  durations  from  0  to  '^.1  ms  but  then  decreased  as  the  pulse 
duration  lengthened.  This  decrease  was  found  to  be  caused  by  an  increasing  quench¬ 
ing  rate  of  N2(A)  by  N  atoms.  The  N2(B)  density  was  found  to  build-up  within  the 
first  ms  pulse  length  and  then  remain  roughly  constant  for  longer  pulse  durations, 
despite  the  changing  densities  of  the  other  species. 

The  N  atom  density,  measured  in  the  microwave  system,  showed  a  modest  in¬ 
crease  with  the  initial  introduction  of  H2  and  then  declined  as  further  H2  was  added. 
The  increase  in  N  atom  density  with  the  addition  of  a  small  %H2  was  suggested  to 
occur  due  to  an  enhancement  of  the  high  energy  tail  of  the  EEDF  which  increased 
N2  dissociation  by  direct  electron  impact.  Experimental  observations  showed  a  rise 
in  N2(C®n„— discharge  emission  and  a  decline  in  Nj(B^S+  X^S+)  emis¬ 
sion  as  the  %H2  was  increased.  These  observations  are  consistent  with  the  theory 
of  an  enhancement  in  the  high  energy  tail  of  the  EEDF  and  a  large  reduction  in 
the  associative  ionization  involving  N2(a'  and  N2(A^S+)  metastables.  Direct 
LIF  measurements  of  [N2(A3e+)]  in  the  rf  cell  with  the  addition  of  H2  showed  an 
rapid  quenching  of  this  metastable  by  H  atoms.  This  quenching  is  also  expected  for 
N2(a  ^E^).  An  increase  in  mid  level  N2(X)  vibrational  population  was  also  deter¬ 
mined  experimentally  for  small  %H2,  which  is  consistent  with  theoretical  predictions. 

The  eventual  decline  in  N  atom  density  while  %H2  was  increased  from  1  to  10% 
was  attributed  to  the  changing  N2(X^E+)  vibrational  distribution  and  the  resulting 
decline  in  the  dissociative  rates  from  processes  involving  N2(XiE+,v).  Experimental 
evidence  of  the  behavior  of  the  upper  v>25  states  of  N2(X)  could  not  be  obtained  due 
to  the  difficulty  in  state  selective  detection  of  these  states.  Theoretical  modeling  in 
previous  studies  of  the  N2(X^E+ )  vibrational  distributions  for  various  N2/H2  mixtures 
supports  the  findings  of  an  increase  in  [(N2(XiE+,v>5)]  for  small  %H2  in  the  rf  cell  as 
well  as  a  more  gradual  depletion  of  [N2(X^E+,v>35)]  as  the  %H2  rises  to  10%.  Peak 


112 


N  atom  densities  in  a  100  W  microwave  discharge  were  found  to  occur  at  3  Torr  with 
less  than  1%  H2  in  the  N2/H2  mixture.  The  same  percent  mixture  of  NH3  in  N2  was 
found  to  affect  the  atomic  nitrogen  density  in  a  very  similar  way. 
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CHAPTER  7 

N  atom  destruction  channels  in  rf  peirallel  plate 

discharge 


7.1  Introduction 

In  a  pure  N2  discharge,  the  net  volume  loss  of  atomic  nitrogen  is  limited  to  a  relatively 
slow  three-body  recombination  process  [22].  Any  additional  recombination  depends 
on  the  reactor  surface  geometry  and  surface  materials  within.  The  optimum  state 
of  operation  for  many  plasma  processing  reactors  requires  the  atomic  nitrogen  flux 
to  a  substrate  surface  to  be  maximized  [70].  To  achieve  the  maximum  atomic  flux, 
care  must  be  taken  to  maintain  the  low  volume  loss  of  atomic  nitrogen,  as  well  as  to 
minimize  the  atomic  losses  due  to  the  interaction  of  N  atoms  with  surfaces. 

In  this  chapter,  the  atomic  nitrogen  volume  and  surface  losses  within  a  parallel 
plate  rf  Na  discharge  are  quantifled  using  the  TALIF(207)  technique  introduced  in 
Chapter  3.  The  recombination  rate,  7,  of  atomic  nitrogen  was  determined  on  stainless 
steel,  aluminum,  silicon,  and  boron  nitride  surfaces  within  the  active  discharge  of 
the  pulsed  rf  cell,  where  7  is  defined  as  the  probability  that  an  atom  impinging  on 
the  surface  will  recombine.  The  metal  surfaces,  aluminum  and  stainless  steel,  were 
investigated  since  these  are  commonly  used  in  discharge  reactor  construction.  A  pair 
of  electrodes  was  made  from  each  of  the  metals  and  incorporated  into  the  reactor.  The 
atomic  loss  was  measured  near  the  electrode  surfaces  facing  the  discharge.  A  silicon 
wafer,  which  represented  a  common  processing  surface,  was  loaded  onto  the  bottom 
aluminum  electrode  and  the  N  atom  loss  rate  was  measured  nearby  in  a  manner 
similar  to  the  metals.  Boron  nitride  wafers  were  investigated  since  this  is  a  material 
that  manufacturers  boast  to  be  chemically  inert.  Boron  nitride  could  potentially  be 
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integrated  into  the  reactor  as  a  replacement  material  or  even  a  coating  on  a  more 
reactive  surface.  For  each  surface,  the  loss  rate  was  determined  from  the  temporal 
decay  of  the  N  atoms  measured  by  TALIF  near  the  surface  being  studied.  The  surface 
materials  within  the  reactor  were  all  subjected  to  a  preparatory  discharge  cleaning  to 
eliminate  foreign  impurities.  No  attempt  was  made,  though,  to  remove  any  surface 
effects  caused  by  prolonged  ex3)osure  to  the  active  nitrogen.  This  means  that  the 
measurements  may  be  applied  as  a  suitable  gauge  of  N  atom  recombination  on  the 
surface  material  as  it  exists  in  an  N2  processing  discharge,  but  it  is  not  necessarily  a 
measure  of  N  atom  recombination  on  the  material  in  its  pmre  state. 

The  N  atom  TALIF  data  from  within  the  rf  cell  was  fit  by  the  output  of  a  two- 
dimensional  numerical  model  which  simulates  atomic  diffusion,  volume  loss  and  sur¬ 
face  loss  within  the  rf  reactor.  By  using  a  least  squares  fit  of  the  model  to  the  N  atom 
decay  data  near  a  surface,  the  recombination  rate  of  that  surface  was  determined. 
Since  this  numerical  modeling  approach  to  solving  for  the  recombination  rate  was 
novel,  the  rationale  and  procedure  is  discussed  in  detail  in  Section  7.3  and  Appendix 
A. 

When  employing  a  low  pressure  N2  discharge  as  an  atomic  nitrogen  source,  the 
use  of  pure  N2  feed  gas  should  limit  the  volume  loss  of  N  atoms  to  the  slow  three-body 
recombination  process.  Any  impurity  in  the  system  may  alter  the  volume  loss  and 
affect  the  resulting  atomic  flux.  The  infiltration  of  air  or  water  vapor  is  the  most 
common  contamination  of  low  pressure  discharge  cells.  The  effect  of  sudi  contam¬ 
inants  was  investigated  by  measuring  the  destruction  of  atomic  nitrogen  due  to  a 
controlled  introduction  of  O2  impurity  into  a  parallel  plate  pulsed  rf  N2  discharge. 
Spatial  and  temporal  TALIF  measurements  found  a  significant  N  atom  loss  from  very 
small  additions  of  O2  in  a  reactor  that  had  otherwise  very  low  surface  and  volume 
loss  rates. 

7.2  Experiment 

The  rf  parallel  plate  discharge  cell  used  in  this  study  is  described  in  detail  in  Section 
3.3.  Electrodes  made  of  stainless  steel  or  aluminum,  or  a  combination  of  each,  are 
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used  in  the  rf  reactor.  The  discharge  can  be  run  with  both  electrode  surfaces  e3q)Osed 
to  the  discharge  or  a  thin  substrate  may  be  placed  on  the  bottom  electrode  and  aflBxed 
to  the  top  electrode.  The  cell  is  fed  with  a  flow  of  molecular  nitrogen  and  maintained 
in  the  pressure  range  of  1  -  5  Torr.  For  temporal  analysis  of  the  evolution  and  decay 
of  the  atomic  nitrogen  concentration,  the  rf  discharge  is  pulsed  for  a  duration  of  2  ms 
at  a  rep.  rate  of  10  Hz. 

The  atomic  nitrogen  concentration  in  the  interelectrode  space  was  monitored  by 
the  TALIF(207)  technique  detailed  in  Chapters  3  and  4.  The  temporal  evolution  and 
decay  of  the  N  atom  density  was  obtained  by  triggering  the  laser  pulse  to  occur  at 
speciflc  times  during  the  pulse  discharge  or  decay  cycle.  The  reactor  translation  stage 
was  used  to  position  the  cell  with  respect  to  the  laser  focus  to  obtain  the  temporal 
behavior  of  the  atomic  nitrogen  at  various  locations  between  the  electrode  surfaces. 
This  data  was  compared  to  the  output  of  a  numerical  model  whidi  simulates  volume 
and  surface  processes  occurring  within  the  reactor. 

7.3  Numerical  model  for  atomic  nitrogen  in  rf  cell 

A  two-dimensional  numerical  model  was  developed  to  simulate  the  volume  and  surface 
processes  occurring  with  atomic  nitrogen  in  the  parallel  plate  rf  reactor.  The  model 
accounts  for  regional  production,  diflFiision,  volume  loss,  and  loss  of  atomic  nitrogen 
on  multiple  surfaces.  Time  varying  distributions  of  N  atoms  can  be  obtained  and 
compared  to  experimental  data. 

The  motivation  behind  developing  this  model  was  to  improve  on  the  analytical 
solution  for  an  infinite  parallel  plate  cell  by  including  the  effects  of  a  finite  radial 
dimension  of  the  electrode  for  this  experimental  cell.  These  effects  include  radial 
diffusion  of  particles  beyond  the  interelectrode  space  and  any  boundary  effects  from 
the  low  loss  Pyrex  surface  which  eventually  bound  the  outer  volume  of  the  cell. 
The  two  dimensional  numerical  model  was  developed  in  cylindrical  coordinates  with 
azimuthal  symmetry  so  that  the  interelectrode  space  as  well  as  the  outer  volume 
bounded  by  Pyrex  can  be  accounted  for  along  with  the  different  surfaces  and  their 
respective  loss  properties. 
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A  complete  description  of  the  numerical  model  can  be  found  in  Appendix  A,  with 
details  on  how  the  processes  of  N  atom  production,  diffusion,  volume  loss,  and  surface 
loss  were  represented  in  the  code.  The  basic  grid  structure  for  the  two-dimensional 
model  is  drawn  in  Fig.  36  with  the  electrode  and  wall  surfaces  shown.  The  actual 
model  uses  a  grid  spacing  of  1  mm  in  each  dimension,  which  is  finer  than  what  is 
drawn.  The  grid  elements,  contain  the  N  atom  density  as  a  function  of  position 
and  time,  where  i  is  the  z  index,  j  is  the  p  index  and  k  is  the  time  index.  The  grid 
elements  are  updated  as  the  program  steps  through  a  series  of  operations  on  the  grid 
elements  representing  the  production,  diffusion  and  loss  processes  that  occur  during  a 
short  time  interval.  Each  update  represents  a  time  interval  on  the  order  of  10  //s,  with 
the  exact  value  depending  on  the  diflnsion  coefficient  as  described  in  Appendix  A. 
When  modeling  a  pulsed  rf  discharge  and  subsequent  N  atom  decay,  the  production  is 
set  to  zero  after  a  certain  pulse  duration  time  while  the  diffusion  and  loss  operations 
continue  beyond  the  end  of  the  pulse.  Most  of  the  N  atom  decay  data  has  been 
recorded  by  TALIF  at  a  position  near  one  of  the  electrodes,  as  shown  in  Fig.  36.  The 
TALIF  decay  data  is  fitted  with  the  decay  of  the  grid  element  at  that  same  position, 
where  the  most  influential  input  parameter  is  the  surface  recombination  rate,  7,  at 
the  electrode  surface. 

Two  example  outputs  of  the  model  are  shown  in  Figs.  37  and  38,  where  the 
full  two-dimensional  N  atom  spatial  distributions  for  a  2  ms  rf  discharge  pulse  at  3 
Torr  are  displayed  for  difierent  times.  Selected  spatial  regions  of  the  model  output 
will  be  compared  with  experimental  data  in  Section  7.4.1.  Figure  37  displays  the 
distribution  immediately  after  a  2  ms  rf  pulse,  while  Fig.  38  shows  the  distribution 
20  ms  after  the  2  ms  rf  pulse  was  tumed-off.  The  input  parameters  for  this  example 
included  7=0.6%  for  electrode  surface  loss  and  7=0.01%  for  the  loss  at  the  outer  pyrex 
boundary.  The  diffusion  coefficient  for  N  atom  diffusion  through  N2  has  a  value  220 
cm^/s  at  1  Torr[71]  and  is  therefore  set  to  73  cm^/s  at  3  Torr  in  this  case.  The  N  atom 
production  zones  were  modeled  as  regions  with  a  radius  the  same  as  the  electrode 
(2.5  cm)  and  boundaries  in  the  z-direction  at  3  mm  and  6  mm  in  firom  ea.di  electrode. 
The  plots  in  Figs.  37  and  38  show  that  the  spatial  distribution  spreads  radially  into 
the  outer  volume  of  the  reactor  with  time.  This  indicates  that  the  temporal  decay  of 
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Figure  36:  Drawing  of  grid  structure  for  numerical  model  where  each  element, 
contains  the  N  atom  density  at  that  point  and  time,  with  i=z  index,  j=p  index, 
k=time  index 
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atomic  nitrogen  is  affected  by  loss  at  the  surfaces  and,  to  a  certain  degree,  by  diffusion 
from  the  production  volume  into  the  outer  volume  of  the  reactor. 

7.4  Results  and  Discussion 

7.4.1  N  atom  recombination  rates  at  reactor  surfaces 

The  N  atom  recombination  rate,  7,  is  determined  for  several  surface  materials  exposed 
to  a  pulsed  rf  discharge.  The  rate  of  decay  of  the  N  atom  density  was  measured  near 
aluminum  and  stainless  steel  electrodes  as  well  as  silicon  and  boron  nitride  wafers 
that  covered  the  electrode  surface.  The  experimental  decay  data  was  fit  by  a  com¬ 
puter  model  whidi  was  a  strong  function  of  7  near  the  surface.  An  example  of  the 
consistency  between  the  model  and  data  is  shown  in  Fig.  39  for  a  2  ms  pulsed  rf 
disdiarge  at  5  Torr  of  N2  with  stainless  steel  electrodes.  Each  set  of  data  is  a  spatial 
scan  along  the  z-direction  between  the  center  of  the  electrodes  at  p=0.  Experimental 
data  points  are  shown  for  three  different  times  in  the  post-discharge  with  the  mesh 
representing  the  output  of  the  model.  The  pulsed  N  atom  production  rate  and  re¬ 
combination  rate,  7,  were  adjusted  in  the  model  to  fit  the  spatial  distribution  data  at 
0.2, 1  and  3  ms  after  the  discharge.  The  resulting  value  of  7=0.48%  was  determined 
for  this  case. 

A  more  eflicient  method  of  determining  7  was  to  fit  only  the  decay  of  a  single 
spatial  point  near  the  surface.  The  decay  rate  of  the  N  atom  density  was  studied 
near  each  surface  by  taking  TALIF  measurements  at  a  point  2  mm  from  one  of 
the  electrodes  and  centered  with  the  electrode.  This  point  corresponds  to  the  p=0 
coordinate  and  would  be  at  either  z=2  mm  or  20  mm  from  the  bottom  electrode, 
depending  on  whether  the  bottom  or  top  electrode  was  being  studied.  The  atomic 
nitrogen  decay  at  these  coordinates  during  the  post-discharge  was  a  strong  function  of 
the  recombination  loss  rate  to  the  nearby  surface.  A  steady  decay  rate  was  established 
at  this  point  near  the  surface  even  as  early  as  1  ms  after  the  discharge.  This  is  in 
contrast  to  the  center  coordinate  halfway  between  the  electrodes,  shown  as  z=ll  in 
Fig.  39,  where  the  N  atom  density  actually  increased  for  several  milliseconds  due  to 
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Figure  37:  Example  output  of  numerical  model  for  N  atom  spatial  distribution  i 
3  Torr  N2  rf  pulsed  discharge  immediately  after  a  2  ms  rf  pulse 
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diSiision  from  production  regions.  It  was  thus  found  to  be  most  effective  to  solve  for 
7  by  fitting  the  data  taken  at  2  mm  from  the  surface. 

The  temporal  decay  of  the  N  atom  density  near  each  of  the  surfaces  of  interest  is 
shown  in  Fig.  40.  Each  plot  shows  the  N  atom  decay  after  a  2  ms  rf  pulse  has  been 
applied  at  5  Torr  of  N2.  The  data  points  are  the  experimental  measurements  and 
the  lines  represent  the  least  squares  fit  of  the  numerical  model  with  the  associated 
recombination  rates  of  the  fit,  7,  also  shown.  As  seen  in  Fig.  40,  all  the  surface 
materials  tested  at  5  Torr  had  a  7  less  than  1%,  with  stainless  steel  having  the  largest 
surface  recombination  rate  and  boron  nitride  having  the  smallest. 

Each  of  the  surfaces  were  found  to  have  recombination  rates  that  were  pressure 
dependent.  This  is  demonstrated  in  Figs.  41  through  44,  which  show  the  N  atom 
temporal  decays  at  various  discharge  pressures  near  each  surface.  The  complete  set  of 
N  atom  recombination  rates  for  each  surface  at  1,  3,  and  5  Torr  are  compiled  in  Table 
10.  All  of  the  surface  materials  were  found  to  have  the  similar  trend  of  an  increasing 
N  atom  surface  recombination  at  lower  pressures. 


Surface  Material 

5  Torr  3  Torr  1  Torr 

Boron  Nitride 
Silicon 

Aluminum 

Stainless  Steel 

0.02%  0.08%  0.09% 

0.05%  0.16%  0.26% 

0.10%  0.18%  0.28% 

0.48%  0.63%  0.75% 

Table  10:  N  atom  recombination  rates  in  N2  pulsed  rf  discharge 


The  trends  of  this  data  can  be  compared  to  previous  studies,  where  N  atom 
recombination  data  was  obtained  for  various  materials  exposed  to  the  late  afterglow  in 
a  fiowing  post-discharge  [22, 25,  71, 72, 73].  Previous  results  for  N  atom  recombination 
rates  on  silver  and  copper[73]  and  sputtered  molybdenum  [71]  in  a  nitrogen  afterglow 
are  particularly  relevant.  The  values  of  7  for  these  metals  fall  in  a  range  between  the 
7’s  for  aluminum  and  stainless  steel  in  Table  10.  More  importantly,  the  silver,  copper, 
and  sputtered  molybdenum  data  showed  a  decreasing  trend  in  7  as  the  pressure 
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Figure  40:  N  atom  temporal  decay  in  5  Torr  N2  pulsed  rf  discharge  with  corresponding 
N  atom  recombination  rates,  7,  shown  for  boron  nitride,  silicon,  aluminum,  and 
stainless  steel 
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Figure  41:  N  atom  temporal  decay  2  mm  from  boron  nitride  surface  after  a  2  ms  rf 
discharge  pulse  at  5  Torr,  3  Torr  and  1  Torr 
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Figure  42:  N  atom  temporal  decay  2  mm  from  silicon  wafer  surface  after  a  2  ms  rf 
discharge  pulse  at  5  Torr,  3  Torr  and  1  Torr 
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Figure  43:  N  atom  temporal  decay  2  mm  from  aluminum  surface  after  a  2  ms 
discharge  pulse  at  5  Torr,  3  Torr  and  1  Torr 
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was  increased  that  is  similar  to  the  trend  for  all  the  surfaces  in  Table  10.  The 
suggested  theory  for  this  pressure  dependence[71 ,  73]  involves  preliminary  adsorption 
of  atomic  nitrogen  occurring  on  the  surface  with  a  protective  molecular  nitrogen  layer 
forming  over  the  saturated  atomic  layer.  The  molecular  nitrogen  layer  would  have 
a  greater  fractional  surface  coverage  at  higher  pressure.  With  this  protective  layer 
preventing  gas  phase  atomic  nitrogen  from  reaching  the  catalytic  surface,  the  atomic 
recombination  rate  will  be  lower  in  the  higher  pressure  environment.  The  pressure 
dependence  of  the  7  values  in  Table  10  for  each  surface  material  is  consistent  with 
the  trend  of  the  previous  surfaces  tested  and  the  suggested  theory  of  a  protective 
nitrogen  layer.  Other  factors,  though,  should  be  considered  in  the  theory,  such  as  an 
increased  atomic  density  competing  for  recombination  sites  at  higher  pressures  or  an 
increased  ion  bombardment  flux  at  lower  pressure  that  may  produce  recombination 
sites. 

With  respect  to  the  design  of  a  discharge  reactor  for  an  atomic  nitrogen  source,  it 
is  clear  in  terms  of  N  atom  loss  that  aluminum  is  a  better  dioice  than  stainless  steel 
for  a  metal  surface  material  that  is  positioned  near  the  disdiaxge.  It  is  also  seen  that 
the  surface  loss  on  aluminum  does  not  greatly  exceed  that  of  silicon,  which  represents 
a  common  processing  material  that  would  be  placed  in  the  reactor.  Therefore,  a 
moderate  amount  of  exposed  aluminum  near  the  discharge  volume  where  a  silicon 
wafer  is  being  processed  would  not  cause  a  large  increase  in  N  atom  loss  rate  above 
that  due  to  the  silicon.  In  contrast,  the  recombination  loss  for  stainless  steel  was 
found  to  be  much  greater  than  silicon.  Thus  stainless  steel  is  more  likely  to  alter  the 
balance  between  N  atom  loss  and  production  when  positioned  near  a  discharge  with 
a  silicon  wafer  loaded.  The  boron  nitride  surface  displayed  the  lowest  7  values  of 
all  the  materials  tested.  An  optimal  N  atom  source,  therefore,  may  be  developed  by 
replacing  or  covering  surfaces  exposed  to  the  discharge  with  boron  nitride.  Although 
boron  nitride  wafers  covered  the  electrodes  in  this  experiment,  commercial  processes 
exist  where  pyrolytic  boron  nitride  can  be  deposited  and  form  coatings  on  objects  of 
various  shape[74]. 
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7.4.2  N  atom  response  to  oxygen  impurities  in  rf  cell 

Small  amounts  of  O2  were  introduced  into  the  N2  rf  parallel  plate  discharge  to  simu¬ 
late  the  infiltration  of  oxygen  impurities  into  the  system.  The  effects  on  the  atomic 
nitrogen  density  were  obvious  with  only  a  few  tenths  of  a  percent  O2  in  N2.  Figure  45 
shows  spatial  scans  of  the  N  atom  density  along  the  z-axis  between  the  center  of  the 
electrodes  for  a  continuous  rf  discharge  with  pure  N2  and  with  0.2%  O2  in  N2.  The 
spatial  distribution  with  0.2%  O2  added  shows  that  the  N  atom  density  was  depleted 
to  approximately  one-third  of  its  value  in  pure  N2.  The  differences  in  the  spatial 
profiles  indicate  that  the  introduction  of  O2  either  caused  a  significant  decrease  in 
the  N  atom  production  or  caused  an  increase  in  the  volume  or  surface  loss  of  N  atoms. 
The  enhanced  density  around  the  sheath  production  regions  of  the  0.2%  O2  scan  is 
an  indication  that  an  increased  loss  rate  is  limiting  the  N  atom  build-up  by  diffusion. 

Oxygen  was  also  added  in  various  small  amounts  to  the  parallel  plate  discharge 
operated  with  a  2  ms  rf  pulse  at  10  Hz.  The  temporal  decay  for  each  mixture  was 
recorded  for  post-discharge  times  up  to  20  ms  after  the  discharge  was  tumed-off. 
Plots  of  temporal  decay  are  shown  in  Fig.  46  for  pure  N2,  0.1%  O2,  0.2%  O2,  and 
0.4%  O2  in  N2.  The  data  demonstrates  the  N  atom  dec»y  becoming  more  rapid  for 
higher  concentrations  of  the  O2  impurity.  This  is  consistent  with  an  increased  loss 
rate  of  N  atoms,  either  at  the  reactor  surfaces  or  within  the  volume  of  the  cell  as  O2 
is  introduced. 

The  volume  loss  rate  of  atomic  nitrogen  would  be  increased  if  the  introduction  of 
O2  in  the  discharge  led  to  formation  of  NO,  which  reacts  quichly  with  N  atoms.  This 
mechanism  must  be  studied  closer,  though,  since  the  formation  of  NO  itself  often 
involves  atomic  nitrogen  as  a  reactant  or  product.  The  major  processes  of  formation 
and  destruction  of  NO  in  an  N2/O2  discharge  are[75] 


N2(v>12)-fO 

NO  -f  N(hot) 

(ki  =  1  X  10“^^cm®/s) 

(7.1) 

N(hot)  +  O2  — > 

NO-t-0 

(ki  =  3  X  10“^^cm®/s) 

(7.2) 

N  +  m 

N2(v  =  3,4)  +  0 

(k2  =  1.8  X  10~^^cm^/s). 

(7.3) 

A  primary  reaction  for  production  of  NO  is  Eq.  7.1  since  the  vibrational  kinetics 
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Figure  45:  N  atom  spatial  profiles  along  the  central  z-axis  between  the  electrode 
centers  for  a  continuous  rf  discharge  with  pure  N2  and  0.2%  O2  in  N2. 
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Figure  46:  Temporal  decay  of  N  atom  after  a  2  ms  rf  discharge  pulse  with  pure  N- 
0.1%  O2,  0.2%  O2,  and  0.4%  O2  in  Nj. 


of  the  discharge  cause  the  N2(v>  12)  states  to  be  well  populated[75].  The  second 
reaction  that  produces  NO,  Eq.  7.2,  has  an  energy  threshold  of  0.271  eV.  This  would 
require  a  significant  concentration  of  “hot”  N  atoms,  represented  by  N(hot),  with 
enough  translational  energy  to  overcome  the  0.271  eV  threshold  for  this  process  to 
contribute  to  the  production  of  NO.  The  hot  N  atoms  come  firom  the  reaction  in  Eq. 
7.1,  where  approximately  one-half  of  the  the  additional  energy  firom  N2  vibrational 
states  with  v>  12  that  are  involved  in  the  reaction  is  channeled  into  translational 
energy  of  the  N  atoms  produced[75]. 

A  numerical  model  of  the  kinetics  represented  in  Eqs.  7.1  through  7.3  was  devel¬ 
oped  to  demonstrate  that  the  depletion  of  atomic  nitrogen  density  with  the  addition 
of  O2  in  Fig.  45  is  consistent  with  an  increase  in  the  volume  loss  rate  within  the 
cell.  The  differential  rate  equations  for  N  atoms,  N2(v>12),  N(hot)  atoms,  NO,  and  0 
atoms  are 

^  =  PN  +  [N2(v>12)][0]ki 
clt 

-[N][NO]k2  -  [N,io„][0,]k5  -  ^[N]  (7.4) 

=  Pn,(v)  -  [N,{v  >  12)][01k, 

-5^[N,(v  >  12)J  (7.5) 

^N2(v) 

^2^  =  +[N2(v  >  12))(0]ki  -  [N(ko.)l[NO|k2 

-[N(bo.)ll02]k3  -  (7.6) 

=  +[N,(v  >  12)l[0]k,  -  [N][N01k2 

+lN(h«t)l[02lk3  (7.7) 

-[N2(v  >  12)]10]ki  +  lN][NO]k2 
+[N(b«,1102]k3  -  ^[0]  (7.8) 


dt 
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where  Pjv  and  Pn2(v)  are  the  point  production  rates  of  N  atoms  and  N2(v>12)  re¬ 
spectively  due  to  discharge  processes  not  included  in  Eqs.  7.1  throu^  7.3,  sudi  as 
direct  electron  impact  and  N2  vibrational-vibrational  interaction.  The  diffusion  loss 
rate  to  the  walls  are  given  by  where  X=N,  0,  and  N2(v). 

A  numerical  model  using  Euler’s  method  has  been  applied  to  the  five  differential 
rate  equations.  The  differential  rate  equations  for  each  species,  Mj,  where  i=l,...,5, 
are  of  the  form 

^  =  (7.9) 

with  the  formula  for  the  Euler  method  being 


Mi,k4.i  —  Mi^k  +  hfi(tk,  Mi,k>  Ms^k) 


(7.10) 


which  advances  a  solution  by  a  time  interval  of  h  from  time  tk  to  tk+i=tk+h.  A  time 
stepsize  interval  value  of  h=10“®  seconds  was  used  here  which  provided  the  necessary 
stability  and  accuracy  for  the  numerical  analysis. 

Since  the  data  in  Fig.  45  was  extracted  from  a  continuous  rf  disdiarge,  the 
densities  of  the  discharge  products  were  assumed  to  have  reached  a  steady-state  value. 
The  numerical  model  was  programmed,  therefore,  to  continue  to  advance  the  solution 
in  time  until  steady-state  conditions  were  reached  for  all  five  species.  The  density 
results  for  the  five  species  were  then  checked  for  accuracy  by  inserting  them  back  into  a 
set  of  analytical  equations  representing  the  steady-state  solutions  for  dMi/dt=0.  The 
analytical  steady  state  solutions  to  the  differential  rate  equations  for  each  species  are 

rxTi  _  Pn  +  [N2(v  >  12)][0]k,  -  [Nn,ot)][02lk3  . 

^  ^  [NO]k2  -h  ^ 


[N2(v  >  12)]  = 


[Ojki  -t- 


Pn2(v) 

r.  _L  °N2(v) 


^  =  [N.(v  >  12)l[0]k, 

‘  [Ozlks  +  [NO]k2  +  ^ 


(NO]  = 


[N2(y  >  12)l[0]ki  +  [N(h„)l[02lk3 


rpi  [N][N0]k2  +  [Nti^)][02]k 

i  J  rKT.  ^  1  o\ii,  1  Do 


[N2(y  >  12)]ki  +  W 


(7.12) 


(7.13) 


(7.14) 

(7.15) 
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In  order  to  execute  the  numerical  model,  values  were  needed  for  the  diffusion  loss 
terms,  and  along  with  the  production  terms,  Pn  and  Pn2(v)-  These 

Aff  Aq  ^N2(v) 

values  were  set  corresponding  to  the  conditions  which  produced  the  data  in  Fig.  45 
whidi  was  a  5  Torr  Nj  discharge  in  the  rf  cell  with  [O2]/[N2+O2]=0%  or  0.2%.  The 
diffusion  loss  terms  for  each  of  the  species  was  determined  either  from  experimental 
observation  of  decay  by  diffusion  loss  or  from  a  theoretical  calculation  of  the  funda¬ 
mental  diffusion  mode  within  the  rf  cell.  For  atomic  nitrogen,  the  diffusion  loss  term 
was  determined  from  the  pulsed  rf  experimental  data  for  %O2=0  in  Fig.  46.  In  the 
absence  of  oxygen,  the  differential  rate  equation  for  [N]  after  the  rf  pulse  simplifies  to 

=  -25  ■  (7.16) 

[N]  dt  A^'  '  ^ 

The  diffusion  rate  with  loss  at  the  surfaces  is  thereby  isolated  as  the  only  major  loss 
term.  The  value  of  W  is  equal  to  the  N  atom  exponential  decay  rate  of  24  s~^  as 
determined  from  Fig.  46. 

The  atomic  oxygen  diffusion  loss  rate  in  a  0.2%  O2/N2  mixture,  though,  can  not 
be  so  easily  isolated  from  the  volume  losses.  In  the  absence  of  experimental  data  that 
directly  compares  the  surface  recombination  rates  of  N  and  O  atoms,  a  reasonable 
assumption  could  be  that  Dq/Aq  «  Dn/A^.  To  check  this  assumption,  a  theoretical 
diffusion  loss  rate  for  Dq/Aq  was  calculated  by  modeling  the  rf  cell  as  an  infinite 
parallel  plate.  The  detailed  diffusion  theory  for  the  rf  cell  presented  in  Appendix  A 
leads  to  the  transcendental  equation 


V  J  TT  .TTZo, 

4D(1-|)  ”  z. 


(7.17) 


where  v  is  the  particle  velocity,  7  is  the  surface  recombination  coefficient,  Zo  is  the 
electrode  spacing  (2.2  cm),  and  Zj  is  the  extrapolation  parameter  representing  an 
imaginary  position  b^ond  the  electrode  where  the  particle  density  goes  to  zero.  The 
parameter  z,  for  an  infinite  plate  is  related  to  the  fundamental  diffusion  length.  A, 
by  considering  only  the  z  dimension  in  Eq.  A.21  to  give 

By  combining  equations  7.17  and  7.18  and  using  Do=40  cm^/s  in  5  Torr  of  N2[75],  the 
value  of  Do/Ao=24  s~^  corresponds  to  7o=0.25%.  This  70  value  is  very  reasonable. 


135 


considering  a  value  of  70  in  a  pure  O2  discharge  at  3  Torr  has  been  measured  at  a 
stainless  steel  surface  to  be  0.4%[12]  while  70  was  tending  to  decrease  with  increasing 
pressure.  With  these  considerations,  the  value  of  Do/A^=24  s~^  was  applied  to  the 
numerical  model. 

The  assignment  of  Dn2(v)/An2(v)  involved  a  significant  uncertainty  due  to  the  lack 
of  experimental  measurements  of  N2(X,v)  difiusion  and  wall  loss.  Using  a  theoreti¬ 
cal  model  that  assumed  a  single-quantum  wall  de-activation  or  N2(X,v)^^N2(X,v-l), 
Loureiro[47]  found  that  a  value  of  7n2(v)=0-7%,  which  was  independent  of  v,  on  molyb- 
denic  glass  resulted  in  the  best  fit  of  experimental  dissociation  rates  verses  vibrational 
temperature  at  5.6  Torr.  The  value  of  7n2(v)=0.7%  was,  therefore,  chosen  as  an  ad¬ 
equate  approximation  for  N2(X,v)  wall  loss  at  5  Torr  in  this  model.  Again  using 
equations  7.17  and  7.18  with  Dn2(v)=10  cm^/s  at  5  Torr[65]  and  7n2(v)=0.7%,  a  value 
of  Dn2(v)/Anj(v)=17  s“^  was  determined. 

The  assignments  of  the  production  terms  were  based  on  a  measured  steady- 
state  N  atom  density  of  1.3x10^®  cm"®  for  the  given  conditions  and  a  density  of 
[N2(X,v>12)]=9x10^^  cm“®  derived  from  observations  in  the  rf  cell  under  similar 
conditions  in  Section  6.4.2.  When  these  values  were  inserted  into  Eqs.  7.11  and  7.12 
in  the  absence  of  oxygen,  the  values  of  Pn=3.12x10^^  cm®/s  and  Pn2(v)  =1.5x10^^ 
cm®/s  resulted  for  the  pure  N2  discharge.  It  was  assumed  that  these  production 
terms  were  not  altered  by  the  addition  of  0.2%  O2  to  the  N2  discharge  and  that  the 
changes  in  N  atom  density  were  due  to  the  kinetics  of  Eqs.  7.1  though  7.3. 

The  resulting  densities  of  the  five  species  that  represent  a  steady  state  solution 
to  the  differential  rate  equations  for  a  5  Torr  continuous  rf  discharge  are  shown  in 
Table  11  for  both  0%  O2  and  0.2%  O2.  As  a  check  for  accuracy,  the  results  from 
the  numerical  model  for  0.2%  O2  were  then  inserted  back  into  the  five  steady-state 
analytical  solutions  and  the  resulting  densities  each  agreed  to  within  0.5%  of  the 
numerical  solutions.  The  results  of  the  kinetic  model  in  Table  11  demonstrate  a 
depletion  in  atomic  nitrogen  concentration  as  0.2%  O2  is  added  to  the  N2  discharge 
and  the  increase  loss  rate  can  be  accounted  for  by  the  volume  processes  in  the  rf  cell. 
The  N  atom  depletion  in  the  model  agrees  well  with  the  experimental  observations 
that  are  plotted  in  Fig.  46. 
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%02 

[N] 

1N2(X,v>12)] 

[N(hot)] 

[NO] 

[0] 

0% 

0.2% 

1.3x10^'’ 

3.77x10“ 

9x10“ 

1.37x10“ 

1.84x10“ 

3.49x10“ 

9.22x10“ 

Table  11:  Steady-state  density  values  in  cm  ®  for  numerical  model  of  N2/O2  rf  dis¬ 
charge  at  5  Torr 


To  summarize  this  volume  reaction  series,  the  hot  N  atoms  created  in  Eq.  7.1 
react  readily  with  O2  (Eq.  7.2)  which  produces  NO.  The  NO  produced  in  Eqs.  7.1 
and  7.2  react  with  the  remaining  atomic  nitrogen  in  the  discharge  which  results  in  a 
significant  N  atom  volume  loss  mechanism. 

7.5  Conclusion 

Surface  and  volume  loss  processes  have  been  investigated  for  a  N2  rf  parallel  plate 
discharge.  Atomic  nitrogen  recombination  rates  were  determined  for  several  surface 
materials  that  were  positioned  near  the  discharge  volume.  Stainless  steel,  aluminum, 
silicon  and  boron  nitride  were  all  found  to  have  values  of  7  that  had  an  inverse  pressure 
dependence.  This  pressure  dependence  is  consistent  with  the  theory  of  a  molecular 
nitrogen  layer  forming  over  an  adsorbed  atomic  nitrogen  layer,  where  the  fractional 
coverage  of  the  molecular  layer  increases  with  pressure.  The  arrival  rate  of  N  atoms 
reaching  the  surface  will  be  lower  at  higher  pressures  and  therefore  7  will  decrease 
with  increasing  pressure.  Other  possible  effects  with  increasing  pressure  that  could 
decrease  the  N  atom  recombination  rate  include  an  increase  in  N  atom  competition  for 
recombination  sites  or  a  decrease  in  the  generation  of  sites  from  ion  bombardment.  It 
was  also  found  that  almninum  had  a  surface  loss  rate  of  3  to  5  times  less  than  that  of 
stainless  steel.  This  indicates  that  aluminum  should  be  considered  over  stainless  steel 
for  a  reactor  material  that  will  be  exposed  to  the  discharge  products  in  an  N  atom 
source.  Boron  nitride  was  found  to  have  the  lowest  N  atom  recombination  rate  of  the 
materials  studied  under  eadi  discharge  condition.  This  low  recombination  rate  could 
be  utilized  in  an  N  atom  source  by  coating  exposed  reactor  surfaces  with  pyrolytic 
boron  nitride. 
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The  effects  of  the  infiltration  of  small  amounts  of  oxygen  into  the  rf  discharge 
were  demonstrated.  A  controlled  addition  of  O2  to  the  N2  discharge  was  found  to 
cause  a  significant  reduction  in  the  N  atom  density.  The  O2  was  then  added  to  a 
pulsed  rf  discharge  and  the  N  atom  decay  was  observed  in  the  post-discharge.  It  was 
determined  from  the  N  atom  decay  that  the  O2  was  responsible  for  the  onset  of  a 
substantial  volume  or  surface  loss  medianism  that  was  reducing  the  net  concentration 
of  atomic  nitrogen  in  the  discharge.  A  series  of  reactions  involving  the  formation  of 
NO  and  its  subsequent  reaction  with  atomic  nitrogen  were  introduced  as  a  volume 
loss  mechanism.  A  set  of  kinetic  equations  involving  the  volume  reactions  were  solved 
numerically  with  the  solution  being  consistent  with  the  reduction  of  N  atom  density 
with  the  introduction  of  oxygen  impurities. 
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CHAPTER  8 


Conclusion 


Detection  of  atomic  nitrogen  by  two-photon  laser  induced  fluorescence  (TALIF)  has 
been  used  to  study  the  ground  N(2p®)^S3/2  atomic  state  in  a  variety  of  discharge 
environments.  Our  flndings  not  only  demonstrated  the  utility  of  N  atom  TALIF  in  a 
variety  of  gas  phase  kinetic  studies,  but  also  promoted  an  alternative  excitation  and 
detection  scheme  for  atomic  nitrogen  TALIF. 

An  alternative  N  atom  TALIF  scheme  was  found  to  have  many  advantages  over  a 
TALIF  scheme  that  was  previously  believed  to  be  the  optimal  technique  for  atomic 
nitrogen  detection.  Utilizing  laser  excitation  at  207  nm,  this  alternative  scheme  was 
shown  to  produce,  in  general,  a  greater  TALIF  signal  strength  than  the  traditional 
scheme,  which  utilizes  laser  excitation  at  211  nm.  This  increased  signal  strength  was 
found  to  be  partly  due  to  a  slightly  greater  two-photon  absorption  rate,  but  primarily 
due  to  a  greatly  reduced  coUisional  quenching  rate  by  N2.  Detection  of  the  alternative 
TALIF  signal,  with  its  near  visible  fluorescence  wavelength  of  750  nm  was  also  found 
to  be  more  convenient  for  PMT  detection  than  the  870  nm  in&ared  fluorescence  of 
the  previous  scheme. 

Photodissociation  of  N2O  by  207  nm  laser  light  and  subsequent  TALIF  detection  of 
N  atom  photofragments  was  investigated  for  the  potential  application  of  determining 
the  spatial  view  factor  within  a  discharge  cell.  The  results  were  not  encouraging  for 
this  application  due  to  a  low  N  atom  photodissociation  yield  and,  more  importantly, 
severe  coUisional  quenching  of  the  upper  TALIF  state.  This  discovery,  though,  led 
to  questions  of  why  was  the  TALIF  signal  being  quenched  so  much  more  than  in 
the  N2  discharge  environment.  Additional  TALIF  measurements  on  N  atoms  in  a  N2 
flowing  discharge  afterglow  showed  that  N2O  addition  did  not  quench  the  upper  state 
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with  the  same  rate  as  that  observed  in  the  photodissociation  results.  A  study  of  the 
N2O  photochemistry  and  the  dynamics  of  the  photofragments  showed  the  possibility 
that  “hot”  N  atom  fragments  were  the  reason  for  the  discrepancy  in  the  quenching 
data.  This  explanation  of  the  greater  TALIF  signal  quenching  duiring  N2O  photolysis 
required  the  N  atom  fragment  to  have  a  kinetic  energy  of  0.6  eV,  which  was  found  to 
be  energetically  allowed  for  the  photoreaction.  We  did  not  obtain  data  to  confirm  the 
hot  N  atom  photofragment  theory,  but  instead  outlined  some  promising  suggested 
experiments. 

A  quantitative  model  capable  of  predicting  the  absolute  concentrations  of  several 
discharge  species  was  developed  and  then  validated  by  extensive  experimental  mea¬ 
surements  in  an  rf  discharge  cell.  The  temporal  and  spatial  resolution  of  the  N  atom 
detection  along  with  the  capability  to  determine  the  absolute  concentrations  through 
titration  measurements,  made  the  TALIF  measurements  ideal  data  from  which  to 
build  the  kinetic  model.  This  model  combined  absolute  N  atom  densities  with  rela¬ 
tive  LIF  and  discharge  emission  measurements  of  various  molecular  nitrogen  species 
and  calculated  the  absolute  densities  for  several  molecular  species.  Absolute  densities 
of  some  species,  such  as  N2(X,v>5),  were  determined  from  the  model  in  spite  of  an 
inability  to  detect  the  species  by  optical  diagnostics.  The  comprehensiveness  of  this 
model  and  its  ability  to  replicate  a  wide  variety  of  discharge  conditions  indicate  a 
capacity  to  serve  as  a  predictive  model  for  various  discharge  reactors. 

The  TALIF  technique  was  used  in  investigating  the  effect  of  N2/H2  gas  mixtures 
in  high  frequency  discharges.  Small  additions  of  H2  to  an  N2  microwave  discharge 
were  found  to  alter  the  heavy  particle  kinetics  of  the  system  and  a  modest  peak  in 
N  atom  density  occurred  at  a  few  tenths  of  a  percent  addition  of  H2  in  N2.  The 
enhancement  and  decline  of  atomic  nitrogen  density  with  increasing  %H2  along  with 
the  observed  behavior  of  several  molecular  species  were  shown  to  be  consistent  with 
the  changes  in  heavy  particle  kinetics  for  small  additions  of  H2  in  N2. 

The  TALIF  spatial  and  temporal  resolution  was  also  used  in  determining  the  N 
atom  surface  recombination  rates  on  various  materials  within  the  rf  reactor.  Data  on 
surface  recombination  rates,  such  as  the  substantially  lower  N  atom  recombination 
rate  on  aluminum  as  compared  to  stainless  steel,  can  be  applied  in  the  design  of  an 


140 


discharge  cell.  The  atomic  nitrogen  recombination  rate  on  boron  nitride  was  found 
to  be  exceedingly  small  with  the  potential  use  as  a  low  loss  coating  in  an  N  atom 
source  reactor.  Likewise,  atomic  nitrogen  density  measurements  quantified  the  extent 
of  N  atom  depletion  that  occurs  when  a  small  air  leak,  virtual  or  otherwise,  is  left 
undiecked.  Oxygen  added  to  N2  on  the  order  of  0.1%  was  found  to  greatly  reduce  the 
steady  state  N  atom  concentration  and  accelerate  the  volume  loss  processes  within  the 
reactor.  In  conclusion,  N  atom  TALIF  has  proven  to  be  a  responsive  and  dependable 
tool  in  the  investigation  of  a  variety  of  low  pressure  discharge  situations.  Although 
the  N2  discharge  chemistry  as  a  whole  still  retains  much  mystery,  the  behavior  of  the 
ground  state  atomic  nitrogen  can  now  be  clearly  observed. 
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APPENDIX  A 


Numerical  model  of  N  atom  behavior  in  rf  cell 


A.l  Description  of  model 

A  two  dimensional  array,  representing  the  z  and  p  coordinates  in  cylindrical  geome¬ 
try,  is  described  here  in  which  each  array  element  contains  the  value  of  the  atomic 
nitrogen  density  at  that  respective  position  within  the  reactor.  Azimuthal  symmetry 
is  assumed.  The  array  evolves  in  discrete  steps  of  time  while  considering  production, 
diffusion,  and  loss  processes  during  each  step.  The  array  models  the  volume  within  a 
right  circular  cylinder  with  the  top,  bottom,  and  radial  boundaries  being  Pyrex.  Par¬ 
allel  disk  shaped  electrodes  are  suspended  within  the  volume,  centered  on  the  z-axis 
of  the  cylinder  but  extending  only  partially  in  the  radial  direction.  This  model  is  a 
close  representation  of  the  ^cperimental  cell.  The  electrode  dimensions  and  spacing 
are  accurate  and  the  diffusion  radially  away  ffom  and  around  the  electrodes  is  accu¬ 
rately  represented,  which  are  both  very  important  to  modeling  the  decay  of  radical 
species  within  the  interelectrode  space.  The  cylindrical  Pyrex  outer  boundary  of  the 
model  is  a  weak  approximation  of  the  true  shape  of  the  six  way  Pyrex  cross  with  3 
inch  ports,  but  the  position  of  these  walls  is  much  less  critical  to  the  model.  Placing 
the  model  Pyrex  walls  at  an  effective  radial  distance  from  the  edge  of  the  electrode 
and  a  specified  z  distance  behind  the  electrode  will  allow  the  important  effects  of 
diffusion  beyond  the  electrode  space  to  be  accounted  for  while  including  the  eventual 
long  term  build  up  and  back  diffusion  of  N  atoms  from  the  outer  wall  reflection. 

The  program  steps  through  a  series  of  operations  on  the  array  elements  repre¬ 
senting  processes  that  occur  during  a  short  time  interval.  These  operations  include; 
production,  in  which  atomic  density  is  added  to  a  set  of  elements,  diffusion,  in  which 
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atomic  density  is  exchanged  between  elements  in  accordance  with  the  diffusion  equa¬ 
tion,  wall  loss,  where  the  atomic  density  is  adjusted  at  each  surface  to  account  for 
surface  recombination,  and  volume  loss,  where  atomic  density  is  lost  according  to 
reaction  rates  from  collisions  with  other  particles.  After  eadi  operation  is  complete 
for  that  time  interval,  the  proems  begins  again  on  the  augmented  array  for  the  next 
timp  interval.  The  increments  in  time  are  determined  by  a  stability  criterion,  to  be 
discussed  later,  that  allows  an  accurate  analysis  to  proceed  as  rapidly  as  possible. 

The  atomic  production  operation  involves  the  addition  of  the  N  atom  array  with 
a  production  array  with  values  of  atoms  per  cc  being  produced  within  the  respective 
time  interval.  The  distribution  and  values  of  the  atomic  production  elements  within 
the  array  can  be  matched  to  the  known  production  regions  within  the  discharge  cell. 
It  is  common  within  an  N2  rf  parallel  plate  discharge  to  have  two  pancake  shaped 
sheath  regions  near  the  electrodes  in  which  dissociation  of  N2  and  therefore  N  atom 
production  is  concentrated.  The  programmer  may  also  set  the  production  array 
to  zero  at  any  time  to  simulate  the  disdfciarge  being  switdied  off,  as  in  a  pulsed  rf 
discharge  system. 

The  diffusion  operation  is  governed  by  the  diffusion  equation 

^  =  -DV^N  (A.1) 

dt 


where  N  is  the  atomic  nitrogen  density  and  D  is  the  diffusion  coefficient  of  atomic  ni¬ 
trogen  in  N2.  In  cylindrical  coordinates,  assuming  azimuthal  symmetry,  the  diffusion 
equation  is 

dN  r.n  d  ,  dN.  ,  d^N] 

This  formula  can  be  applied  to  our  temporally  evolving  N  atom  array  by  the  develop¬ 
ment  of  a  finite-difference  numerical  recipe.  The  N  atom  array  is  expressed  as 
where  i=z  index,  j=p  index,  and  k=time  index.  The  spacing  between  the  z  or  p  array 
elements  is  given  by  h  while  At  is  the  time  step  increment  between  array  updates. 
The  numerical  recipe  is  based  on  a  two-point  forward  difference  slope  model  of  the 
time  derivative  and  a  three-point  difference  model  of  the  spatial  second  derivative 
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[27].  The  radial  term  on  the  RHS  of  Eq.  A.2  can  be  expressed  in  numerical  form  as 

id,dN,  _  ,,  1 


1  a  ,  m.  1  d  \  p  ,  j 

pdp^^Tp^  ~ 


"I"  "I"  +  •^iJ+l,Jfe)j 

~  Ijfc)  ~  •  (A. 3) 


The  radial  index,  j,  is  equal  to  p/h  so  that  the  numerical  expression  is  simplified  to 

1  a ,  9iV,  1  ri  ^  1 

pdp^  dp'  ~  L2j  +  {^ij-i,k  -  +  ■^tj+i,ik)J  (A.4) 

and  a  divide  by  zero  error  is  ignored  when  j=0  since  the  term  multiplied,  Ni^+i^k  ~ 
Ni,-i,k,  is  also  zero  because  of  symmetry  around  p=0.  This  is  rewritten  for  later 


convemence  as 


dp^^  dp^  /i2  ^•^*>7.*’  •  (A-5) 


The  z  term  on  the  RHS  of  Eq.  A.2  is  expressed  in  three  point  numerical  form  as 

1  \.r  J 

Q^2  ~ /j2  ~  "f"  A^*-ij,fc)J-  (A.6) 

The  complete  numerical  formula  representation  for  the  difiiision  equation  is 

^  At  ~  ~  Ifi  27^  ~in)~  ^ij,k  +  • 


(A.7) 

Solving  the  expression  for  the  N  atom  array  element  of  the  next  time  interval,  Nij^k+i, 


gives 


=  «'v^(l 


4DAt\  DAtl 


^r  Ar  7,  4L/At\  DAt[^^  1  1  1 

Aij,fe+i  -  Nijj,  1^1  j,fc+Arij+i,fc(l+— )+iVij_i,fc(l-— ) 

(Ais) 

which  essentially  computes  the  effect  of  diffusion  on  the  array  element,  Nijj,,  by 
comparison  with  its  four  nearest  neighbors.  A  stability  criteria  must  be  met  when 
operating  with  this  formula  that 


DAt  1 
“4 


(A.9) 
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In  this  program,  D  and  h  are  set  by  the  programmer  and  At  is  computed  so  that 
the  stability  criteria  factor  is  just  less  than  This  allows  the  program  to  run  as  fast 
as  possible.  The  program  loops  through  and  applies  this  diffusion  operation  to  each 
array  element,  except  those  elements  representing  positions  adjacent  to  a  surface. 
The  operation  on  these  surface  elements  will  be  described  in  the  next  section. 

The  surface  loss  operation  accounts  for  the  loss  of  N  atom  particles  at  a  surface 
by  recombination.  This  operation  proceeds  by  adjusting  the  values  of  array  elements 
adjacent  to  a  surface  according  to  solutions  derived  from  the  boundary  conditions.  A 
general  form  for  the  boundary  condition  from  Chantry  [76]  is 


N\dq)b~  X 


(A.10) 


where  q,  in  this  case,  represents  either  the  z  or  p  dimension  and  the  subscript  b  im¬ 
plies  that  the  quantities  are  evaluated  at  the  boundary.  The  quantity  A  is  a  linear 
extrapolation  length  which  indicates  the  distance  beyond  the  boundary  that  the  den¬ 
sity  AT (5)  of  diffusing  particles  would  extrapolate  to  zero.  A  surface  loss  coefficient, 
7,  is  now  defined  as  the  ratio  of  the  flux  of  particles  lost  at  a  surface  to  the  total 
flux  of  particles  impinging  on  a  surface.  The  linear  extrapolation  length,  A  may  be 
written  as  function  of  the  surface  loss  coefficient  as[76] 


(A.n) 


where  v  is  the  paxticle  velocity  given  by  ^S/cT/tttti  and  e  is  a  coefficient  depending 
weakly  on  the  surface  coefficient,  7.  Chantry  recites  a  possible  range  of  e  from  |  to 
0.71,  with  I  appropriate  for  small  7.  Since  7  is  relatively  small  for  all  of  our  cases 
involving  ground  state  atomic  nitrogen,  we  will  use  e  =  |.  This  makes  the  expression 
for  A 


AD{l-l) 


(A.12) 


and  the  boundary  condition  becomes 

O.-S^y  “■•> 

Since  any  array  element  adjacent  to  a  surface  will  be  the  last  element  in  an  array 
row  or  column,  the  surface  gradient  is  best  represented  by  a  multiple  point  numerical 


f-l  =-■ 

\  dq  )h 
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formula  that  evaluates  the  derivative  at  the  last  point.  An  appropriate  formula  for 
this  situation  is  the  four  point  recipe  [27] 


(dN\ 

1  ■ 

\dqj, 

b  6/1  . 

-  im  +  ISNb-i  -  9Nb.2  +  2Nb-3  . 


(A.14) 


By  combining  Eqs  A.  13  and  A.14,  an  expression  in  one  dimension  for  the  boundary 
array  element,  Nb,  is  found  to  be 


~  ~  +  2Nb-3,) 

3hv'y  +  22D{l-l) 


(A.15) 


The  two  dimensional  expression  used  in  the  numerical  operation  on  a  boundary  array 
element  adjacent  to  a  surface  in  the  z-direction  is 


^  _  2D{1  -  -  m-2j  +  2Nb-sj) 

Shvj  +  22D{1-^) 

and  when  adjacent  to  a  surface  in  the  p  direction 

ff  . .  20(1  -  ?)(18JV<p-i  -  9iVip-;  +  2JVip_3) 
3/iu7  +  22D(1-|) 


(A.16) 


(A.17) 


The  program  applies  this  operation  to  each  surface  element  during  each  time  interval. 

The  volume  loss  operation  accounts  for  loss  of  atomic  nitrogen  through  gas  phase 
reactions.  An  example  of  this  in  a  pure  N2  discharge  is  the  three-body  recombination 
reaction 


N  +  N  +  N2 N2(B)  +  N2.  (A.18) 

The  N2  concentration  is  assumed  to  remain  spatially  uniform  and  constant  during  the 
program  and  is  determined  from  the  total  pressure  provided  by  the  programmer.  The 
loss  of  N  atom  density  in  an  array  element  during  a  certain  time  interval  is  expressed 
as 

{ANij)voi  =  X  [Aa]  x  k  x  At  (A.19) 

where  k  is  the  recombination  rate  for  Eq.  A.18  with  a  value  of[48]  fc=4.4  x  10"^^  cm® 


The  volume  loss  operation  subtracts  eadi  resulting  ANij  value  from  the  respective 
array  element.  Under  the  conditions  in  this  effort,  the  three-body  recombination 
volume  loss  has  been  found  to  have  a  very  little  effect  on  the  evolution  of  the  N  atom 
distribution.  A  similar  operation  may  be  applied  for  the  much  faster  reaction  of  NO 
with  atomic  nitrogen  if  oxygen  impurities  exist  within  the  discharge  system. 
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A. 2  Accuracy  of  Model 


Some  degree  of  uncertainty  exists  when  implementing  a  numerical  approach  to  mod¬ 
eling  the  diffusion  and  wall  loss  of  atomic  nitrogen.  The  three-point  method  used 
to  approximate  the  second  derivative  in  the  diffusion  equation  carries  an  uncertainty 
that  scales  as  h^,  where  h  is  the  spacing  between  array  elements.  The  four  point 
method  used  to  approximate  the  first  derivative  boundary  condition  at  a  surface  was 
therefore  chosen  because  its  uncertainty  likewise  scales  as  h^.  One  is  tempted  to 
decrease  h  in  the  program  to  reduce  the  uncertainty,  but  unfortunately  the  speed 
of  the  program  scales  as  well  with  h^.  This  reduction  in  program  speed  is  due  to 
the  combination  of  At  decreasing  as  from  Eq.  A.9  and  the  total  number  of  array 
elements  increasing  as  1/h^.  Since  the  slowing  of  the  program  speed  can  eventually 
become  prohibitive,  a  measure  of  accuracy  should  be  done  with  an  h  value  that  gives 
an  acceptable  speed. 

The  model,  as  constructed,  represents  a  right  circular  cylinder.  A  good  test  of 
the  diffusion  and  wall  loss  modeling  would  be  to  see  if  it  conforms  to  the  exact 
analytic  solution  of  the  fundamental  diffusion  mode  within  a  right  circular  cylinder 
with  partially  reflecting  walls.  To  test  the  accuracy  of  the  numerical  model,  the 
fundamental  diffusion  mode  and  corresponding  decay  rate  is  analytically  determined 
for  a  right  circular  cylinder  with  a  specified  wall  loss  coefficient.  This  result  is  then 
compared  to  the  steady  state  decay  for  similar  conditions  using  the  numerical  model 
with  zero  production  and  no  volume  loss.  Given  an  initial  particle  distribution,  the 
computed  distribution  and  decay  rate  should  converge  in  time  to  the  analytic  results, 
with  any  variance  indicating  an  uncertainty  in  the  model. 

The  anal3rtic  test  case  involves  a  right  circular  cylinder,  defined  with  the  origin 
of  the  coordinates  z  and  p  at  the  center  of  the  cylinder.  The  cylinder  has  a  radius 
of  po  and  a  height  of  Zo,  which  places  wall  surfaces  at  p  =  po  and  z  —  Zo/2.  The 
fundamental  diffusion  mode  distribution  of  a  right  circular  cylinder  with  reflecting 
walls  is  given  by  [77] 

N  {z,  p)  =  N{0, 0)cos{Trz/zs)Jo{2A05p/  Ps)  (A.20) 
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where  Zg  and  axe  parameters  representing  extrapolated  diamber  dimensions  where 
the  particle  density  distribution  would  extrapolate  to  zero  beyond  the  chamber  boimd- 
aries.  These  dimensions  are  determined  by  the  boundary  conditions  and  related  to 
the  fundamental  diffusion  length,  A,  by 


1 

A2 


The  net  particle  loss  frequency  is  then  given  by 


(A.21) 


1  _ 

T~N\dt}~  A2‘ 


(A.22) 


By  combining  Eqs  A.IO  and  A.13,  the  boundary  conditions  can  be  written  as 


1  /  dN\  V  7 

N\d^)b^~AD{l-iy 


(A.23) 


The  extrapolation  parameters,  z,  and  p,,  can  now  be  determined  for  a  given  wall  loss 
coefficient,  7.  For  the  z  dimension,  using  the  derivative  of  the  fundamental  mode 
distribution  from  Eq.  A.20,  the  boundary  condition  is 


1 

,p)\dz)  b=Zo/2  Zs  2Zs 

and  for  the  p  dimension. 


(A.24) 


Y 


dN\ 


2.405  Ji(2.405po/p,) 

Ps  •^o(2.405po/pa) 


(A.25) 


N(z,po)\dpJb=:po 

Combining  Eq.  A.23  with  Eqs.  A.24  and  A.25  results  in  a  set  of  two  transcendental 
equations 

V  •y  'IT  7rr_ 

(A.26) 


4D  (1  -  a  -  2^;) 


and 


V  7  _  2.405  Ji(2.405po/p«) 


Ps  Jo(2.405po/p,) 

which  can  be  solved  for  a  given  7  to  determine  z^  and  p^.  Subsequently,  these  values 
can  be  inserted  into  Eqs  A.21  and  A.22  to  determine  the  decay  rate  of  the  fundamental 
diffusion  mode. 

The  analytical  and  numerical  models  are  solved  for  a  cylinder  of  height  Zo=  2.2 
cm  and  radius  po  =  5  cm  with  all  surfaces  having  a  loss  coefficient  of  7  =  1%.  The 
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diffusion  constant  for  N  atoms  through  3  Torr  of  N2  is  D  =  73  cm2/s[71]  and  assuming 
T  =  300  K,  a  velocity  of  v  —  73,000  cm/s  is  used. 

For  the  numerical  model,  the  production  of  particles  was  set  to  zero  as  was  the 
rate  of  volume  loss.  The  grid  element  spacing  was  set  to  h  =  1  mm,  at  which  the 
program  speed  is  sufficiently  fast.  The  value  of  At  corresponding  to  this  grid  spacing 
is  34.1  //s.  The  approach  taken  to  test  the  numerical  program  was  to  begin  with 
some  arbitrary  initial  N  atom  distribution.  After  running  for  a  period  of  time,  the 
calculated  distribution  and  decay  rate,  1/r,  should  reach  a  steady  state  mode  that  can 
be  compared  to  the  analytic  solution.  The  decay  rate  of  an  arbitrary  array  element 
was  monitored  as  the  program  ran  in  order  to  determine  the  degree  of  steady  state 
reached.  The  final  distribution  was  then  fitted  to  Ekj.  A.20  to  determine  the  z,  and 
Pa  of  the  numerical  solution.  Regardless  of  the  initial  distribution,  the  calculated  N 
atom  distribution  was  found  to  reach  a  the  same  steady  state  condition. 

The  analytical  expressions  were  solved  directly  by  insertion  of  the  given  parameter 
values.  The  results  of  both  the  analytical  and  numerical  methods  are  shown  in  Table 
12. 


2^8 

Ps 

l/r 

Analytical 

2.9561  cm 

5.4116  cm 

97.307  S-" 

Numerical 

2.9562  cm 

5.4081  cm 

97.368  s-^ 

Table  12:  Numerical  and  Analytical  Model  Comparison  for  7  =  1%. 


In  the  comparison,  each  parameter  is  found  to  have  less  than  a  0.1%  error  com¬ 
pared  to  the  exact  analytic  solution.  This  example  proves  that  the  model  with  h=l 
mm  reaches  the  steady  state  solution  with  good  accuracy. 
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